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THE SPECTRUM OF NEUTRONS FROM THE IRRADIATION 
OF URANIUM BY 14-MEV NEUTRONS! 


R. L. CLARKE 


ABSTRACT 


The spectrum and absolute yield of neutrons resulting from 14-Mev neutrons 
incident on natural uranium was observed at 90°+16 to the incident neutron 
beam by a time-of-flight method. Further observations of the spectrum from 
9 to 13.5 Mev were made for scattering angles at 30° and 60°. In addition to the 
expected neutrons resulting from evaporation and fission, strong evidence is seen 
for direct interaction inelastic scattering. The (,m’) cross section is estimated 
to be 150+70 millibarns. 


INTRODUCTION 


The spectrum of neutrons resulting from 14-Mev neutrons incident on 
natural uranium has been given by Zamyatnin ef al. (1957) and by Vasilev 
et al. (1957). The experiment of Zamyatnin was done with photographic plates 
and covered the region from 0.5 to 6.0 Mev. The work of Vasilev was done 
with a time-of-flight system and covered the region from 0.05 to about 1.1 
Mev. These workers did not give any normalization of their yields either on 
an absolute scale or relative to each other.* 

The object of the present experiment was to extend the spectrum measure- 
ments to a higher energy and to give the yields on an absolute scale. The 
shape and angular distribution of the high-energy part of the spectrum strongly 
suggest the presence of direct interactions. The shape and intensity of the 
lower-energy part of the spectrum are expected to depend sensitively on the 
competition between neutron evaporation and fission following compound 
nucleus formation, which has been discussed by, for example, Jackson (1956) 
and Leachman (1958). Hanna and Clarke (1961) have made numerical 
estimates of the contributions from fission and evaporation and have com- 
pared these with the experimental results to give some information about 
the nuclear processes involved. 

1Manuscript received January 16, 1961. 


Contribution from Atomic Energy of Canada Limited, Chalk River, Ontario. 

Issued as A.E.C.L. No. 1240. 

*NOTE ADDED IN PROOF: Recently Vasilev et al. (Vasilev et al. Soviet Physics J.E.T.P. 11, 
483 (1960)) have measured the spectrum of neutrons emitted in coincidence with fission. 
Their data suggest a higher cross section in the region of 5 Mev than was observed in the 
present work. The reason for the discrepancy is not clear. 


Can. J. Phys. Vol. 39 (1961) 
957 








958 CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


APPARATUS 

The spectrum of neutrons resulting from 14-Mev neutrons incident on 
natural uranium was measured by a time-of-flight technique, using neutrons 
from the TD reaction. These neutrons were produced by a 120-kev deuteron 
beam which struck a tritium-zirconium target. Observations were carried 
out at right angles to the deuteron beam; hence the energy of the incident 
neutrons was 14.1 Mev. 

A schematic view of the apparatus, as seen looking toward the accelerator, 
is shown in Fig. 1. The recoil alpha particles from the reaction are counted 
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Fic. 1. The time-of-flight apparatus. The deuteron beam is normal to the plane of the 
figure and strikes the T—Zr target just above the a-particle counter. The shaded areas represent 
plastic scintillators. 


by means of a thin plastic scintillator in direct contact with an R.C.A. 6342 
photomultiplier. The scattered neutrons were detected in a plastic scintillator, 
7 in. diameter by 5 in. thick, coupled to a 5 in. diameter Dumont photo- 
multiplier. The measured times, from which the neutron flight times were 
deduced, were taken from the occurrence of a pulse in the alpha counter to 
the occurrence of a pulse in the neutron counter. This delayed coincidence 
requirement between the neutron and alpha counters in effect defined a beam 
of neutrons whose width was determined by the collimator in front at the 
alpha counter. The full width of this beam at half maximum was 32°. The 
center of the phosphor was 200 centimeters from the target on a line approxi- 
mately at right angles to the deuteron beam axis. 

The sample consisted of a sheet of natural uranium 10 cm X10 cm X0.7 cm 
supported on a light brass frame at an angle of 45° to the neutron beam, 12 
centimeters from the target. The major portion of the neutron beam defined 
by the alpha-particle collimator thus passed through the sample. The thick- 
ness of the sample was so chosen that about 28% of the neutrons passing 
through it suffer an interaction. Of these interactions, half are due to elastic 
scattering which, because it is strongly peaked in the forward direction, has 
little effect on the experiment. 

A copper shadow bar was placed between the target and the neutron counter 
to reduce the background due to direct neutrons from the target. Additional 
shielding was placed between the neutron counter and those portions of the 
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accelerator, such as diaphragms in the beam path, which gave rise to signi- 
ficant background. With this system only a few counts per minute were 
observed, distributed over all flight times. Thus it was necessary to count 
for long periods to get adequate statistics. 

This system was used for the high-energy portion of the work. 

The time-to-pulse-height converter resembled, with only minor modi- 
fications, that of Neilson and James (1958). The distribution of flight times 
was displayed on a 30-channel pulse-height analyzer on which each channel 
corresponded to about 3.510 seconds. The time-to-pulse-height converter 
was calibrated in three ways: by the insertion of known delays in artificial 
pulses put in at the preamplifiers, by changing the length of the flight path 
with the neutron beam aimed directly at the neutron counter, and by observ- 
ing the flight time of neutrons scattered by hydrogen through known angles. 
The satisfactory agreement between these methods showed that the equip- 
ment was operating as required. The resolution of the system was found to 
be ~10X10~ seconds. This was limited principally by the 5 in. photomulti- 
plier. 

The energy sensitivity of the neutron detector was measured in two ways. 
In the first a cylindrical piece of polyethylene 2.5 cm in diameter and 5 cm 
long was supported 25 centimeters from the target with its axis parallel to 
the deuteron beam axis. The neutrons scattered from both the hydrogen and 
carbon in the sample can be seen as peaks in the time spectrum in Fig. 2. 


1000 
800 
600}- 


400 


YIELD ARBITRARY UNITS 


200 





° 3 8 2 16 20 24 28 32 36 
CHANNEL NUMBER 


Fic. 2. The spectrum of neutrons from 14-Mev neutrons scattered by polyethylene. The 
scattering angle is 40°. The peaks due to elastic scattering from carbon, inelastic (Q = —4.43) 
scattering from carbon, and scattering from hydrogen are labelled I, II, and III, respectively. 
The areas under these peaks, less the background, were used to calibrate the detector. Part 
of the y-ray peak can be seen at the right of the figure. 


The observed yields, that is the areas under the peaks, were used to calibrate 
the energy sensitivity of the detector in terms of the known scattering cross 
sections for hydrogen (Hughes and Schwartz 1958) and carbon (Anderson 
et al. 1958). The results from hydrogen as the angle @ (Fig. 1) was varied from 
25° to 55° gave a calibration from 3.9 to 9.7 Mev. The elastic scattering from 
carbon gave a point at 13.5 Mev. 
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In the second method the polyethylene scatterer was replaced by a cylinder 
of plastic phosphor 2.5 cm in diameter and 5 cm long, connected to an R.C.A. 
6342 photomultiplier. This scattering counter was used instead of the alpha- 
particle counter to give the zero time signal in the time-of-flight system. With 
a bias on the scattering counter set low enough, all events due to neutrons 
scattered from the hydrogen in it to the neutron counter could be identified 
from their flight times. Because of the relatively low counting rate in the 
scattering counter and because the neutrons scattered from carbon did not 
give rise to any background, this method of calibration was considered more 
accurate than the first. It gave a calibration from 2.2 to 12.0 Mev. 

In both methods the scattered neutron energy was calculated from the 
geometry rather than from the position of the neutron peak in the time 
spectrum. The results from the two methods agree within the limits of experi- 
mental error. The energy sensitivity curve finally used is shown in Fig. 3(a). 
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Fic. 3. The sensitivities of the neutron counters to neutrons of various energies, (a) for 
a plastic scintillator 17.8cm diameter X 12.7 cm thick, (6) for a stilbene crystal 2cmX2 cmX 
4 cm. In Fig. 2(a) the crosses refer to data taken by scattering from a scintillator (see text), 
the triangles by scattering from polyethylene, the circles to an average of an earlier cali- 
bration using both the above techniques. The units of sensitivity used are the same for both 
counters. 
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This counter could not be used for the lowest-energy part of the spectrum 
because of the rapid increase of background as the bias was reduced below 
3 Mev. For this part of the spectrum a neutron counter consisting of a stilbene 
crystal 2 cm X2 cm X4 cm in contact with an R.C.A. 6342 photomultiplier was 
used with a flight path of 1 meter. This counter was calibrated in the same 
manner as the large counter and was found to have a resolution of 6X10~° 
seconds. The sensitivity curve is shown in Fig. 3(0). 

Both these calibrations gave relative efficiencies as a function of energy; 
the absolute value of the uranium neutron yield was obtained in terms of the 
known carbon cross section by comparison with the yield from a carbon sample 
of known mass and the same geometry as the uranium. 

The angular distribution of high-energy neutrons is important in deter- 
mining the processes which govern neutron emission, as discussed below. For 
those neutrons of energy greater than 9 Mev the angular distribution was 
measured using a later version of the instrument. A 10 cm diameter by 10 cm 
long cylinder of plastic phosphor was mounted on an R.C.A. 7046 photo- 
multiplier and surrounded by heavy shielding. This system gave a resolution 
of 2X10-* seconds. In this case the flight path was 2.76 meters. 

It was found that there was a substantial dependence of the background 
in the neutron counter on the beam intensity. Hence, the random coincidence 
background was approximately proportional to the square of the beam current. 
In spite of this, an accurate subtraction of the accidental coincidence rate 
could be made by assuming that this background was proportional to the 
product of the total alpha-particle and neutron counts. This was possible 
because during the rather long periods of observation required, the beam 
current was automatically controlled so as to hold the alpha-particle counting 
rate, and hence the neutron production constant. The background which 
was used for each channel of the time spectrum was found from a run made 
with the sample removed, but with other conditions unaltered. The back- 
ground spectrum thus found was weighted by the product of the neutron and 
alpha-particle counting rates and subtracted from the spectrum observed with 
the sample in place. 

RESULTS 

A typical time-of-flight spectrum after background subtraction but before 
correction for counter sensitivity is shown in Fig. 4. This is a composite spec- 
trum, the result of several runs taken with the large phosphor. The peaks 
due to gamma rays, elastically scattered neutrons, and inelastically scattered 
neutrons can all be seen. The uncertainty in the shape of the elastic peak 
introduced an error of a factor of about 2 at 12 Mev. This error becomes 
negligible at 10 Mev. In subtracting the elastic peak from the observed 
spectrum, the peak was assumed to be symmetrical about its maximum. This 
procedure may overestimate slightly the yield of inelastic neutrons. 

No correction to the spectrum shape was made for multiple scattering in 
the sample. Over the range of neutron energies of interest, from 3 to 14 Mev, 
both the total and non-elastic cross sections show only a small variation 
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Fic. 4. Time-of-flight spectrum observed at 90° for 14-Mev neutrons scattered from natural 
uranium using the apparatus of Fig. 1. Peak (1) corresponds to y-rays, peak (2) to elastically 
scattered neutrons, and the continuum (3) to the inelastically scattered neutrons. 


(Hughes and Schwartz 1958; and Almqvist and Clarke 1955). Since the effect 
of non-elastic scattering is in general to degrade the neutron energy to less 
than 2 Mev, such a scattering event will remove a neutron entirely from 
observation by the present experiment. Since elastic scattering is predomi- 
nantly forward, its effect is to spread out slightly the angular resolution of the 
experiment. This effect will be most pronounced at high energies, where, as 
is shown in the present experiment, the angular distribution is most anisotropic. 
At energies below 8 Mev, where the distribution is isotropic, there should be 
no distortion of the spectrum due to elastic scattering. A rough estimate 
shows that at the highest observed energy, 12 Mev, the distortion is less than 
2% so that to the accuracy of the present experiment, it may be neglected. 

Because the uranium sample was about four times thicker than the carbon 
sample in terms of neutron mean free path, a correction was made for self- 
absorption in making the comparison between carbon and uranium absolute 
yields. The result was to increase the uranium cross section by 10+3%. 

The energy spectrum obtained from the time-of-flight spectrum after the 
corrections discussed above have been applied, is shown in Fig. 5. The figure 
is a combination of the data from the three detectors. The data of Zamyatnin 
(1957) and Vasilev et al. (1957) are shown for comparison. Their data have 
been normalized to the present work only by eye. Although none of the three 
curves shown overlap each other perfectly, there is general agreement as to 
the shape of the spectrum. The most conspicuous features are the low energy 
at which the maximum occurs and the flattening of the curve at energies 
above 6 Mev. 
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Fic. 5. The spectrum of neutrons observed at 90° from the scattering of 14-Mev neutrons 
from natural uranium. Curve | is the data of Vasilev et al. (1957), curve II that of Zamyatnin 
et al. (1957). The experimental points and curve III are the data of the present work. The 
curve is drawn only to aid the eye. The errors shown are due to counting statistics only, as 
discussed in the text. 


Sources of systematic error may be summarized as follows: 

(1) In the calibration of the time-to-pulse-height converter, there is an 
uncertainty of about 5% in the number of nanoseconds per channel and an 
uncertainty of about 0.5 channels in the position of the channel corresponding 
to the origin of the time scale. These uncertainties produce errors both in the 
energy scale and in the cross-section values. The latter arises in the calculation 
of the yield per Mev from the observed yield and the calculated channel 
width in Mev. The error thus introduced in the cross section is a maximum 
near 5 Mev, where it reaches a value of about 20%. 

(2) For energies near 14 Mev there is an error due to the uncertainty in 
subtracting the elastic peak. This introduces an uncertainty in the last few 
points, rising from zero for the point near 9 Mev to a factor of 2 for the point 
at 12 Mev. 

(3) The error in the normalization of the uranium yield relative to the carbon 
cross section of Anderson et al. (1958) arises from the self-absorption cor- 
rections discussed above, from statistics, from uncertainties in placing the 
carbon in the same position as the uranium, and from the limited stability 
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of the equipment over the considerable periods of time involved. This error 
may be taken to be 5%. The error in the work of Anderson e¢ a/., which appears 
to be less than 5%, is negligible. 

(4) The spectrum and the corrections to it become increasingly unreliable 
at low energies due to uncertainties in the sensitivity curve. Allowance for 
this has been made in assigning errors at the lower end of the spectrum. 

Because the yield of neutrons for energies greater than 6 Mev is difficult to 
account for on the basis of either evaporation or fission (Hanna and Clarke 
1961), it was felt that it might be due to direct interaction. Spectra were 
measured at scattering angles of 30°, 60°, and 90° using the 10 cm X10 cm 
detector. The results after subtraction of the elastic peak are shown in Fig. 6, 
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Fic. 6. The yield of neutrons resulting from 14-Mev neutrons incident on natural uranium. 
Symbols refer to angles of observation; circles, 90°; stars, 60°; triangles, 30°. In each case the 
contribution of elastically scattered neutrons has been subtracted. The curves are drawn 
only to aid the eye. 


where it can be seen that the yield is higher at small angles, the effect being 
most pronounced as the energy approaches 14 Mev. These data are quali- 
tatively similar to the data of Rosen and Stewart (1955) on bismuth and 
suggest that the forward peak in the neutron yield is a result of direct inter- 
action processes. 

Using these data, an estimate of the (,n’) cross section can be made. The 
thresholds for the (y,7) and (y,a) reactions on U*8 are 6.0 and 5.8 Mev, re- 
spectively (Cameron 1956). Thus, it is reasonable to assume that those neutrons 
which, on inelastic scattering, leave the target nucleus, U*’, excited to less 
than 6 Mev are the outgoing neutrons from the (,n’) reaction. In estimating 
the number of these neutrons from the angular distribution data, the total 
yield from 8 Mev to 14 Mev is required for all angles. The present data were 
extrapolated by assuming (1) that differential cross section is constant from 
90° to 180° angle of observation, and (2) that the cross section from 8 Mev to 
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9 Mev is the same as that at 9 Mev. The former assumption probably over- 
estimates the cross section; the latter probably constitutes an underestimate. 
In integrating over energy to get the total o(m,n’) the region from 13.5 Mev 
to 14 Mev was neglected. There is really no way in this experiment of dis- 
tinguishing scattered neutrons which have lost only 0.5 Mev from elastic 
neutrons, nor is there any way of estimating their total yield. Subject to these 
uncertainties, a value of o(u,n’) or 150+70 millibarns is obtained. 


DISCUSSION 

The shape of the observed spectrum, with its maximum at 0.5 Mev, and 
the relatively large yield at energies up to 8 Mev strongly suggest that both 
evaporation and fission neutrons contribute significantly to the spectrum. An 
analysis of the shape in terms of multiple neutron evaporation from U** 
followed in some cases by fission has been given by Hanna and Clarke (1961). 
From the ratio of the total area under the curve, 7.8 barns and the non- 
elastic cross section, 2.90 barns, a value of 2.7 is obtained for 7, which is in 
satisfactory agreement with experiment (see Hanna and Clarke). Because of 
the difficulty in normalizing the data of Vasilev et al. (1957) and Zamyatnin 
et al. (1957) to the data of this experiment, the value given here should not be 
considered an accurate determination. However, it does indicate that the cross 
sections found are consistent with other experimental determinations of 7. 

In Fig. 6 it can be seen that between 11 and 14 Mev an increase in yield is 
seen, especially at small angles. This is not seen in Fig. 5 because of the large 
uncertainties in the data near 14 Mev. This increase in differential cross 
section resembles the increase in cross section seen at small angles for indi- 
vidual levels in light nuclei (Anderson et al. 1958; Clarke and Cross 1959, 
1960). For those neutrons from U having energies greater than 12.5 Mev, 
corresponding to levels in U of low excitation, the ratio of the yield at 30° to 
that at 90° is of the order of 10:1, which is greater than the ratio of yields seen 
for the first excited states of S*® (W = 2.24 Mev) and Mg*4(Q = 1.37 Mev). 
The angular distributions seen for neutrons in the region 9 Mev to 12 Mev, 
corresponding to Q = —2 to —5 Mev, appear to be comparable to those seen 
in Sand Mg for individual states having similar Q values. Thus it is reasonable 
to suppose that the angular distribution seen in U scattering arises from causes 
similar to the direct interactions seen in light nuclei. This hypothesis is further 
supported by the large intrinsic quadrupole moment of U**, which would 
be expected to give rise to strong surface direct interaction according to the 
models of Inopin (1956), Ui (1959), and Blair (1959). 

As stated above, the total cross section for neutrons of energy between 8 
and 13.5 Mev is about 150 millibarns, which is 15% of the total non-elastic 
cross section. The portion of the cross section corresponding to excitations 
less than 0.5 Mev was not measured, and the fraction of this cross section 
which may be due to direct interaction is not known, consequently the true 
direct interaction cross section is not known. However, the ratio of this high- 
energy yield is similar to the ratios of direct to total non-elastic seen for Mg, 
C, and Be, which lie between 10% and 30% (Anderson et al. 1958; and Ball 
et al. 1958). 
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NEUTRON EVAPORATION IN THE 14-MEV 
NEUTRON FISSION OF URANIUM! 


G. C. HANNA AND R. L. CLARKE 


ABSTRACT 


The numbers of prefission and postfission neutrons expected from the 14-Mev 
neutron fission of uranium have been calculated from the experimental data on 
partial cross sections and neutron multiplicities. It has been assumed that the 
competition between neutron emission and fission does not change with increasing 
energy beyond threshold. The neutron energy spectrum measured recently by 
Clarke gives a smaller intensity in the energy range 4-6 Mev than is expected 
from the calculated number of postfission neutrons. The measured intensity 
suggests that only 2.1+0.5 neutrons are evaporated from the moving fission 
fragments. 


INTRODUCTION 

As described in the preceding paper, Clarke (1961) has measured the energy 
spectrum of neutrons resulting from the interaction of 14-Mev neutrons and 
natural uranium. He has also measured the angular distribution of the highest- 
energy neutrons and concluded that these are probably the result of direct 
interaction processes. Below about 6 Mev the direct interaction neutrons can 
be neglected in comparison with those emitted following compound-nucleus 
formation. These arise from two sources. The lowest-energy region is domi- 
nated by neutrons evaporated from the compound nucleus either before 
fission takes place or when there is no fission. At somewhat higher energies, 
especially between 4 and 6 Mev, most of the observed intensity is due to 
neutrons evaporated from the moving fission fragments. 

Enough is known of the course of the reactions to predict with reasonable 
accuracy the total number of prefission and postfission neutrons, and, for the 
latter neutrons, the energy spectrum as well. The present paper calculates 
these expected effects and compares them with what is observed experi- 
mentally. 

EXPECTED INTENSITIES 

The capture of a 14.1-Mev neutron in U** forms U*® with an excitation 
energy of 18.86 Mev (Foreman and Seaborg 1958), which may undergo 
fission or evaporate a neutron (its de-excitation by gamma radiation is very 
improbable). In the latter case a U** nucleus is formed with an excitation 
energy that is less than that of the parent U* by the sum of the neutron 
kinetic energy, which is on the average 1 Mev,* and the neutron binding 
energy; the average excitation of the U®®* is therefore 13.1 Mev. This U?* 


1Manuscript received January 17, 1961. 
Contribution from the Nuclear Physics Branch, Atomic Energy of Canada Limited, Chalk 


River, Ontario. 
Issued as A.E.C.L. No. 1241. 
*That is, twice the temperature of the evaporation spectrum, which is close to 0.56 Mev 


(see below). 
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nucleus can in turn undergo fission or evaporate another neutron, and so on 
according to the scheme of Fig. 1. The competition between the different 


ia,t)} tee gO ts) 
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ni 
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Fic. 1. Nuclear reactions produced by 14-Mev neutrons on U2*, 


processes is expressed in terms of the various partial cross sections, on which 
there is fairly complete experimental information (see Appendix). The finally 
adopted values are given in Table I, which also lists N;, the number of neutrons 


TABLE I 


Partial cross sections and neutron multiplicities in the interaction of 
14.1-Mev neutrons with U2% 


elas ) ; “wll o( ss sai o(n,xn) 


Excitation o(n,xnf) 
x Nucleus (Mev) (barns) Nt x Ne: X(Ni+x) (barns) <x 
0 U-239 18.86 0.55 4.22, 2.32; 2.323 0 0 
1 U-238 13.1 0.45 3.407 1.533 1.98; 0.15 0.15 
2 U-237 6.0 0.17 2.40 0.403 0.745 0.84 1.68 
3 U-236 ~ 0 0 0 0 0 0.74 2.22 
Totals Lag 4.26, 5.054 1.73 4.05 

Note: » = number of postfission neutrons per fission = 4.264/1.17 = 3.644. 

» = number of pre- and post-fission neutrons per fission = 5.054/1.17 = 4.320. 

n = number of secondary neutrons per interaction = (5.054 + 4.05)/2.90 = 3.14. 


emitted following the different types of fission; the basis of these estimates 
is given in the Appendix. 

Table I shows that the average number of postfission neutrons per fission, 
u, is 3.64 compared with the average total number per fission (v) of 4.32, 
which includes the prefission evaporation neutrons. The average number of 
neutrons produced per interaction, 7, is calculated to be 3.14, which may be 
compared with recently reported direct measurements of 2.80.25 by Lebedev 
et al. (1958) and 2.99+0.15 by Flerov and Talyzin (1958) and with an un- 
published value of 3.29+0.15 by Allen (1960). The value obtained from the 
curve of Fig. 2 is 2.7, but this result is sensitive to the fitting of the low-energy 
Russian data to the rest of the curve, and Clarke’s (1961) experiment cannot 
be regarded as a satisfactory measurement of 7. 
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It is therefore not legitimate to discuss his experimental data in terms of 
u/n, although this would have the advantage of eliminating uncertainties in 
absolute intensity measurements. Rather, the calculated value of u is to be 
compared with the value expected from the observed intensity in the 4- to 
6-Mev region. It should therefore be noted that yu is calculated from the 
published values of the directly measured quantity v by the relation 


on on watt 20n ,2nt 
SF 
Ot 
which is independent of the non-fission cross sections. Estimates of the un- 
certainties in the individual terms are given in the Appendix; they lead to 
a value of about +0.18 for the uncertainty in uw. No precise meaning can be 
attached to this estimate of the uncertainty, but, at +5%, it is clearly smaller 
than that of the experimental intensity measurements. 


THE EXPECTED FISSION-NEUTRON SPECTRUM 

The total area of the postfission-neutron spectrum, yo;, is expected to be 
4.26 neutron barns (Table 1). With the assumption that all these neutrons are 
evaporated from the moving fission fragments, the spectral shape is expected 
to be rather similar to that of the slow neutron fission of U2*.* Wahl (1954) 
has shown that the fragment velocities in 14-Mev fission are essentially the 
same as in thermal fission. However, the 14-Mev fission-neutron spectrum 
is expected to be slightly more energetic because the fragment temperatures 
are higher. The spectrum measured in a particular direction will be further 
modified as a result of the angular correlation of the fission neutrons and 
fission fragments since the fragments are not isotropic with respect to the 
14-Mev neutron direction. 

The correction for fragment temperature was calculated following Terrell 
(1959), who showed that the fission-neutron energy spectrum can be adequately 
represented by a Maxwellian distribution (2/./x) (E/é)! exp(—E/#) (dE/®), 
where @ is a fictitious temperature defined by 30/2 = E, the mean neutron 
energy in the laboratory system. From a study of the available data Terrell 
has concluded that FE = 0.78+0.621 (N;+1)? Mev, where N; is the number 
of postfission neutrons (and is equal to v in systems where only (,f) reactions 
are possible). Terrell’s scale of »v values is based on 2.46 for U?® thermal fission ; 
for the present analysis, with v equal to 2.43 (Hughes et a/. 1960), his relation 
becomes E = 0.78+0.624 (N;+1)!. The postfission-neutron spectrum expected 
from the 14-Mev neutron fission of U*** is therefore the sum of three compo- 
nents with (see Table I) intensities of 2.323, 1.533, and 0.408 neutron barns 
and fictitious temperatures corresponding to, respectively, N; values of 4.224, 
3.407, and 2.40, i.e. 6 values of 1.470, 1.393, and 1.287 Mev. 

Table II gives the ratio, as a function of energy, of the intensity of this 
complex spectrum to that of the simple spectrum with @ equal to 1.287 Mev. 
These will be referred to as the “hard’’ and “soft’’ spectrum respectively. 


*This is true provided that the angular distribution in the fragment system, generally 
assumed isotropic, is the same in both cases. 
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TABLE II 


The ratio of the intensity of the expected postfission-neutron spectrum in 14-Mev fission 
to that of the thermal fission ageeTnaG. 


E( Mev) 0.5 1 1.5 2 3 4 5 6 8 10 


Ratio 
hard/soft 0.901, 0.925, _0.960 0.996 1.073 1.157 1.247 1.347 1.59 1.84 


It is difficult to estimate the uncertainty in this hardening correction, which 
is appreciable at high neutron energy. Fortunately its reliability is not of 
crucial importance for the purposes of the present discussion. 

In the 14-Mev neutron fission of U** the fragments are emitted preferentially 
in the direction of the incoming beam; according to Henkel and Brolley 
(1956) the ratio of intensities at 0° and 90° is 1.40+0.14. In the present experi- 
ment the secondary neutrons are observed at 90° to the 14-Mev neutron 
beam axis and thus preferentially at 90° to the direction of motion of the 
fission fragments. The neutrons evaporated from the moving fragments are 
peaked in the fragment direction, and the present experiment therefore 
observes fewer fission neutrons than would be expected if the fission process 
were isotropic. The magnitude of this effect was calculated by using a simpli- 
fied model of the fission process due to Milton (1958), which assumes a single 
typical mass splitting and, from each fragment, an equal neutron evaporation 
characterized by two nuclear temperatures, so that in the fragment system 
the energy spectrum from each fragment is 


(—s), (en fs) 
ne 1 T; =? Ts 

a, 7, and 7», are adjusted to fit the energy spectrum observed in the laboratory 

system. Although this model is evidently oversimplified it is acceptable for 


the present purpose since the calculated effect is quite small (see Table III). 





N(n) = = 2 EXp —— 


TABLE III 


The effect of angular correlation on the fission neutron flux at right angles to the 14-Mev 
neutron beam 


Laboratory energy (Mev) 0.025 1 2 4 6 8 10 
(Average intensity —intensity at 90°) ae as ; re Soni 
Average intensity X100 0 2.3 3.5 4.2 48 6.2 7.0 7.3 








Norte: These numbers were caine for a 14-Mev neutron beam of zero angular divergence; for the ae 
(~ 30°) of the present work they should be multiplied by about 0.9. 


COMPARISON WITH THE EXPERIMENTAL DATA 


The ey data of the preceding paper aan 1961) are reproduced 
in Fig. 2. The smoothed data, curve III of Fig. 2, are reproduced as curve 
A of a 3, which also shows two calculated fission spectra, B and C, both 
corrected for angular distribution effects. Curve B is the ‘‘hard”’ spectrum 
corresponding to » = 3.64 postfission neutrons per fission. The total neutron 
emission is expected to be the sum of this spectrum, a lower-energy evaporation 





HANNA AND CLARKE: NEUTRON EVAPORATION 971 


1000 


100 


mb 
Mev x STER) 
3 


*( 





oO 2 4 6 8 10 12 14 
E (MeV) 


Fic. 2. The spectrum of neutrons observed at 90° from the scattering of 14-Mev neutrons 
from natural uranium. Curve I represents the data of Vasilev et al. (1957), curve II those of 
Zamyatnin et al. (1957). The experimental points are the data of Clarke (1961) and the errors 
hon are due to counting statistics only. Curve III is drawn through these data only to aid 
the eye. 


spectrum, and a weak higher-energy component resulting from direct inter- 
action processes. However, it appears that in the neighborhood of 5 Mev the 
calculated intensity of the fission spectrum alone is greater than the total 
observed intensity. 

This suggests that not all the postfission neutrons are evaporated from the 
moving fragments. Curve C has been calculated on the assumptions that the 
moving fragments are at no higher temperature than in thermal-neutron- 
induced fission, and that the ‘“‘excess’’ postfission-neutron emission occurs 
between the saddle point (where the nucleus is committed to fission) and the 
scission point (where it actually ruptures). Curve C has a total intensity of 
2.81 neutron barns (2.40 neutrons per fission) and a correspondingly ‘‘soft”’ 
shape. 

If curve C is subtracted from A the difference spectrum in the region below 
about 5 Mev has roughly the shape expected for evaporation from a nucleus 
having a single temperature 7, i.e. N(E) « Ee~®/". Curve D is N(E)/E, and 
it will be seen that it is approximately exponential in shape with a temperature 
of 0.55+0.07 Mev. The trend of curve D in the region above 4 Mev suggests 
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Fic. 3. A comparison of the smoothed experimental data (A) with fission-neutron spectra 
calculated on the assumption that all postfission neutrons are evaporated from the moving 
fragments (B) and that only 2.4 neutrons per fission are so evaporated (C). Curve D is the 
difference (A—C) divided by the energy (see text). 


that too large a postfission contribution has been subtracted. However, an 
increase of about 15% in curve A in this region, which is within the experi- 
mental uncertainty, would remove the downward trend in D. 

This value of 7 agrees with that of 0.48+0.05 Mev given by Zamyatnin 
et al. (1958) for U8 and 14-Mev neutrons, but it is substantially lower than 
the temperatures of 0.75 to 1 Mev found by Graves and Rosen (1953) in the 
region of lead and bismuth. The difference is greater than expected from the 
combined effects of mass number on the level-density parameter and of the 
greater importance of the (”,3”) reaction in uranium. Although it is not 
reasonable to expect that the complex evaporation phenomena occurring here 
can be accurately represented by a single nuclear temperature, a more elaborate 
analysis has not been tried because of the sensitivity of the upper part of the 
evaporation spectrum to the absolute amplitude of curve A and the assump- 
tions made about the fission-neutron intensity. 

To return to the fission neutrons, Batzel (1954) was apparently the first 
to mention the possibility of a substantial neutron emission between the 





HANNA AND CLARKE: NEUTRON EVAPORATION 973 


saddle point and scission. It may also be recalled that the angular distribution 
measurement of Harding and Farley (1956) implies that, out of the 13 neutrons 
emitted in the 147-Mev proton fission of uranium, the number emitted by the 
moving fragments (u’) is only 2.5+1. In the present experiment y’ is deter- 
mined essentially by the absolute intensity measurement in the 4- to 6-Mev 
region, which is accurate to about +20%. Accordingly, uw’ = 2.4—6+0.5, 
where 6 represents the correction for the tail of the evaporation spectrum; 
extrapolation with JT = 0.55 Mev gives 6 = 0.33, ie. w’ = 2.140.5. 

It should, however, be pointed out that there is another possible inter- 
pretation of the data. An alternative explanation of the observed low intensity 
near 5 Mev is that the neutrons are evaporated predominantly at right angles 
to the direction of motion of the fission fragments. Terrell (1959) has shown 
that an anisotropy of (1—0.4 cos’@) reduces the intensity at 5 Mev by about 
5%. Thus an exceptionally large anisotropy would be required to explain the 
observed effect. 
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APPENDIX 
The Partial Cross Sections of U?*® for 14.1-Mev Neutrons 

The total non-elastic cross section at 14 Mev has been measured as 3.00+ 
0.04 barns by Almqvist and Clarke (1955) and by Lebedev ef al. (1958) as 
2.91+0.14 barns. Beyster, Walt, and Schrandt (1956) give 2.89 barns, and an 
extrapolation of the data of MacGregor, Ball, and Booth (1957) gives 2 2.8+0.1. 
A value of 2.90 barns will be used here. 

The fission cross section of U*** has been measured by Smith, Henkel, and 
Nobles (1957) as a function of neutron energy; at 14.1 Mev the total fission 
cross section is 1.17+0.05 barns.* Their results show that the competition 
between fission and neutron emission, which have similar ‘‘thresholds”’ in 
all the nuclei involved, is not very sensitive to an increase in excitation energy 
beyond the threshold. It is therefore possible to extrapolate with some con- 
fidence from lower energies to obtain values of 0.55+0.1 barns for o(n,f), 
0.45+0.1 barns for o(u,nf), and 0.17+0.05 barns for o(m,2nf) at 14.1-Mev 
neutron energy. 

Of the non-fission cross sections, o(n,y) is negligible (Hughes and Schwartz 
1958), o(n,n’) has been measured by Clarke (1961) as 0.15+0.05 barns,t and 
a(n,2n) as 0.84+~0.15 barns at 14.1 Mev by Graves, Conner, Ford, and 
Warren (1958). o(n,3”) may then be obtained by difference, as gne—o;—a(n,n’) 
—a(n,2n) = 0.74+~0.16 barns. 


The Number of Neutrons Emitted in Fission 

The number of neutrons evaporated from the different uranium isotopes is 
given directly by the values of the partial cross sections. The number of 
neutrons following fission can be calculated from the experimental data on 
v for neutrons of different energies. 

Probably the best values for vy in the thermal neutron fission of U* and 
U8 are 2.43 and 2.50 respectively (Hughes et al. 1960). After subtracting the 
delayed neutron contributions and correcting for changes in excitation energy 
by taking dv/dE = 0.14 Mev (see below), values of 2.359 and 2.375 for 
U** and U** at 6-Mev excitation are obtained. There are a large number of 
measurements for v for the fission of U**® by neutrons of a few Mev energy 
(see Leachman 1958) which give a value of 2.503 for U®® with an excitation 
energy of 6 Mev. A value of vy appropriate to the fission of U**’ at this excitation 
energy has been obtained by weighting the three values for U**, U?%*, and 
U*4 inversely as the difference in mass; it is 2.40+ ~ 0.05. 

Values of v for the 14-Mev fission of U8 have been reported as follows: 

4.55+0.15, Gaudin and Leroy (1958); 
4.50+0.32, Flerov and Talyzin (1958); 
4.45+0.35, Flerov and Tamanov (1958); 
4.12+0.25, Johnstone (1956). 


*A rather lower value 1.13+0.05 barns at 14.6 Mev has been reported by Flerov, Berezin, 
and Chelnokov (1958). 

tThis value does not include the contribution from inelastic scattering events in which the 
neutron loses less than 0.5 Mev. It is the appropriate number to use in conjunction with the 
value given earlier for the total non-elastic cross section, since this also fails to include such 
low-energy-loss events. 
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The first of these measurements needs correcting for the revision of the 
value of v in the thermal neutron fission of U***, The second needs a correction 
for the effect of the (m,n’) and (n,3n) cross sections, which were ignored. 
With relative weights of 4, 1, 1, and 2, the weighted mean is 4.32+0.12, 
which agrees well with an unpublished value of 4.31+0.12 due to Allen (1960). 

If it is assumed that the number of postfission neutrons per fission NV; is 
independent of mass number and increases linearly with excitation, then the 
values of N, for U?*® at 18.86 Mev and U2 at 13.1 Mev can be calculated from 
the value of v(= N,) for U7 at 6 Mev and the value of v (which includes 
prefission evaporation) for the 14-Mev neutron fission of U***. The value of 
ON,/dE deduced from the change of v between 6 and 18.86 Mev is0.142 Mev—. 
This value, which has been used above to correct v values to a nuclear excita- 
tion of 6 Mev, is only slightly less than the value of dv/dE since prefission 
evaporation requires about 6.5 Mev per neutron (see Table I). Values of 
0v/0E in other nuclei are similar (Leachman 1958). 








VERTICAL DISTRIBUTION OF SODIUM 
IN THE UPPER ATMOSPHERE! 


W. R. BuLLock AND D. M. HUNTEN 


ABSTRACT 


The sodium twilight has been observed over the period October 1958 to 
December 1959 with a birefringent photometer. Abundances and vertical dis- 
tributions have been found for 120 twilights. The abundance variation is similar 
to that found previously but is more nearly symmetrical about the solstices. 
A small seasonal variation in peak height is found, 3 to 5 kilometers with maxi- 
mum in March. In general, changes in abundance cannot be correlated with 
changes in distribution; however, a few exceptions to this were found in June 
and July, when a large overnight increase in abundance accompanied a large 
decrease in the height of the whole distribution. 


1. INTRODUCTION 


There are many unsolved problems in the attempt to understand the 
behavior of atmospheric sodium. These have been discussed by Hunten 
(1956b), Omholt (1957), and Blamont, Donahue, and Stull (1958). The 
observations show a large abundance variation with its maximum in winter 
and minimum in summer; it would not be surprising if this were accompanied 
by a change in height or vertical distribution, but Hunten (1956a) found 
that no significant change in height could be observed at Saskatoon. This 
puts severe restrictions on any proposed explanation of the abundance varia- 
tion. Since greatly improved methods of observation and of treating the 
observations are now available (Rundle, Hunten, and Chamberlain 1960; 
Hunten 1960), a further investigation seemed justified. This is all the more 
true since Montalbetti and McEwen (1959) have observed a considerable 
height variation at Churchill, which is only 7° further north than Saskatoon. 
Moreover, their abundance variation has a somewhat different shape from 
the one observed at Saskatoon, and a still different shape has been found by 
Blamont and Donahue (1958) at a latitude of 44°. The new measurements 
show a small but apparently significant height variation of 3 to 5 kilometers 
with maximum in March; this is exactly out of phase with the smaller variation 
found by Hunten (1956a), which was not considered significant. The abun- 
dance variation is now found to be nearly symmetrical about the solstices 
and thus resembles the results of Montalbetti and McEwen. 


2. OBSERVATIONAL METHOD 

The photometer which was used has been described by Rundle, Hunten, 
and Chamberlain (1960). It includes a birefringent filter of the type now 
widely used for airglow studies, mounted in a large vertical telescope which 
observes near the zenith. Electrical compensation is used to eliminate the 

‘Manuscript received February 10, 1961. 

Contribution from the Department of Physics, University of Saskatchewan, Saskatoon, 
Saskatchewan. Supported by the Geophysics Research Directorate of the Air Force Research 
Division, Air Research and Development Command, under Contracts No. AF 19(604)-1831 
and AF 19(604)-7265. 
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signal caused by the Rayleigh component of the twilight. The field of view 
is rotated around the zenith to permit a direct measurement of the derivative 
of sodium intensity versus shadow height; we will call this AB. The derivative 
may also be calculated numerically from the ordinary record of total intensity; 
this version will be called B. In principle, as explained by Rundle, Hunten, 
and Chamberlain, AB is preferable to B, but we have found in practice that 
AB is less useful because its signal-to-noise ratio is lower. 

B was found from the chart records in the following way. Intensity (B) was 
measured at 1-minute intervals and punched into computer tape. The computer 
calculated the geometrical shadow height (z,) for each of these times, and then 
the ratio [B(t)—B(t+3 min)]/[z:(t) —2:(t+3 min)] which can be taken as B 
(that is, dB/dz,) for the average z;. To find the distribution, the procedure 
given by Hunten (1960) was used: B was plotted against (2:+@) with # taken 
in the range 26.3 to 27.0 km according to season, and a smooth curve drawn.* 
The self-absorption correction tabulated by Hunten was applied when the 
abundance was large. Then the Bracewell chord-correction method was 
applied to unfold the distribution from the curve, using a chord appropriate 
to the season. The AB data were plotted in the same way, after removing the 
geometrical factors discussed by Rundle, Hunten, and Chamberlain (1960), 
and the Bracewell treatment applied. 

A typical result is shown in Fig. 1 where the results from B and AB are seen 
to agree well. Both curves show negative densities around 70 km; this effect 
was general in all the runs made and was also found by the matrix method 
used by Rundle, Hunten, and Chamberlain. The Bracewell procedure does 
produce small oscillations as shown by Hunten (1960), but these are much 
smaller than the one in Fig. 1. The error was finally traced to an instrumental 
defect: during the bright part of twilight, the noise level from the white light 
is very large and was overloading the phase-sensitive detector. The result was 
to reduce its sensitivity to the sodium signal; therefore, the derivative goes 
to zero too quickly and may even go negative. The effect was demonstrated 
experimentally with an artificial interfering signal to simulate the noise, and 
found to be of the right order of magnitude. One or two twilight runs were 
made at reduced sensitivity; they gave good distributions as expected with 
no negative densities. 

Since the reason for the error was not discovered until after all the obser- 
vations had been made, it was necessary to find a method of correcting them. 
The overload is serious only when the white light is bright, that is, for de- 
pression angles less than about 7° corresponding to a shadow height of 48 km 
or a sodium height of (48+) = 75 km. Therefore, the part of the distribution 
above this should be little affected, and the part below, as observed, too small 
or even negative. A reasonable correction can be made by extending the curve 
from 75 km smoothly to the axis at 55 or 60 km, since there is good evidence 
from earlier work that there is little or no free sodium below this height. Such 
corrected curves are shown in Fig. 1. 

The density scale was found by measuring the “plateau” brightness corre- 


*x can be regarded as a screening height, that is, the effective height of the shadow of the 
atmosphere and ozone layer. 
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Fic. 1. Sodium distributions measured January 24, 1959, in morning twilight. The one 


labelled B was found from the derivative calculated from the total intensity; the one labelled 
AB, from the directly measured derivative. As described in the text, the dotted lines are 
believed to give a nearly correct distribution from 50 to 75 km. 


sponding to a solar depression angle of 6.5°; at this angle, the whole sodium 
layer is fully illuminated. Because of the instrumental error just discussed, 
the values recorded were too low and were multiplied by a factor 1.15. This 
was estimated from several distributions similar to Fig. 1 by comparing the 
areas under the corrected and uncorrected curves, counting the part to the 
left of the axis as negative. 

The resulting brightness was converted to an abundance by Fig. 3 of 
Chamberlain, Hunten, and Mack (1958). The calculations of Blamont, 
Donahue, and Stull (1958) are probably more accurate, but the difference is 
less than the uncertainty in the solar intensity; we preferred to use the results 
of Chamberlain, Hunten, and Mack to facilitate comparison with our earlier 
measurements given there. 

RESULTS 
Average Distributions 

These were found by averaging the data for B and applying the chord cor- 
rection to the result. It would perhaps have been preferable to average the 
distributions themselves, but this would have been much more difficult 
because they were available in graphical form only. The average distribution 
for each quarter is shown in Fig. 2; a strong seasonal variation of abundance 
is present and there is indication of a small height variation. The shape of the 
distribution does not change appreciably and it seems to be symmetrical 
about the peak. 
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Fic. 2. Average distributions for each season of the year. 
Seasonal Variations 


The abundances found are plotted in Fig. 3 along with a sine curve found 
by harmonic analysis of the monthly means. The results observed previously 


cS on oa 


SODIUM ABUNDANCE (10° atoms/cm’) 
w 





J F M A M J J A Ss 0 N D 


Fic. 3. Annual variation of abundance. Large squares denote monthly means, and the 
sine curve is from their harmonic analysis. The other two curves are the ‘‘extremes”’ trans- 
ferred from Fig. 4 of Chamberlain, Hunten, and Mack (1958). 
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at Saskatoon (Chamberlain, Hunten, and Mack 1958) are indicated by the 
two “‘extreme”’ lines given in Fig. 4 of that paper; actually, about half a dozen 
points lay outside these ‘‘extremes’”. The agreement of general magnitudes 
in the two sets of measurements is surprisingly good, considering the difficulty 
of standardizing a calibrator containing a sodium lamp. The large “spring 
hump” which was so evident in the earlier results has disappeared, and the 
sine curve drawn through the monthly means represents them well. Its 
extremes are 4.710% atoms/cm? and 1.010%, about 3 weeks after the 
solstices. These times cannot be considered accurate to better than 3 or 4 
weeks. Part of the ‘‘random”’ variation about the average may be due to 
observational difficulties, but most of it is believed to be real. The largest 
abundance was about 710° atoms/cm? and the smallest about 0.6 X 10°. 

A quantity called ‘70% peak height’ was found from each distribution; 
it is the average height of the two points having densities 70% of the peak 
value. Such a measure is expected to be less affected by noise than the position 
of the peak itself. The harmonic analysis gave a maximum height of 88.0 km 
in March and a minimum of 85.2 km in September. The average distributions 
shown in Fig. 2 suggest a slightly larger variation, from 89 to 84 km. Individual 
heights were mostly between 82 and 92 km, with a few outside these limits. 

The ‘20% thickness’ of the layer was defined as the height difference 
between the two points at 20% of the peak density. There was some sign of 
seasonal variation, but it is not considered significant; the mean value is about 
29 kilometers. 

“Random” Variations 
All the quantities measured showed large variations about the means for 
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Fic. 4. One of the summer events in which a sudden large increase in abundance seems 
to be confined to the lower part of the sodium layer. 
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a given time, a behavior which has been reported in many earlier papers. 
The variations are not really random, since there are many occasions when 
two or three successive results agree closely, but they are certainly unpredict- 
able. Almost without exception, no correlation could be found between the 
sudden changes in the various parameters described above. Large overnight 
increases in abundance occurred once in June and four times in July; on all 
five occasions, a proportional decrease in peak height was found, so that the 
density above the original peak was almost unchanged (see the example in 
Fig. 4). It is interesting that a similar occurrence was observed by Rundle, 
Hunten, and Chamberlain (1960), but in May and over an interval of 3 days. 
In these few cases, the change in abundance appears to have been caused by 
a lowering of the height at which sodium compounds become almost com- 
pletely dissociated, but this mechanism does not explain the changes observed 
during the rest of the year, nor the seasonal variation. 


4, DISCUSSION 


The seasonal variation of abundance at Saskatoon (52° N.) has been given, 
from earlier measurements, by Chamberlain, Hunten, and Mack (1958). 
The present results are similar except that they do not show the pronounced 
hump at the end of March and give an average variation which is much more 
nearly symmetrical about the solstices. At Churchill (59° N.) McEwen and 
Montalbetti (1959) have found a similar variation whose amplitude may be 
slightly smaller. But at Haute-Provence Observatory (44° N.) Blamont and 
Donahue (1958) find a considerably different curve with two very high maxima 
in December and March and a low, almost constant value for the rest of the 
year. (If their peaks were plotted on Fig. 2, they would be off-scale by factors 
of 2} to 4.) And recently Donahue (1960) has found still different results at 
23° N. in the Sahara. We thus have indications of a strong latitude variation, 
much stronger than the hemisphere-to-hemisphere effect which has been 
postulated up to now. 

This suggestion is reinforced by the height measurements. The Churchill 
results give a variation of 10 kilometers with minima in February and No- 
vember, maxima in May and December. At Saskatoon we find a smaller 
variation, 3 to 5 kilometers, with only one maximum in March. Earlier results 
(Hunten 1956a) gave a still smaller variation of just the opposite phase; it 
was not considered significant since the measurements then available were 
much poorer than the present ones. At Haute-Provence, Donahue (1960) 
gives a variation of the same magnitude as at Churchill but with a deep mini- 
mum in June when Churchill is near its major maximum. The behavior in 
the Sahara is even more violent and does not resemble any of the others. 

It thus appears that a detailed study of the seasonal and geographical 
variations will require a network of many stations and several years’ obser- 
vation. Particular attention should be paid to the southern hemisphere. The 
effort required for such a program is considerable, but the results may shed 
light on the meteorology of a part of the atmosphere which is difficult to study 
by other means. 
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RADIO-FREQUENCY EREAKDOWN CONTROLLED BY DRIFT OF 
ELECTRONS IN AN INHOMOGENEOUS MAGNETIC FIELD! 


L. T. SHEPHERD AND H. M. SKARSGARD 


ABSTRACT 

A study has been made of r-f. breakdown in which the controlling loss 
mechanism arises from the drift of electrons in an inhomogeneous magnetic 
field. The study was carried out using a toroidal system with parallel r-f. electric 
and steady magnetic fields. An approximate average-electron theory of drift- 
controlled breakdown is presented. Experimental measurements of breakdown 
r-f. electric field versus magnetic field were made at various pressures from 1.25 
to 6.0X10-% mm of Hg, using hydrogen and helium gases. A radio frequency of 
8 Mc/sec was used. Magnetic fields up to 2000 gauss were employed. The r-f. 
breakdown field was found to vary as the inverse square root of the magnetic field 
as predicted by the theory. 


I. INTRODUCTION 


Breakdown in a gas involves a sudden, large increase in the density of 
electrons and positive ions. In r-f. breakdown the production phenomenon is 
the primary ionization caused by relatively few electrons, initially present 
in the gas, receiving power from the r-f. field. At the same time electrons are 
lost from the gas by one or more loss mechanisms. The competition between 
the production and loss mechanisms determines the conditions required for 
breakdown. When the production rate exceeds the loss rate a sudden build-up 
of the degree of ionization occurs. 

Various types of breakdown in gases have already been studied by several 
workers. A review of this work has recently been made by Brown (1959). 
Most of the previous work has been done on diffusion-controlled high-frequency 
breakdown in which the dominant loss mechanism is diffusion of the electrons 
through the gas to the walls of the chamber. 

In this paper a study is reported of r-f. breakdown controlled by the loss of 
electrons due to drift in an inhomogeneous magnetic field. So far as is known no 
previous study of this type of breakdown has been reported. The work was 
initiated in connection with an experiment being planned at this laboratory 
on betatron acceleration of electrons in plasma (Skarsgard 1959). In this 
experiment the plasma is to be formed within a toroidal vacuum chamber by 
striking a r-f. discharge in a gas at low pressure. In previous experiments of 
this kind (Reynolds and Skarsgard 1958, 1959) a similar scheme of plasma 
formation was successfully employed without much understanding of the 
process of r-f. breakdown involved. 

The study of r-f. breakdown has been carried out with a system shown in 
Fig. 1. A discarded Allis-Chalmers betatron donut was employed as the 
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to pump 










3 of 18 equally spaced 
coils providing B 


discarded betatron donut 
coated on inside with thin 
palladium conducting layer 


probe 


gas inlet 


2- turn winding for inductively 
coupling r-f. field E into gas 


Fic. 1. Toroidal system used for study of r-f. breakdown, 


vacuum chamber. Breakdown at low pressure was facilitated by a “‘steady”’ 
magnetic field B parallel to the direction of the r-f. electric field E. The coils 
used to provide these fields are indicated in Fig. 1. 


II. THEORY 


An accurate theoretical description of r-f. breakdown can only be given in 
terms of the velocity distribution function for the electrons. Here an average- 
electron theory is presented in an attempt to predict the broad features to be 
expected in r-f. breakdown controlled by drift of electrons in an inhomogeneous 
magnetic field. The rates of ionization and loss of electrons are determined in 
terms of average properties of the electrons. Since the degree of ionization is 
extremely small at the beginning of breakdown, conditions for breakdown 
can be examined without considering such plasma phenomena as ambipolar 
diffusion, plasma oscillations, or instabilities. 


II.1. Average-electron Theory of Drift-controlled Breakdown 

The magnetic field B of Fig. 1 is inhomogeneous, increasing toward the axis 
of the system as shown. When an electron moves in an inhomogeneous field 
of this nature it experiences a drift in the direction of VBXB (Spitzer 1956) 


with a speed 


9 


9 
Ww WwW) 
1 Ww = —-4+—L 
(1) -  Seagkt * anae’ 
where w, = gB/m is the cyclotron frequency for the electron, 
R is the radius of curvature of the magnetic lines of force, 
W,,w, are the components of electron velocity perpendicular and 


parallel to the magnetic lines of force, and 
Wa is the speed of drift of the instantaneous center of gyration. 
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Mixed Gaussian units are employed unless otherwise noted. 

Equation (1) is derived from a first-order theory and its validity demands 
that B be large enough so that wg < w, the particle velocity. The first term 
on the right of (1) is the so-called “gradient drift’’ arising from VB. The 
second term is a consequence of the centrifugal force experienced by the 
electron moving along the curved magnetic lines of force. 

The electron drift given by (1) leads to a loss of electrons from the gas 
inside the chamber. The name “‘drift-controlled breakdown”’ is given in this 
paper to r-f. breakdown in which the controlling loss mechanism is this drift 
of electrons toward the walls of the chamber. Breakdown occurs when the 
rate of production of electrons due to ionization exceeds the rate of loss of 
electrons due to drift. The condition for equality of the rates of production and 
loss is called the breakdown condition and may be expressed by the equation 


(2) y= Varitt 


where v, is the average rate of production of electrons per electron due to 
ionization and vg, is the average rate of loss of electrons per electron due to 
drift. 

First consider v;. In the absence of collisions an electron oscillating in a 
r-f. electric field gains no energy on the average. However, elastic collisions 
with the neutral gas atoms alter the phase relationship between the r-f. field 
and the electron velocity so that power is received by the electrons. Kinetic 
energy of directed motion gained from the r-f. field between collisions is 
transformed into kinetic energy of random motion upon collision. In this way 
the kinetic energy can be increased in successive collisions until the electron 
has acquired sufficient energy to ionize a gas atom. The average power received 
by an electron from a r-f. electric field (Brown 1959) is given by 





(3) pat, 


where v,, is the frequency of collision for momentum transfer, 
w is the angular frequency of the r-f. field, 
E is the r.m.s. field. 

If all the power goes into ionization, the ionization rate per electron is 
P/qu, where u, is the ionization potential for the gas under consideration. In 
general, part of the power gained from the r-f. field will be spent in excitation 
of the gas atoms resulting in subsequent emission of optical radiation. 
Generally, also, part of the power goes into heating up the gas atoms through 
the elastic collisions with the electrons; however, for the low pressures con- 
sidered here power dissipated in this way can be shown to be negligible. In 
any case 
P on: - me. 








(4) — Kiqu; s Kymu, Vato 


where K;, is a factor, greater than unity, which accounts for the fact that only 
part of the power goes into ionization. It is assumed that K, is constant for 
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any given gas. In equation (4) », should be averaged over all the electrons. 
However, for hydrogen and helium, the gases considered in this work, vm is 
approximately independent of electron energy, at least for energies above a 
few electron volts (Brown 1959). 
Now consider vgn. If Z is the average distance to the chamber wall in the 
direction of the electron drift 
(Wa )average 


(5) Varitt = 3 


Averaging the drift velocity given by (1) over all the electrons gives 


a: 
Kov 


(6) Varitt = woRE 


where v is the root mean square electron velocity and K, is a factor, assumed 
constant, which depends on the electron velocity distribution. For a purely 
random velocity distribution K, = %. However, for the low pressures studied 
in this work there is significant peaking of the velocity distribution in the 
direction of the r-f. electric field. This leads, furthermore, to a non-zero 
average magnetic force on the electrons due to the r-f. magnetic field produced 
by the r-f. winding. This magnetic force affects the average electron drift 
velocity and consequently the value of Kz to an extent which is difficult to 
determine for the system used. 

Substituting for w, and putting v? = 2qiu/m gives 
(7) Varitt = 2K a : 
Equations (2), (4), and (7) then determine Ey, the electric field required for 
breakdown. 
(8) Ey = ae ae 

gGBRZ wm 

where K = (K;K.)}. This result is valid only if B is large enough to satisfy 
the criterion for validity of the first-order drift theory employed. Thus, the 
apparent singularity at B = 0 has no significance. 

For the system used (Fig. 1) R is 19 cm and Z is approximately 2 cm. For 
hydrogen, taking @ = u, = 15.4 volts gives 


4 st 2 2\ 1/2 
(9) es eee Keo =t#*) 





cm Vm 


The value of vy, depends on the gas pressure according to the approximate 
relationship », = 5.93X10° p where p is in mm of Hg. For helium, taking 


&@ = u, = 24.5 volts gives 
. { volts 1.34X107° _. v2 ty?\ 2 
wes z,( 2) , Se a 7 


For helium », = 2.37 X 10° p. The ionization potentials and collision frequencies 
for hydrogen and helium are taken from Brown (1959). 











% 
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It is notable that the predicted breakdown field for both gases varies 
inversely as the square root of the magnetic field. This relationship appears 
to be unique for drift-controlled r-f. breakdown. 


II.2. Possible Importance of Other Loss Mechanisms 
Here the relative importance of loss mechanisms other than electron drift 


in the inhomogeneous magnetic field will be examined. 

Diffusion 

If diffusion of electrons to the walls of the chamber is the dominant loss 
mechanism the loss rate can be found from a solution of the diffusion equation. 
When a magnetic field is present the diffusion loss rate per electron, vai, can 
be expressed (Brown 1959) in the average-electron approximation by 


(11) Yau = D/A. 

where D = (v?/3vm)average iS the average diffusion coefficient, assuming a 
random velocity distribution for the electrons, and A, is the equivalent dif- 
fusion length for the system. For a cylindrical chamber of length Z and radius 
r, with a longitudinal magnetic field, A, is given by 


1 vy, (2.405\" (z) 
9 >= =5>-5 . 
(12) Ae za ( r + L 


The magnetic field effectively increases the chamber dimensions perpendicular 


to the field, reducing the diffusion to the walls. 
For the present purposes the toroidal chamber of Fig. 1 may be represented 
by an endless cylinder of mean radius 7 between 2 and 3 cm. Equation (12) 


then gives 





1 a3 
13 > = a5- 5 ae 
( ) Ae Vin + We A 
with A approximately equal to 1 cm. If », is independent of electron energy 
equations (11) and (13) give 


‘ : 
v Te 1 
a. 


(14) Yau = 3 To A? 
From (6) and (14) 


(15) Vari, _ 3K2A* 1+ (we/%m)” 
Vaitt RZ We /Vmn 


In obtaining this result the possibility that one loss mechanism affects another, 
through its influence on the electron velocity or space distribution for example, 
is ignored. It should be noted, however, that both vai, and var have the 
same dependence on the mean square velocity so that the ratio given by (15) 
is not explicitly dependent on the mean electron energy. 

A rough estimate of the ratio of drift and diffusion loss rates for the system 
of Fig. 1 may be obtained by putting in equation (15) A = 1 cm, R = 19cm, 
Z = 2cm,and K, = 1. For hydrogen, taking B = 200 gauss and p = 2.5X10™ 
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mm Hg, equation (15) then gives varir/vair = 19. Similarly for helium, with 
B = 200 gauss and p = 6X10-* mm Hg, equation (15) gives varite/vair ~~ 20. 
For lower pressures or higher magnetic fields these ratios increase. Thus for 
the conditions under which breakdown measurements have been carried out, 
loss of electrons due to drift in the inhomogeneous magnetic field should be 
roughly an order of magnitude more important than loss by diffusion. 


Electron Attachment 

The possibility of electron attachment to neutral gas atoms providing a 
competing loss mechanism has been examined. Oxygen has a relatively large 
cross section for electron attachment. Estimates for the concentration of air 
impurity in the gas together with measured electron attachment cross sections 
for air (Harrison and Geballe 1953) were used to calculate the rate of loss of 
electrons to be expected from this phenomenon. It was concluded that electron 
losses due to attachment would be approximately two orders of magnitude less 
than drift losses. 


III. EXPERIMENTAL DETAILS 


Observations of r-f. breakdown have been made in hydrogen and helium. 
Commercial grades of these gases were used; the specified impurity content 
was less than 0.2 mole per cent for hydrogen and less than 0.01 mole per cent 
for helium. Measurements of breakdown r-f. field, Ey, versus magnetic field, B, 
were made at various pressures in the range 1.25 to 6.0X10-* mm Hg. A 
r-f. frequency of 8 Mc/sec was used. Magnetic fields up to 2000 gauss were 
employed. 

A photograph of the apparatus is shown in Fig. 2. Initially gas is allowed 
into the vacuum chamber (donut) until the desired pressure is reached. The 
breakdown field is found by using pulsed r-f. and magnetic fields. Each pulsed 
operation is initiated by a push button, controlled by conventional electronic 
trigger and delay circuits (Fig. 4), and proceeds as follows. The push button 
activates a trigger circuit, which, in turn, initiates the time sequencing and 
also fires the B supply, connecting a charged condenser bank across the coils 
used to provide the magnetic field. The current through these coils reaches a 
maximum after ~1 msec, at which time the r-f. oscillator is turned on for a 
period of ~200 usec. The magnetic field B is nearly steady during the time the 
r-f. ison (Fig. 3). The r-f. field is increased in successive pulsed operations until 
breakdown is observed. Then the donut is flushed out with the gas under study, 
pumped down to a pressure of ~10~‘ mm Hg, and the cycle is repeated for a 
different value of B. 

The breakdown can usually be observed visually. However, a probe method 
was found more reliable. The unshielded tip of the probe extends into the 
donut (Fig. 1) and is connected through a high impedance to the grounded 
coating of the donut. At the moment of breakdown the probe (and presumably 
the plasma in contact with it) swings negative by as much as several hundred 
volts, providing a clearly recognizable signal on an oscilloscope trace. The 
generation of this large voltage is not at present well understood. 
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Fic. 2. 
(donut). 


Fic. 3. Oscillogram showing the time scale for the magnetic field and the r-f. electric field. 


The envelope of the pulsed r-f. is superimposed on a Rojofsky coil measurement of the current 
through the coils providing the magnetic field. 


Photograph of part of the apparatus including the toroidal vacuum chamber 
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Fic. 4. Block diagram of electronics. 


The B supply employs a 20,000-uf bank of electrolytic condensers with a 
voltage rating of 450 volts. A 5550 ignitron switch is used to connect the 
condenser bank to the parallel connection of 18 coils, of 70 turns each, equally 
spaced around the donut. High current silicon diodes (International Rectifier 
type 45LR70) are used to protect the electrolytic condensers against reverse 
voltage. Calibration of the B supply was carried out by Rojofsky coil measure- 
ments of the maximum current versus initial condenser voltage. A check was 
made on this calibration by direct measurement of the magnetic field inside 
the donut, using a small test coil. 

The r-f. circuit (behind the donut in Fig. 2) consists of a Colpitts oscillator 
followed by an amplifier whose tank circuit is tuned to the frequency of the 
Colpitts oscillator. The two turns wound around the donut (Fig. 1) provide 
the inductance of the tank circuit and inductively couple the r-f. electric field 
into the gas. The duration of the r-f. is controlled by the high voltage pulse 
providing the plate voltage for the Colpitts stage. The value of the r-f. electric 
field is controlled by the adjustable voltage on a 5-uf condenser providing the 
plate voltage for the amplifier stage. Plate voltages of up to 16 kv were 
employed on the Eimac 6C21 triode used as the amplifier tube. Corresponding 
r-f. fields of up to 6 volts/cm (r.m.s.) were induced inside the donut. Calibration 
of the r-f. field versus amplifier plate voltage was based upon measurement of 
the r-f. voltage induced around a single-turn coil inserted around the interior 
of the donut. This determination therefore gives the average electric field 
around the donut. Since the voltage around the single turn was inconveniently 
large it was compared, at low amplifier plate voltage, with the r-f. voltage 
picked up by a small coil coupled loosely to the r-f. turns; this coil was then 
used to calibrate the oscillator. 

The vacuum system included a mercury diffusion pump with liquid air trap 
and an isolation and baffle valve. Pressures were measured by a pirani gauge 
calibrated against a McLeod gauge and an ionization gauge. 

The various calibrations involve uncertainties, especially in the absolute 
magnitude of the quantities involved. Estimated probable errors are: 
E~ 15%, p ~ 10%, B ~ 15%. 


IV. RESULTS AND DISCUSSION 


Figures 5 and 6 show the experimental results together with the normalized 
results of the simple average-electron theory of drift-controlled breakdown. 
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Fic. 5. Breakdown r-f. electric field versus magnetic field for hydrogen. Points are 
measured values. Straight lines are theoretical with K(q,.) = 7.0 in equation (9). 

Fic. 6. Breakdown r-f. electric field versus magnetic field for helium. Points are measured 
values. Straight lines are theoretical with K(ue) = 4.8 in equation (10). 


On a log plot the predicted dependence of breakdown field on the magnetic 
field appears as a straight line of slope minus one-half. The normalization of 
the theory for hydrogen was made by arbitrarily drawing a line of slope minus 
one-half through the experimental points for a pressure of 2.5 microns. This 
gives Kiy,) = 7.0 in equation (9). Similarly, for helium, fitting the theory to 
the experimental results for a pressure of 6.0 microns gives Kye) = 4.8. 
Within the accuracy of the measurements, indicated by the scatter of the 
points, the breakdown field appears to vary as the inverse square root of the 
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magnetic field as predicted by the theory. Also the predicted increase of 
breakdown field with decreasing pressure is in fair agreement with the experi- 
mentally observed dependence, although with measurements at only two 
pressures for each gas a conclusive test of the theory in this respect cannot be 
claimed. It should be pointed out also that the dependence of breakdown field 
on pressure, indicated by (9) and (10), is not unique for drift-controlled 
breakdown; it only describes the effectiveness of energy transfer to the electrons 
from the electric field through collisions and is therefore not directly connected 
with the loss mechanism. 

Without more knowledge of the electron velocity distribution the significance 
of the values obtained for Ki.) and Kye) cannot be determined conclusively. 
The larger value obtained for K iq.) is possibly connected with greater excitation 
losses of power in hydrogen. In hydrogen there are excitation levels well below 
the ionization potential whereas in helium the lowest excitation level is at 
19.8 volts, not much below the ionization potential. 

The maximum electric field available with the present system determines 
the lowest pressure at which useful breakdown observations can be made. 
Measurements were attempted at higher pressures but results were irrepro- 
ducible and inconsistent with the results at lower pressure. The reasons for 
this are not clear. There is evidence, at pressures as high as 10-' mm Hg, that 
breakdown occurs at localized regions within the donut. This suggests regions 
of stronger r-f. electric field due, perhaps, to nonuniformities in the conductive 
coating on the inside of the donut. Such effects could lead to localized dis- 
charges provided the oscillation amplitude of the electrons is small compared 
to the geometrical scale of the nonuniformities. At pressures and r-f. fields 
indicated in Figs. 5 and 6 the estimated oscillation amplitude is several 
centimeters, and decreases with increasing pressure (and decreasing r-f. field). 
Since this is small compared with the circumference of the donut, difficulties 
may be expected due to irregularities in the r-f. field and could possibly 
account for the curious behavior at higher pressure. 

The considerable spacing of the coils providing the magnetic field (Fig. 1) 
leads to inhomogeneities in the field along the direction of B. The system 
resembles somewhat the “bumpy torus” investigated by Gibson, Jordan, and 
Lauer (1960). Under certain conditions and over restricted regions of the torus 
cross section these workers found closed precessional orbits for the single 
particle drift motion perpendicular to the magnetic field. Such effects of 
‘“‘bumpiness”’ in the system of Fig. 1 are estimated to be small but they never- 
theless add to the uncertainty in the value of the constant Ke, introduced in 
the expression for vari, the rate of loss of electrons due to drift. 
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THE PRESSURE PULSE PRODUCED BY A LARGE EXPLOSION 
IN THE ATMOSPHERE! 


V. H. WEsToN 


ABSTRACT 


The pressure pulse produced by a large explosion in the atmosphere is investi- 
gated. The explosion is represented in terms of the excess pressure and normal 
velocity on a closed surface, outside of which the hydrodynamical equations are 
linearized. The pulse is represented in terms of a Fourier transform of the 
associated harmonic frequency problem, for which a ring-source Green’s function 
is obtained in terms of an expansion of the discrete modes. It is shown that the 
excess pressure may be represented in terms of an integral (containing the 
Green's function) over the surface surrounding the source. The gravity wave 
portion of the pressure pulse at the ground is computed for various ranges from 
the source, which is located at various altitudes, and for three models of the 
atmosphere. In calculating the head of the pulse a new asymptotic technique is 
introduced which gives very good results for intermediate and long ranges. 


1, INTRODUCTION 


There have been many recent efforts devoted to the study of ‘‘free’’ waves 
in the atmosphere and the pressure pulse produced by large explosions such 
as the Krakatoa volcanic eruption of 1883, the great Siberian meteorite of 
20 June, 1908, and, recently, nuclear explosions. As an example of the latter, 
the Japanese have recorded such atmospheric waves from nuclear explosions 
(Yamamoto 1957). These waves travel great distances around the earth. 

The pressure pulse produced by such explosions can be decomposed into 
two portions, a low-frequency ‘gravity’? wave train, plus high-frequency 
“‘acoustic”’ wave trains. Theoretical analysis has been mainly devoted to the 
“gravity”? wave portion of the pressure pulse produced by large explosions on 
the ground for simple models of the atmosphere (Pekeris 1939, 1948; Scorer 
1950). 

Hunt et al. (1960) have calculated both the gravity wave and acoustic wave 
portions of the pressure pulse produced by explosions on the ground. 

Here the interest will be mainly concerned with the head of the gravity 
wave portion of the pulse produced by explosions not only on the ground but 
at various heights in the atmosphere. The results presented here are discussed 
in more detail by Weston (1960). 


Initial Nomenclature 
(1) Spherical polar co-ordinates (r, 0, ¢) 
(2) Radius of earth r = a 
(3) T = absolute temperature 
R = gas constant 
7 = ratio of specific heats 
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p= 
f 


= gas pressure 
= gas density 
Co = YRTo = ypo/po 
u = velocity = (u,, U9, tg) 
(4) Unperturbed state denoted by subscripts zero, i.e., To, po, etc. 
(5) Perturbed state is represented in terms of the unperturbed state plus 
a perturbation term (which has no subscript) i.e., po+p, To+T,... 


2. CHARACTERIZATION OF THE SOURCE OF THE PRESSURE PULSE 


A large explosion in the atmosphere is usually detected at large distances 
by changes in pressure. Hence the excess pressure p(r, ¢) will be obtained. At 
a sufficient distance from the source of the explosion it is known that the 
excess variables, pressure, density, and temperature are small in comparison 
with the corresponding unperturbed variables for the atmosphere. 

Let the source of the explosion be at the point (Ro,0,0) where Ro > a 
(the radius of the earth). It will be assumed that the explosive effects are 
independent of the azimuth angle ¢. The source can be enclosed by a surface 
of revolution Sp with axis of revolution § = 0 (i.e. the z axis), such that outside 
this surface, the excess pressure, density, etc. are small in comparison with 
their respective unperturbed values. The surface Sp may be a small sphere 
with center (Ro, 0, 0). However, its actual form at present will not be specified. 

In the region outside Sy and outside the earth (r = a), it will be assumed 
that the excess variables are small in comparison with their respective un- 
perturbed values. In addition, viscosity and the effects of winds will be 
neglected. 

The source of the pressure pulse will be characterized by boundary conditions 
on the surface So. To specify the problem uniquely p(r, ¢) and neu for r on So 
must be specified. n represents the unit outward normal to So. 


3. THE EQUATIONS OF MOTION 
To begin with, note that the explosive effects are independent of ¢, hence 


all the variables will be independent of ¢. 
From Ramsey (1949), we have the equations of motion 


du 1 

) qt UV U+K = ~~~ v(p+pu)+(—g, 0, 0) 
where 

K = , (—[ws+ug], [uv9e—ug cot 6], [u,tg+u9ug cot 6]), 


the equation of continuity 


(2) = (pt po) + (U+¥) (e+ px) + (p+p0)¥ *U = 0, 


the adiabatic energy equation 











ee 
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3) 2 ptpe+ta-wyiptpe = cf 2 ote+ (wry (otm) | 
and the equation of state 


(p+po) = (pt+po)R(T+T >). 


If use is made of the hydrostatic equations for the unperturbed atmosphere 
and second-order pressure and density terms neglected, equation (1) becomes 


(4) oO + (u-v)u+K ~ —1yp+(—# ge , 0, 0). 
Also equations (2) and (3) become 
(5) 2 + urpi-+pwW *U = 0 
0 , 2| 0 ’ 
(6) P+ urp ™ ¢ | 224. s| 


where the prime indicates differentiation with respect to r. 


4. THE UNPERTURBED VARIABLES 
Before proceeding further, several relations are needed involving the 
unperturbed variables po, po, and 7». From the hydrostatic equation 
po = —£po0 
and equation of state 
Po = RpoTo 


we obtain, on integration, the following 


(7) po(r) = po(a) exp| —4 “(£ ) ar| . 
pola) 9-2 “(£) | 
(8) po(r) = RT o(r) T, dr}. 
The following relation can be obtained from (8) 
é To wae 
9 Ee . r 
' Po Ss 0 Co 


5. REDUCTION OF THE EQUATIONS FOR HARMONIC TIME DEPENDENCE 
The pulse problem is reduced to that of a harmonic time-dependent problem 

by representing the excess pressure pulse in terms of a Fourier transform. 

For the present the harmonic time-dependent problem will be considered. 
For harmonic time dependence of the form e~“‘ equations (4), (5), and (6) 


become, on neglecting second-order terms in u, 


(10) —iwu+— 1ep4(t 0,0) =0, 
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(11) —iwp+u,potpoV eu = 0, 
(12) —iwp+u,ph = co[—iwp+u,pil. 
If the hydrostatic equation is used, equations (12) and (11) become 
(13) —iw[p—cop] = u,[cops+ go] 
(14) ov-u—“ p—gu, = 0. 
Po 
Using equation (13), eliminate p from equation (10) to give 
(15) —iwupy” = —L(po"p)—t,ph po” A 
where 
(16) h = 1+; (#48), 
@ po Co 
ee - 1 po 
i so [eate]. 
( ‘) ata Po 
and the operator L has components 
1010 1 a.) 
i) ae (212,453 


The velocity vector u may now be eliminated from (14) and (15) giving 


m) vL (oi *p)+2P (oa""p) = 0 
with 
(20) q(r) = ea 244'—A h’ 2 i) 
ae) = [ty \~Ai4a'-A p+, 4) J 
Let LZ, and Le represent the differential operators 
ete av) 
(21) Ly = rv (: 7 +y¥q(r) 
t wy. av) 
99 ee re) ay 
wes La ee (si O30): 


Then (19) becomes 


(23) (L,+Le)(po'"p) = 0, 
which is separable. 
The required solution of (23) must satisfy the boundary condition that the 
vertical velocity vanish at r = a, and that the total kinetic energy be finite. 
Let ¢(r) be the r component of p,;'(r) p(r), and w?—} the separation 


constant. The separated equation involving ¢(r) is then 
Lo + (4—u)o = 0 


or more explicitly, 
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om FFB) a)+e EHH -0 


The next problem is to find the particular modes resulting from the separation 
of (23), subject to the appropriate boundary conditions, and then expand 
po'*p(r, 8) in terms of these separated solutions. 


6. FREE-WAVE MODES 


In discussing the various modes (particular solutions of equation (24)), 
the following restrictions will be placed upon the models of the atmosphere that 
will be considered. First, it will be assumed that above some altitude, the 
temperature is isothermal. Secondly, since the basic problem under considera- 
tion is the investigation of the form of the pressure pulse in the lower atmos- 
phere produced by a source which is at not too high an altitude (say less than 
80 kilometers) the flat earth approximation will be made, whereby in (24) the 
terms 7? will be set equal to a’, giving 


h d d¢) Ss > 4? @ @) ae 

9° -5-—-—— — —_— j— = 
(25) & € BIA Tsar eT 
where 
(26) we — } = w*d’a2. 


Both terms uw and A will be used later in the remaining analysis depending 
on which is the more convenient. Of the two, A, being the inverse phase 
velocity of the wave, has the most physical significance. 

As was indicated above, the appropriate solution of (23) must be such that 
the total kinetic energy be finite. For the flat earth approximation, this reduces 
to the condition that po'”(r) p(r, 8) be squared integrable over the range 
aqreg ©. 

The required boundary condition at r = a is that the vertical velocity must 
vanish. This reduces to the condition that 


(27) ¢ + Ag = 


The problem represented by equation (25) and (27) is self-adjoint, but 
because of the infinite upper limit of the independent variable, it is the singular 
case. If we define the constant A by 


A= —A*/k,—-w +0 /e for h, > 0 
A = A’/(wh,) +N —-1/c3 for h, < 0 


where the subscript s refers to the values in the upper isothermal layer, it can 
be shown (Weston, to be published) that for high frequencies such that 
h > 0, and low frequencies such that h < 0, for all r in the rangea <r < o, 
there exist a discrete spectrum A; < 0, which is bounded below, and a con- 
tinuous spectrum A > 0. The modes corresponding to the discrete spectrum 
are such that they attenuate exponentially in the upper isothermal layer, 
whereas the modes of the continuous spectrum are oscillatory there. 


(28) 
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From Titchmarsh (1946) the excess pressure can be represented in terms of 
a summation over the discrete set of modes plus an integral over the continuous 
set. However, in the problem on hand, that of propagation of a pressure pulse 
around the earth, the discrete modes will play the dominant role. Henceforth 
the effect of the continuous spectrum will be ignored, and only the contribution 
of the discrete set of modes to the excess pressure will be considered. In addition 
it will be assumed that the expansion holds in the frequency domain wherein 
at some 7, h vanishes. 

For expediency, ¢,(7) will be defined as the particular discrete mode of 
equation (25), corresponding to the eigenvalue A; To each value A,, there 
correspond two values \,;, namely +{A,!. However, the positive value will 
always be taken. The negative value corresponds to the same mode, but 
represents propagation in the opposite direction. In addition the same value 
of \, holds for both positive and negative value of w. 


7. THE GREEN’S FUNCTION 
Instead of working directly with the excess pressure, we shall consider first 
a ring-source Green’s function, defined to be a solution of 


(29) vLG410G = a Bo) 
r 2rsin6r- 
such that at r = a 
dG st A 
0) S+AG = 0 


where 6 is the Dirac delta function. 
Making the flat earth approximation, equation (29) becomes 


afi2.\. A 2. a0 he #4 ey 
(31) Be 4 Feel +6)2-5 +3 h f 
_ 5(r—ro) 0) 6 (0 — Go) 
2m sin ba- 


Since we require that G be squared integrable with respect to r we can expand 
G in the form 


= $5(7) b5(70)X 5(9, 80) 


(32) G(t,r0) = > : 
Paar 


where the continuous spectrum is neglected. The denominator can be expressed 
explicitly in the form 


> 4... «ele rS@ [ 490)-+9'0) || 
(33) cj $;(r)dr = 20d, d ‘2 i se 


T=a 





Setting expression (32) in equation (31), and using equation (25) together 
with the orthogonality property of ¢,; we obtain 
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Oe a ices. ee 
(34) sin 0 00 (sino 2x,)4a'a MX3 = So aind’ 


By usual techniques the solution is found 





1 {oes 6) P_1)2+4;(—Ccos 90) 0<A< % 


5 r : = 
8) Sean P_1)24; (COS 60) P_1/24y;(—cos @) o<d<t 


4 cos TH; 


with wy = (w’a2Aj+43)!4 and P,(cos @) is the Legendre function of order n. 

A note of caution is due here. There may exist a set of values of w for which 
one or more y, is equal to half an odd integer, in which case the denominator 
of (35) vanishes. For present we will restrict ourselves to the values of w not 
contained in this set. This problem will be circumvented later. 

Using the Green’s function we may now express the excess pressure at a 
point (7, 6) in terms of an integral representing the source. Multiply equations 
(19) and (29) by G(r, ro) and po'”p(r) respectively, subtract the two resulting 
equations, and integrate over all space outside the earth’s surface and the 
surface So surrounding the source. Using the divergence theorem, the volume 
integral may be reduced to a surface integral, the result being 


(36) J n+{GLo(0"p) — po pLvG} dSo = po'*(r)p(r). 

s 
The surface S comprises of So, y = a, and the limiting surface r = ©. Because 
both po!p and G satisfy the same boundary on r = a, the surface integral 
over r = a vanishes. Since both G and pq!” are squared integrable for the 
flat earth approximation, the integral over r = © vanishes. Thus the surface 


integral in (36) reduces to an integral over So. 
Upon applying equation (15), the excess pressure has the form 


(37) po (r)p(r) = iw J. Gp (ro)m* u(ro) dS 


“ | po (ro) p (to) {(4) Gn-ietn-Let dS 


where n is the unit outward normal to the surface So. The variables of integra- 
tion are denoted by the subscript zero. The vector Lo is the vector L with the 
variables r,@ replaced by ro and 69, and fo is the unit vector in the direction 
of increasing fo. 

Using the representation of G(r,ro) given by (32), formally rewrite it in 
the form 


- ae f(f, Yo) 
am) OAR, Be) » 2 COS T py; 
and then express p(r) by the following 
(39) ope) =D 5 — 

° 7 2COS Ty; 


te ate 
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where J, is the equivalent of the right-hand side of (87) with G replaced by 
f,(r, Yo). These relations will simplify further analysis. 


8. THE EXCESS PRESSURE PULSE 
With the associated harmonic time-dependent problem established, the 
corresponding pulse problem can now be treated. The explosion is characterized 
by known values of p(r, ¢) and neu (r, ¢) on the surface So, namely 


(40) p(r, t) = 0 t<0 
= m(r, t) t>0 
(41) n-u(r,t) = 0 t<0 
= nr, b) t>0 


for r on So, with m(r, t) and n(r, t) being known functions. Let the Fourier 
transforms of m and n be M(r, w) and N(r, w) respectively. 

Then the excess pressure pulse generated by initial conditions (40) and (41) 
is given by 





(42) pte, t) = so) fe De — dw 


on 
j 2 COS T py; 


where J, are given by the relation 
(43) I; == ia f f(t, I) po (ro) N (fo, w) ds 
So 


- foo M (fo, w) {(4) a *fotn- Lal f dS. 


However, as was pointed out earlier, there may exist a set of w such that 
there exists one or more yu, for which cos z yw; vanishes. Hence the line integral 
in (42) must be deformed to pass around these values of w. 

Perform the following decomposition 


1 — 2 : (—1)" exp(i2 rny;) 
ee. ee = a 9 Nd ee Ree ees, 
rt) 2 COs 7 pb; 2, (i ep ener 2 cos 7 p; 
and set 
n—1 
(45) p(r,t) = > p™(r, +P (r, t) 
m=0 


where p”(r, t) and P(r, t) are given by the right-hand side of (42) with the 
denominator in the integrand replaced by 


(—1)" exp 2{4(2m+1)u,} and (—1)" exp (124ny;)/(2 cos my;) 


. . . 2 . 
respectively. In the exponents involving wy = (w’a?Aj+})!”, we will make the 
approximation 
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which is only inaccurate for an extremely small and negligible range of o, 
roughly given by w < 107‘ sec. Using (46), the only difference between the 
integral expressions for p(r,?f) and P™ (r,t) is that the former contains 
exp(—twt) and the latter (—1)" exp—iw[t—2zad,n]. Thus it immediately 
follows that 

(47) P(r, t) ~ (—1)"p(r, t—2xnay,). 


Since \; is the inverse phase velocity, the time interval represented by 2rna), 
is the time taken for the mode to encircle the earth n times. Hence P(r, t) 
represents the portion of the pulse that has encircled the earth at least ” times. 
Finally it may be shown that p‘”(r, ¢) represents the portion of the pulse 
that has encircled the earth m times. 

Henceforth we will be concerned with the portion of the pulse that has 
arrived directly without encircling the earth. This portion is given by 


(48) pr, 2) = po eo " z Ie patel ads] gs 


This splitting up of the pressure wave into the above components removes 
the problem of poles on the path of integration. 


9. THE PRESSURE PULSE PRODUCED BY A SIMPLE SOURCE 


The relation given by (48) together with (43), for the directly received 
portion of the pressure pulse is for a general explosion. Unfortunately, 
analytical solutions are not available for large explosions. An analytic solution 
has been found for an intense spherically symmetric explosion (Taylor 1954) 
but it holds down to only overpressure of about 20 atmospheres. Hence for 
good source models, values of the excess pressure and normal velocity must 
be obtained from observational data. For numerical evaluation, the problem 
will be simplified by letting the surface 2 be a small sphere with center (Ro, 0). 
Then the limit is taken on letting the radius shrink to zero. In addition, the 
time dependence taken will be such that the source is an instantaneous velocity 
source. Taking the limit, one obtains (Weston 1960) 


(49) I; = tof ;(t; Ro, 0) Vp" (Ro) 


where V is the total volume of gas introduced. 
The case where the observer is on the ground will be considered, namely 
r = a = (a, 6). Then from (32), (33), (35), and (38) one obtains 


w'Ay4(Ro)P- 1/2+nj (— COs 8) 


50) fi; a; Ro, 0 = 
lar h ed; 


r=a 
and except for a very small negligible interval around w = 0, the following 
approximation may be made 


(51) P14p;(—cos 0) ~ (42|wlad, sin 0)~'” cosy |wlad (e~¢)—= 
} \ J 4f 


for0 <6 < x. 
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Of course when @ = = (i.e. the observer at the antipodes of the source) the 
left-hand side of (51) is unity. 

Hence for the simple source model, the directly received pulse on the ground 
is given by 


(52) p(a, 6, t) 
1/277 fo f 
f —«{ eetRadoa(o) ) rs 7 =e 9 
—o f 


27 asin 6 
. 51 Q;(w) dw 


where for convenience 


—w'";(Ro) : 
1¢+¢' “}| ' 
heh J 


h 


fa 
lan 


Since A, and Q;(w) are even functions of w, the integral in (52) may be rewritten 
(54) 4 J Le Ved; 0,(u) cos] wt —ar)+5 | 
0 j 
\ 


00s] w(t—an,(2x—6)) +35 | f dw. 


(53) Q;) = 7; [ae 


T=a 


The first integral in (54) corresponds to the wave that has taken the shortest 
route, whereas the second integral in (54) corresponds to the wave that has 
travelled by way of the antipodal route. 


10. EVALUATION OF THE GRAVITY WAVE PORTION OF THE PULSE 

We will be concerned with the contribution to (52) of one particular mode, 
the low-frequency gravity wave mode. For particular models of the atmosphere 
there exist other modes which have a low-frequency cutoff (usually called 
acoustical modes). The observed pressure pulses produced by large explosions 
(Yamamoto 1957; Scorer 1950; Hunt et al. 1960) indicate the low-frequency 
gravity wave superimposed on which there may exist the high-frequency 
phenomena. For some models of the atmosphere (mainly those with the upper 
isothermal layer the coldest) the gravity wave mode has a high-frequency 
cutoff (Weston, to be published). 

We shall consider the portion of the gravity wave that has taken the shortest 
route. (The antipodal wave is discussed in Weston (1960).) The excess pressure 
has the form 


(55) p(a,0,t) = 4/ ee6Reeo(o) Fl Vex Q(w) cos] w(t—ane) + | dw 


2rasin@ r 


where the subscript j is omitted, and it is understood that the values of 
\, Q(w) correspond to the gravity wave mode. If the gravity wave in question 
has a high-frequency cutoff, w,, then it will be understood that the integrand 
in (55) vanishes for w > we. 
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The integral in (55) is equal to the real part of 


(56) Jove Q(w) exp if ot—ar) += | dw. 


For large ranges (a@ sufficiently large) the above integral may be evaluated 


by the method of stationary phase. 
If wo is the value of w for which 


(57) t—a6(wa)’ = 0 


where the prime indicates the derivative with respect to w, then the stationary 
phase technique gives the asymptotic approximation to (56) 


(58) 4/2 Q(ws) 4/ 


Of course if there is an upper frequency cutoff, w,, this will hold provided 
that wo < w,. In addition the stationary phase technique for (56) is not 
accurate unless (waA)’””/[(waA)’”]22 <+/a6 and (wad)'*/[(wA)” |? < a@. It has been 
shown by Hunt et al. (1960) that this critique fails for most models of the 
atmosphere under consideration when the range a@ is of the order of 1000 km 
or less, and for some models of the atmosphere, even when the range is of the 
order of 3000 km. 

For these cases it is best at present to perform numerical integration to 
evaluate (56) in the manner that is suggested by Hunt et a/. (1960). 

In addition, even at very large ranges, the stationary phase technique fails 
for the initial part of the received pressure pulse. We shall give here an 
asymptotic technique for calculating the initial part or head of the pulse for 
the gravity wave portion. This method will hold for models of the atmosphere 
where the gravity wave mode does not have a low-frequency cutoff, i.e., the 
wave exists at w = 0. 

Define the following series 





oe cos[wot — woAa8]. 
0A, 


(59) A(w) = > rw”, 
(60) (x)7Q(w) = 0370(0) z dgo™ 


where dy = 1. It can be shown that the above expansions are valid for a 
sufficiently large range of w around w = 0, for the models of atmosphere under 


present consideration. 
Define the time interval 


(61) r = t—aOXo. 


By the stationary phase technique, the condition + = 0 or ¢t = a@Ao would 
indicate the beginning of the pulse. However, for |7| small, the stationary 
phase technique fails and we shall show that + = 0 does not represent the 
beginning of the pulse, but a point near the peak of the first compression. 

The integral (56) may now be written in the form 
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(64) q 


o(Aya8)-¥3 
7(A,a9)~¥8 


Il 


and defining the integral 


(65) F,(q) = { gf Me ort 4 
v7 


expression (64) has the asymptotic form 








(66) 4/ Q(0) e “nF 1(q) a8) a oe ams F3(q)+ . 





16 
at ie 120) (A = 


It is shown (Weston 1960) that the F,(g) can be expressed in terms of 
products of Airy integrals. Thus expression (66) may be represented in terms | 
of Airy integrals (which are tabulated by Miller (1946)). After much algebraic 
manipulation (Weston 1960), the real part of (66) has the following asymptotic 
form 
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where A ;,(y) is the Airy integral, with A‘ being its derivative. The argument 
y is given by 


. = Ey 
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The functions C;(7) have the forms 
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11. NUMERICAL RESULTS 
With the theory set up in the previous sections, one can now compute the 
pulse forms. The first problem is to evaluate the radial eigenfunctions and 
eigenvalues corresponding to the gravity wave. For some simple models of 
the atmosphere, good analytical approximations to the eigenfunctions may be 
made using an iteration technique. This is discussed by Weston (1960). How- 
ever, for the models of atmosphere that will be considered here, it is best to 


use numerical techniques. 
Instead of solving for (7) we will investigate the function ¥(r) given by 


(70) v(r) = ¢(r) exp | A dx. 


Using equations (25) and (70), the resulting equation obtained for y is 
" ” ’ h’ 2 1 2 
(71) vy’ +y'| —2A—— |+hw'] 3—d* | y = 0. 
h Co 


The boundary condition at r = a reduces to y’ = 0. 
In the upper isothermal layer, the required eigenfunction y will have the 


form 
(72) v(r) ~ exp[A,—8l]r 


where 


a 2—y ra 1 Z . g° 1/2 
(73) Seis) ala § Reto 


and the subscript s refers to the values in the upper isothermal atmosphere 


where in particular 


9_ 
_ (Sere 
Aa = ( 9 ) 2 


The models of atmosphere that will be considered below will be such that 
the values of temperature are given at a set of altitudes (or data points), 
up to a certain height /, above which the temperature will be taken constant. 
For simplicity in numerical integration the temperature between two adjacent 
sets of data points will be taken linear, the gradient of temperature determined 
by the values at the end points. Hence for numerical integration, the model is 
approximated by a set of straight elements, the points of discontinuity of the 
gradient of temperature being at the data points. 

In transversing these points of discontinuity of 76 in the numerical integra- 
tion of (71), use is made of the conditions that both y and y’/h be continuous. 
These results are obtainable from the condition that the pressure and vertical 
velocity be continuous. 

At the lower boundary of the upper isothermal layer (namely at r = a+/), 
equation (72) and the upper continuity conditions are utilized to obtain 
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2— t 
v ~[-0+(252) 5] 


where y is normalized so that at r = a+/, y is unity, and h, and h, are the 
values of h below and above the discontinuity at r = a+. 

Using the results of Weston (to be published), upper and lower bounds may 
be placed on the required values of \. Hence initially 4 can be guessed quite 
accurately. 

Two models of the atmosphere are given in the following tables, where the 
altitudes are given in geopotential kilometers and the temperature in degrees 
centigrade. 


Atmosphere I Atmosphere II 


Altitude Temperature Altitude Temperature 


Surface 29 Surface 17 
1.5 i8 L.@ 8 
3.1 10 3.0 — 3 
5.9 — 5 5.6 —-19 
9.7 —3l 9.2 —46 

12.4 —53 11.9 —49 
14.2 —67 13.8 —47 
16.6 —81 16.4 —47 
20 —66 20 —48 
25 —65 25 —46 
30 —45 30 —42 
35 —34 35 —34 
40 —20 40 —18 
45 — 5 45 — 2 
50 7 4 


50 


For the sake of comparison we will include Scorer’s model (1950), such that 
the lower non-isothermal atmosphere has a constant temperature gradient 
with temperatures 286.91° K and 229.53° K at the surface and at a height of 
9.6137 km, respectively. 

The calculated inverse phase velocities for the gravity wave are given in 
Table I. Scorer’s model of the atmosphere has a high-frequency cutoff at 
around 5000 w? = 15.9. The gravity wave for the other two models exists for 
much higher frequencies than shown in the table. 

The pressure pulse forms are computed for Scorer’s model of the atmosphere 
for explosions both on the ground and at a height of 9.6137 km, and plotted 
in Fig. 1. The pressure pulses were computed using the stationary phase 
technique for r > 160 sec, and the asymptotic expansion involving Airy 
integrals for +r < 160 sec. There is good matching of the results obtained 
by both methods in the vicinity of r = 160 sec. The complete pulse forms are 
not given, but there is a cutoff at tr about 10} minutes and 17 minutes for the 
ranges 3600 km and 6000 km respectively. The pulse forms calculated here 
do differ appreciably from those calculated by Scorer. As pointed out by Hunt 
et al. (1960), Scorer’s pulse forms are in error. The first compression as com- 
puted by the asymptotic expansion (67) is the most realistic yet calculated. 
It agrees remarkably well with observed results. 
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TABLE I 


Inverse phase velocities \ (km sec) 


5000 w? + AtmosphereI Atmosphere II  Scorer’s model 








0 











3.15283 3.16379 3.18947 
1 3.15819 3.16773 3.19148 
2 3.16277 3.17072 3.19355 
3 3.16695 3.17319 3.19570 
4 3.17091 3.17534 3.19793 
5 3.17472 3.17728 3.20025 
6 3.17844 3.17907 3.20266 
7 3.18209 3.18074 3.20517 
8 3.18569 3.18233 3.20778 
9 3.18923 3.18384 3.21052 
10 3.19274 3.18530 3.21338 
ll 3.19620 3.18672 3.216388 
12 3.19962 3.18809 3.21953 
13 3.20301 3.18943 3.22286 
14 3.20635 3.19073 3.22637 
15 3.20965 3.19201 3.23006 
16 3.21290 3.19327 — 
17 3.21611 3.19450 — 
18 3.21928 3.19571 — 
19 3.22239 3.19690 _— 
20 3.22545 3.19807 — 
21 3.22846 3.19922 _ 
22 3.23142 3.20035 _— 
23 3.23433 3.20147 — 
24 3.23718 3.20257 — 
25 3.23998 3.20365 —_— 
4 ! 
3| | (a) 
2+ 
1+ 
1 + 
l | 
! | 
| i 
| 
0 5 10 ° S.. 10 15 
(Minutes) (Minutes) 
37 | (d) 
2+ 
! 
eros 
‘ thn " 
| 
! 
5 10 0 5 10 15 
(Minutes) (Minutes) 


Fic. 1. The pressure pulse at the ground at a distance of 3600 km (a) and (c), and 6000 km 
(b) and (d) for Scorer’s atmosphere, for an explosion on the ground (a) and (0), and at a height 
of 9.6137 km (c) and (d). One unit of amplitude corresponds to 0.614 wbar per 1 km of gas 
released. 


In addition, pulse forms for atmospheres I and II are computed and given 
in Figs. 2 and 3, for a range of 7000 miles. The tail of the pulse extends beyond 
the section calculated for a very long period. There is much less dispersion 
for atmosphere II than atmosphere I. 
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20 25 (Minutes) 


20 25 (Minutes) 


Fic. 2. The head of the pressure pulse at ground level at a distance of 7000 km from an 
explosion on the ground (a), and at a height of 39 km (0), for atmosphere model I. One unit 
of amplitude corresponds to 1 wbar per 1 kmé of gas released. (The tail of the pulse which 
extends for a very long time interval beyond the broken lines is not given.) 


5 


t Minutes 


Fic. 3. The head of the pressure pulse at ground level at a distance of 7000 km from an 
explosion on the ground for atmosphere model II. One unit of amplitude corresponds to 1 wbar 
per 1 km of gas released. 


12. COMMENTS 
The asymptotic methods given here give good results for the gravity wave 
portion of the pulse. To agree completely with some of the observed data, the 
higher modes must be included, although in some instances they do not appear. 
Previous authors have been limited to the investigation of simple point 
sources at ground level, where their derivation of the pressure pulse forms are 
based more on a restricted intuitive rather than a general rigorous approach. 
Although the calculations here have been confined to a simple source model, 
the effects of realistic source models in the atmosphere may be studied by 
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using equation (48) together with (43). In addition, the technique used here 
to obtain the pressure pulse form represented by equations (48) and (43), is 
not limited to models of the atmosphere for which the upper atmosphere is 
isothermal. Winds play an important role in determining the phase velocity 
of the modes. 

The effects of winds on the gravity wave are discussed by Weston and 
van Hulsteyn (to be published). 

More realistic models of the upper atmosphere are being investigated at the 


present time. 
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GAMMA-ACTIVE FALL-OUT IN ONTARIO 1958-60! 


kK. G. McNEILL, O. A. D. Trojan, AND D. J. DAwson 


ABSTRACT 


Measurements have been made of the Cs'*7 content of man and of milk and 
of the Zr® and Cs'*’ activity of air-borne dust. They show that the Cs!87 content 
of man went through a maximum of 67 wuc/g K in the fall of 1959, and. that the 
Cs'87 content of milk reached a low of 7.5 uuc/g K in September 1960. The 
results suggest stratospheric and tropospheric residence times of 2 years and 
2 months respectively. 


INTRODUCTION 

During the last decade, many measurements have been made of the radio- 
activity caused by fall-out. Most of these measurements have been made of 
the total beta activity of the fall-out (for example, by Wolfson and Terentiuk 
1960), or of one particular beta-active nuclide, particularly Sr®°°. Relatively 
few groups have studied those nuclides which emit gamma rays, largely 
because of the greater complexity of the apparatus necessary. 

During the last three years, measurements of the gamma activity of the 
human body and its environment have been carried out using the “steel 
room” at the University of Toronto. In this paper results obtained are reported 
and commented upon. 


CAESIUM-137 IN THE HUMAN BODY 

The method by which whole-body measurements of the gamma activity of 
people are made using a steel room and a 5-in. Nal crystal was described in 
a previous publication (McNeill and Green 19596). It is, in essence, the same 
as that employed at the Argonne National Laboratories (for example, by 
Miller and Marinelli 1956). 

Two groups of test subjects have been used in the Toronto work, each 
of about 25 persons without any history of contact with radioactivity. Not 
one of the subjects has shown any gamma activity except that due to K*° 
and Cs'87, Figure 1 shows the results of measurements of the Cs? content 
of the body made at the University of Toronto during the period fall 1958 
to December 1960, and at the Argonne National Laboratories since 1955. 
Following normal procedure, the caesium burden of the body is expressed in 
relation to its potassium content. The absolute values of the figures are based 
on a calibration made in 1958 with K® in the body. Substandard measurements 
were made to check the reliability of this calibration at subsequent times. 

The figure shows that the Cs'*7 burden passed through a maximum in the 
fall of 1959. This agrees with measurements made at other centers (for example, 
by Miller 1960), and it is expected as a result of the reduction of Cs!*? in 
foodstuffs consequent to the moratorium on bomb-testing. 
1Manuscript received January 31, 1961. 
Contribution from the Department of Physics, University of Toronto, Toronto 5, Ontario. 
Can. J. Phys. Vol. 39 (1961) 
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Fic. 1. The variation of the Cs!’ burden in the human body as a function of ‘time, using 
measurements made at the Argonne National Laboratories (ANL) and University of Toronto 
(U of T). The results are expressed in terms of the potassium content of the body, in caesium 
units of wuc of Cs!7/g of K. 


For 21 male students (of average age 22 years) on whom measurements 
were made in December 1960, the average Cs? burden was 7.5 X10~ uc, the 
average potassium content 0.207% of the total body weight, and the average 
Cs/K ratio 48.5 wuc/g K. The average weight of the subjects was 159 Ib. This 
Cs'*7 burden is 2.5X10-* the maximum permitted (ICRP 1959). 

A level of 7.5X10-* uc in the body is reasonably consistent with a daily 
intake of 32 wuc of Cs? via food and with a 115-day half-life for caesium 
in the body (McNeill and Green 1959a). The daily intake figure is based on 
values found in the spring of 1959 (McNeill and Trojan 19605) reduced by 
an amount equal to the reduction in Cs'*7 levels in milk between the spring 
of 1959 and December 1960. 


CAESIUM-137 IN MILK AND MEAT 


It has been shown that milk and meat are the main contributors of Cs! 
to the normal Ontario diet (McNeill and Trojan 1960). Using a method 
discussed earlier (McNeill and Trojan 1960a), we have measured the Cs}? 
content of samples of dried milk produced near Toronto since the spring 
of 1959. Figure 2 illustrates the results. 

Reference to the figure shows that the Cs!*7 content rose sharply from a 
winter plateau in May 1959, fell during the summer until October, and then 
rose to a winter plateau. In 1960 the Cs!*7 level dropped in May, continued 
to fall during the summer to a low in late September, and then rose to a 
winter plateau. These changes are associated with the feeding habits of the 
cows and the deposition of Cs'*?7 onto pasture. May is the month in which 
cows in Ontario leave the barns for summer pasture; in October they return 
to barns for the winter. In the barns part of their food is hay collected during 
the early summer, and it is therefore reasonable that the winter level of the 
Cs'87 content reflects the Cs/K ratio of the early summer. The sharp peak 
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Fic. 2. The variation of the Cs!*7 content, with time, of dried milk from the Toronto milk- 
shed. 


in May 1959 is considered to be due to the deposition of relatively large 
amounts of fall-out from the bomb tests of 1958. May 1960 did not show 
any such peak. 

The lowest value recorded was in September 1960, at which time the 
specific caesium activity was 7.5 wuc/g K, and the highest value was in May 
1959, when the specific caesium activity was 51.5 uuc/g K. 

Disregarding the transient effect in May 1959, it can be seen from Fig. 2 
that there are three factors with different time constants governing the varia- 
tion with time of the Cs'*’ content of milk. 

One of these factors is the biological half-life of caesium in the cow. When 
the animal changes diet suddenly, its total Cs? content, and consequently 
its Cs'87 output, will change at a rate governed by this half-life. (This con- 
clusion assumes that the caesium content of the body as a whole does not 
alter radically, and that the only effect is a change in the proportion of radio- 
caesium.) The rate of change of Cs'*7 in milk after the return of the cows to 
barns is consistent with a value of 10 to 20 days for the biological half-life, 
in agreement with the results of Hood and Comar (1953). 

Evidence will be presented later that the major introduction of fission 
products into the troposphere from the stratosphere occurs during the winter 
and spring, and that during the summer the rate of fall-out decreases at a 
rate governed by the residence time for tropospheric fall-out. This decrease 
in rate of fall-out should be reflected in a decrease in the Cs'*7 content of 
milk if, as is assumed to be the case, most of the Cs!” entering the cow comes 
from dust on the vegetation eaten rather than by absorption of Cs'*7 from 
the soil by plants and subsequent ingestion by animals. Figure 4, in which 
we attempt to correlate the data on milk with that obtained from air sampling, 
discussed in the next section, shows that the rates of decrease of long-lived 
fission products in milk and air during the summers of 1959 and 1960 are 
similar. During 1959, the Cs'*?7 content in milk decreased by a factor of 2.7 
in 3 months, and in 1960 the rate was approximately the same. The corre- 
sponding half-life of 2 months for tropospheric fall-out is longer than that 
found by other groups. Anderson et al. (1960) found a value of 20 days, and 
Gustafson et al. (1960) a value of 1 month. Both these values are based on 
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air sampling, and indeed the results at Toronto on air sampling, presented 
later, suggest a lower value for 1959. The difference may be due to the fact 
that not all the Cs'*’ in cows is received via pasture feeding. 

During a period free from substantial bomb-testing, the long-term variation 
of Cs!87 in milk will be governed by the rate of transfer from the stratospheric 
reservoir to the troposphere. In the period under review, any nuclear decay of 
the nuclide is insignificant. It would be expected then that the relative con- 
tents of Cs'8?7 in milk from year to year would give the half-life for residence 
of fission products in the stratosphere. Comparing the three winter plateaux, 
they are in the ratios 37:26:17. This corresponds to a half-life of 1.8 years. 

Measurements of the Cs'*?7 content of meat have not been so detailed as 
the milk determinations, but they show a trend similar to that exhibited by 
the results for milk. The Cs/K ratio in meat is approximately 1.3 times the 
Cs/K ratio in milk, showing a discrimination against caesium by the mammary 


gland. 


GAMMA ACTIVITY IN AIR-BORNE DUST 

The activity of gamma-active fission products in dust has been studied 
using a Staplex Hi-volume air sampler, with pleated filters TFA type S. It 
was found that these filters had a collection efficiency identical with that of 
filters 934-AH, which have approximately 100% collection efficiency for 
atmospheric dust (Wolfson and Terentiuk 1960). The air sampler was placed 
on the roof of the Physics Building of the University of Toronto, 80 feet 
above the ground. Toronto has a latitude of 43° 40’ and a height above sea 
level of 350 feet. The standard sampling period was 3 days, during which 
time a total volume of air of 5000 cubic meters passed through the filter (for 
comparison, it may be noted that the weekly intake of a man is 140 cubic 
meters). 

After removal from the air sampler, the filter papers were placed in individual 
plastic petri dishes, and then left for 1 week to allow radium products to 
decay. The dishes were then individually placed under a 5-in. Nal crystal 
inside the steel room and the spectrum of the gamma rays emitted by the 
samples was found. By comparison with standard sources the activities of 
the various components of the gamma rays could be obtained. It should be 
noted that all figures quoted subsequently in this section refer to ‘‘micro- 
microcuries of gamma activity” i.e., in accord with normal practice in this 
field, the decay scheme has not been taken into account. 

Some preliminary work on the activity of the air was carried out in the 
spring of 1959, and complete records were kept since the fall of that year. 
The only clear lines in the spectra can be ascribed to gamma rays resulting 
from the decay of Zr®, Ru!®, and Ru'®, and Cs'*7, 

The average monthly results for Zr® and Cs'*7 are shown in Fig. 3. The 
interpolated values for Zr® activity during the summer of 1959 (for which 
period no data are available) have been substantiated by direct comparison 
of spring and fall samples. Furthermore, the indicated rate of decrease during 
the summer agrees with that found by Gustafson et al. (1960). 











CANADIAN JOURNAL OF PHYSICS. VOL. 39, 1961 


10-0 
ke o Zr 
#€s 
: x ANL 
UO 
0-10 
S 
O-Ol 
-OOl 


FM: OA NOR OMA: GN 
1959 1960 


Fic. 3. The variation of the Zr® and Cs! activity, with time, of air-borne dust in Toronto. 
The results labelled ANL come from the work of Gustafson et al. (1960). 


Some question has been raised (Stewart et al. 1957 and Martell 1959) about 
whether the rising levels of activity found each spring from 1955 to 1960 are 
caused by the introduction of old fission products from the stratospheric 
reservoir into the troposphere through a break in the tropopause or by the 
fact that until 1959 bomb tests were carried out in the fall of each year. As 
1959 was free of bomb tests, the spring maximum in 1960 is strong evidence 
that a meteorological phenomenon is responsible at least in part. The French 
tests in February 1960 introduce complication; however, the levels were rising 
before the tests, and at least in the case of a long-lived nuclide such as Cs'*? 
the added amount of activity is small compared with the total amount already 
present. 

Since the major part of the decrease in Zr® activity is due to radioactive 
decay rather than to a reduction in the rate of fall-out, it is of interest to 
correct the results for nuclear decay using the start of the experiments, in 
February 1959, as a convenient base time. In this way the variations in the 
activity of the shorter-lived isotope may be compared directly with those of 
Cs}87, 

The corrected figures corresponding to measurements subsequent to the 
French tests of 1960 must be treated with reserve as the correction would 
greatly overemphasize any Zr® from these tests in comparison to Cs!7 from 
them. The corrected values for Zr® activity in air-borne dust from February 
1959 to February 1960, and Cs!” values from November 1959 to December 
1960, are shown in Fig. 4. 
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Fic. 4. The variation of Zr® activity of air-borne dust corrected for nuclear decay (to 
February 1959) compared with the variation in Cs!" activity of milk and air-borne dust. 


A simplified diagram of the variation of the Cs/K ratio in milk is also 
drawn in Fig. 4. If we take into account complications introduced by the 
feeding habits of cows, there is a considerable measure of agreement between 
the two types of data based on air and milk sampling. For the purpose of 
making quantitative comparisons, only the decreases in activities in the 
summers of 1959 and 1960 and the decrease in Cs"? activity between December 
1959 and December 1960 may be employed. (A comparison employing Zr® 
data before and after the French tests in February 1960 must be suspect.) 
From the middle of July 1959 to the middle of October 1959, the air-sampling 
data gave a decrease in long-lived fall-out by a factor of 7, corresponding to a 
tropospheric residence time of approximately one month. This agrees with 
the results of other workers previously quoted. In 1960, the corresponding 
figures are different, a decrease by a factor of only 2 occurring in the same 
3 months. The Cs"? figures for milk in the summers have already been given, 
i.e., a decrease by a factor of about 3 in 3 summer months in both years. 

From December 1959 to December 1960, the Cs'” content of air-borne 
dust fell by a factor of 1.4, corresponding to a 2-year half-life for residence 
of fission products in the stratosphere. This agrees very well with the value 
derived from the results for milk. 

Clearly there is no absolute correlation between the two sets of data, but 
considering the complications in this field, including, for instance, the possi- 
bility of soil-bound Cs'” being splashed onto grass and subsequently ingested 
by animals (Van Wijk and Braams 1960), the agreement is quite reasonable. 
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In summary, it may be concluded that the measurements show a gradual 
decrease in long-lived fall-out in Ontario, at a rate corresponding to a strato- 
spheric residence half-life of approximately 2 years. The tropospheric residence 
half-life would appear to be between 1 and 2 months. This decrease has 
resulted in a decrease in the Cs'*7 burden in man. 
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THE ELECTROMAGNETIC FIELDS OF A HORIZONTAL DIPOLE 
IN THE PRESENCE OF A CONDUCTING HALF-SPACE! 


JaMEs R. WalItT 


ABSTRACT 


The problem considered is a horizontal electric dipole which is located above or 
below the plane surface of a conducting half-space. Expressions for the fields are 
obtained using three different approaches. The formulas developed are quite 
simple and, taken together, the whole range of distances from the far-zone to 
the near-zone is adequately covered. 


INTRODUCTION 


It is well known that a horizontal wire antenna will radiate horizontally 
polarized electromagnetic waves along the surface of the earth. It is not quite 
so well known that such an antenna will also radiate vertically polarized 
radiation. In fact the statement is sometimes made that the field of a horizontal 
electric dipole in the direction of maximum field is proportional to the field 
of a vertical magnetic dipole. Actually, at low radio frequencies, the field 
radiated from a low horizontal antenna is predominantly vertically polarized 
and is proportional to the field of a vertical electric dipole. 

At very low frequencies (less than 30 kc/s) it would be expected that the 
radiated field of a horizontal antenna is almost purely polarized in the vertical 
direction. However, if the range is comparable to the wavelength it is not 
immediately obvious that this situation prevails. It is the purpose of this 
paper to investigate this general problem. Expressions are derived for the 
electromagnetic fields of a horizontal electric dipole which is over, on, or 
below the surface of a homogeneous conducting half-space. 

To facilitate discussion of the past work in this field three overlapping 
distance ranges are considered. These are: 

(i) p>Xo, 

(ii) p~ Xo, 

(iii) p<KXo, 
where p is the horizontal separation between transmitter and receiver and Ao 
is the free-space wavelength. 

Case (i), the radiation field, has been treated in detail by Sommerfeld 
(1926) and Norton (1937) for the case where the dipole and the observer are 
in the upper half-space (i.e., the air). The extension to a spherical earth has 
also been developed (Wait 1956). The radiation field of the horizontal dipole 
immersed in the half-space has been treated by Moore (1951) in a doctoral 
thesis and by Banés and Wesley (1953). The problem was also discussed by 
Wait (1953a) in a published note. 

Case (ii), the near field, has been considered by Norton (1937), who states 
that his formulas are valid for distances down to a wavelength. His results are 


1Manuscript received March 23, 1961. 
Contribution from the Central Radio Propagation Laboratory, National Bureau of Stan- 


dards, Boulder, Colo. 
Can. J. Phys. Vol. 39 (1961) 
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restricted to the upper half-space. This intermediate distance region has also 
been treated by Bands and Wesley (1953) for the dipole in the half-space. 

Case (iii), the quasi-static case, has been treated by Foster (1931) and 
Fock (1933) for the horizontal dipole on the surface and by Lien (1953) and 
Wait (1959) for the dipole inside the conducting half-space. 

In the present paper simple field expressions are obtained for the three 
distance ranges mentioned above. No attempt was made to adapt specialized 
results developed by others; rather, self-contained derivations are given for 
each case. However, some of the methods are very similar to those employed 
by the author in earlier papers treating the fields of a magnetic dipole in and 
over a conducting half-space (Wait 1951, 1952, 19536, 1954, 1959; and Wait 
and Campbell 1953). 


FORMAL SOLUTION 


The earth is assumed to be a homogeneous medium with conductivity o, 
dielectric constant ¢;. With reference to a Cartesian co-ordinate system 
(x,y,z) it occupies the half-space z < 0. The air with dielectric constant eo 
occupies the space z > 0. The magnetic permeability of the whole space is 
assumed to be constant and equal to yw. The electric dipole, of infinitesimal 
length ds, is situated at (0,0,4) and is oriented in the x direction. 

The formal solution of this problem was given by Sommerfeld (1926). The 
method he used is now considered to be straightforward. Expressions for the 
components of the Hertz vector are derived by matching tangential fields at 
the interface z = 0. He showed that for z > 0 and hf = O, 





II e100 e ort roo —uo (2 +h) 

' Be ee af Sa Sure, 

(1) Co Ro R, 0 Ut o(Xe) 
II, = 0, 
I ss A "eee ee 

4 — = — 2 — | re et eee a hi r r dy, 

(2) Co Ox Jo Yolityitto oe) 

where 


Ro = [(g—h)? +p}, Ri = [(2+h)*+e"3, pp = (x*+y°)}, 
Yo = —equw, Yi = toyyw—eyw', 
uo = +5, ui = N+}, and Cy = Ids/(Amwiew). 


Time enters according to the factor exp(iwt). 
The fields are then to be obtained from: 


(3) E = —yilt+grad div Ml, 
(4) H 


dey w curl II. 


When z < Oand hk = O it is found that 
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Il, eo ete ue 
5 = 2 
(5) C. J aaa Jo(Ap)d dx, 
Il, = 0, 
; i 0 J ° (wo~ me mF 
et ne ee 
(6) ci Ox Veli tye “ Ja(o)d an, 


where C, = Ids/[4a(o1+7€.w)]. 
FAR-ZONE FIELDS 
In the present work it is assumed consistently that |y1| > |vo|. This situation 
is quite typical for low radio frequencies. Exceptions may occur when the 
ground is frozen. Because of the largeness of yj it is permissible to replace 
“1 by y1. Thus 





Il, e 1oRo e oF co e uotetn) 
7 eer ones ereear ee 2 
(7) C= Rh R, +2 f en Jo(Xp)d dd 
and 
Il, nN 2a [ae a oe, = —uo(z+h) 
(8) Co ™ ¥1 Ox uotyo/ 1 Uotyn P er 


These integrals are of a type which has already been evaluated by Sommerfeld 
(1926). They may also be treated by a modified saddle-point method (Wait 
1957). Therefore, if kp > 1 


Ts e V0Ro eo tee eR 
eg See ss SS 
where 
p’ = —(ikR,/2)(y1/10)", 
, sth ve) 
= 1 
™ »( “ee Ri ¥1 
and 
oy ae 1 
Ff (w ) —= "Dag . 
Similarly, 
em 
(10) aS [1 i(ap)he™ erie iw!) 
“0 1 
where 
= — (tkRi/2) (yo/v1) 
and 


z+h u)' 
w = (1 1+ a ar 
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It is immediately evident that the fields derived from II, are very similar 
to those of a vertical magnetic dipole whereas the fields derived from II, are 
related to those of a vertical electric dipole. This is particularly evident when 
these expressions for II, and II, are written in the form 


IL, 


(11) G = 2F(p’) (1 +12) (I +yih)e 7" /p 
~ [2 yo/ (vip) ] +r) (I+yhje™, 
and 
i. . 2 —Yop 
(12) C, = —[2 vo/ (yw) FP) (cos ¢)(1+-yv02/71) (1 +-yoh/yi)e™, 


which are valid for small heights (i.e., 22 and h? < p?). It is seen that the 
factors which depend on height have the form (1+7kAz) in which A = y1/yo or 
A= 0/71 

Using equations (11) and (12) it is a relatively straightforward matter to 
obtain the following expressions for the fields: 


! ; 
mw 7, —, uwlds Yo 


(13) E,=> a F(p)e*°” (cos ¢), forz = h = 0, 

(14) Ezy anates Yo 9-1 (sin $), forz =h=0, 
Vi Tr, 1 

(15) E,= inode as F(p)e**” (cos )Gy, for zandh = 0, 
2mp v1. 

(16) H,= a, Yo e ’° (sin @)Gu, forzand h 2 0, 
“7p Y1 

where F(p) = 1—i(mp)*e-erfc(ip?), 

Gy = (1+-y02/11)(1+70h/71), 
and Gy = (A+yi2z)(1+y14). 


Equations (13) to (16) are valid subject to the following approximations: 
(i) |yvop| = kp > 1 (i.e., far-field), 
(ii) A? and 2 « p? (i.e., small heights), 
(iii) [yi] > {vol (i.e., highly conducting ground). 





NEAR-ZONE FIELDS 
The use of the saddle-point method to reduce the exact integral forms for 
II, and II, to the forms given by equations (9) and (10) requires that kp > 1. 
When kp becomes comparable with unity it is necessary to use a different tack. 
In the case h = z = 0 it is shown in the Appendix that 


IT, 2 = = 
as [(1-+-v0p)e 7” — (1+ yap)". 


17 = ae ay 
(17) Co (vYi-—‘Yo)p 


WAIT: ELECTROMAGNETIC FIELDS 1021 


This is an exact result. Now under the restriction that |y1p| > 1, this simplifies 
to 


€ 


? 
(18) 3 (1+-yop)e***. 


For |yop| > 1, this agrees with equation (11). Furthermore, since |yi| > |v6| 
it is reasonable to assume that the factor yiwo+you1 in the exact integral 
form for II, can be replaced by yji#o for the near zone. Thus 


II, ap 2 oO . —1 —uo(z+h) 
sy ox Jo Uo e Jo(Ap)dr 


19 ! 
(19) C, 


9 
a (¢*" /R;). 
Y1 Ox 
For z = h = 0, this becomes 


Il ° ai 
(20) C= 7 (2/y1p )(1+vop)e””’ cos ¢. 


When |yop| > 1 and provided p <K 1, this agrees with equation (12). 


The corresponding forms for II, and II, when z and h > 0 are 
‘ I, ~ (9 /n2,3 —Y0P 
(21) G= (2/vip') (1+-yop)e *’Gu, 


Il, : cai 
C. > —(2/y:p )(1+yop)e "”’ (cos $)Gy. 


The resulting field expressions are then given by 
Lw iuw Ids 


(23) E,= a A=: Oapiy: 2 (1+-yee+vip" ye”? (cos $) forz =h 
1 2 


TH By, MOTE (1 +-yep)e™ (sin 6), for s = h 
v1 ™p'y 


tuw Ids 
Qap 71 
Ids 


& 504 5 (34+3yo0p+vyip Je”? (sin ¢)Gu, for zandh = 


(1+yop)e ’”’ (cos ¢)Gy, for z and h 


These expressions are valid subject to the following approximations: 


(i) |p| <1, (i.e., small ‘numerical’ distances) 
(ii) A? and 2 < p’, 
(iii) ly > vol. 
It is seen that the range of validity of equations (13) to (16) and of equations 
(23) to (26) overlap when, simultaneously, |p| K 1 and |yop| > 1. 
The preceding discussion has been restricted to the case where the observer 
is in the upper half-space. The analysis for the situation where the observer 
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is in the lower half-space (i.e., z <0) is very similar. The exact integral 
expressions for II, and II, then contain a factor exp (wz) in the integrand. 
If |y:| > |yo| and |y:p| > 1, this factor may be approximated by exp (712) 
and thus taken outside the integral. The resulting expressions for the fields 
are identical with equations (23) to (26) for the field components E,, Hg, 
E;, H,, and H, but now Gy and Gy are defined by 


(27) Gy & (1+-y0h/11) exp(v12), 
(28) Gy = (1+y1h) exp (712), 


and (23) and (24) have a factor exp(y1z) on the right-hand sides. The £, field 
component is not continuous at the interface. In fact, 

E,(2 + 0—) & (y0/vi)E.(2 + 0+), 
which indicates that equations (15) and (25) for EZ, are applicable to z < 0 
if Gy is replaced by Gy where 
(29) Gy = (1+-yoh/71) (v0/7i) exp(v18). 

The other obvious extension is to the case where the horizontal dipole is 
lowered into the conducting half-space (i.e., h < 0). A factor exp (uh) is now 
contained within the integrand of the exact integral expressions. This can be 
approximated by exp (712) subject again to the limitations that |y1| > 1 
and || > |yo|. For z > 0, the resulting expressions for the fields are again 
given by equations (13) to (16) and (23) to (26) if 


(30) Gry= exp(yih) (1+y02/71), 
(31) Gy = exp (y1h) (1 +712), 


and (23) and (24) for z = 0 are now multiplied by exp (71h). 
For 2<0 


(32) Gy & exp (y1h) exp (y12) = Gy 
except for the E, component which requires that 


(33) Gy & (75/73) exp(y1h) exp(y12). 


Furthermore, (23) and (24) are still valid provided the right-hand sides are 
multiplied by exp (yh) exp (712). 


QUASI-STATIC APPROACH 
At very short distances where the condition |y:p| > 1 becomes violated 
it is necessary to use a quasi-static treatment. If |yop| «1 it is permissible 
to regard the fields in the upper insulating half-space as solutions of Laplace's 
equation. This is equivalent to setting yo equal to zero in the exact integral 
expressions for IJ, and II,. For 2 < 0 and A S O these may be written 


EME PR I EE EA OA TNE A Se 


| 
| 
| 
: 
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ee) 1s u(z+h) 
(34) II, = cos | -42 Sa o Jo(dp)d an | 
and 
lisa pre 
2m(c+we) dx Jo “utr 
where 7? = yi = ipw(o+iew) and u = u, = (A2+7?)4. (The subscript 1 is 


dropped for the sake of brevity in the remainder of this section.) 
As z and h approach zero (from negative values) it easily follows, from 


equation (17), that 


(35) i, = Jo(Xp) dd 


Ids Si 
Reo te oth ™ 
Fulo-+ios) yp 3 [1—(1+yp)e"’]. 
The integral involved in II, is not quite so simple. However, it may be reduced 


as follows 





(36) I, = 


= P Jor) 4 am r=4f° a 
(37) Jo ek 7 (u%—d)Jo(Ap) dX. 


Since 
J AJo(Ap)dr = 0 


and using 


ue? £ +h *) 70000) = 0, 
p Op 


it is seen that 


ps of hie > ( 2 19 weir Todo) 
(38) 7 J aT ee ly See et a. 


According to Watson (1944) the integral on the right-hand side is equal to 
To(yp/2)Ko(yp/2) where Jo and Ko are penny Bessel functions of the first 
and iAeitt kind, both with argument yp/2. Thus, for s = h = 0— 


j Idscos@ @ ( 2 10 é) 
e y , = : ae f 0 ~ 0 2) . 
39) " 2a(ot+tew)y dp ’ p ap? dp [Zo(yp/2)K o(yp 


In a similar fashion it is found that, for s = h = 0—, 


aul Ids (! aa (: aa :) : 
40 = - : 5 > . ——y } [Lo 2)Ko 2) 
(49) Oz 2r(a+tew)y \p Op e Op/ \p Op ? Op 7 ; [Zo(ve \ve 


and 


’ lds cos @ 
(41) divi = elo-biew)e™ 
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The field components in the interface are then obtained from 


(42a) E, = —vy'Il, cos b+5, div II, 
(42b) E, = 7'Il, sin Pe cA div I, 
p 0g 

; Ol, .; Olle. 
A? = eek ig aoe 
(42c) H, = (o+tew) | 2 ag on 6| , 

i ale 
(43a) Hs = (o+%ew) | az ©O8 o ap | P 
(43d): H, = —(ottew) sin ¢. 


Within the limits of quasi-static theory E, vanishes inside the conducting 
medium for h = 0 and gz tending to zero from negative values (i.e., z ~ 0—). 
On the other hand, if kh = 0 and gz tends to zero from the positive values 
(i.e., z + 0+) it is found that 


‘ _tuwlds 0 "© Jo(Ap) 
(44) E, 24 OxXTo UtrA 


_ _inwlds a 21a | (22) (22) | 
(45) = one et Y pa? ap I 9 Ko 2 : 


In the limiting case where |yp| > 1 it is seen that 


dr 





(46) E, = > cos ¢. 


The field components E,, Ey, and Hy in this quasi-static approximation 
p ¢ ¢ 
may then be written concisely in the form 








; ee ill i 
(47) E, a 2a (o+iew)p° 3 [1+ (1+vyp)e ] cos ¢g; 
: Si ce a nee 
(49) Hy = ses B— (B+3y0+7' eM sin g. 


These results are valid for z = h = 0 and |yop| «1. The value of yp is 
unrestricted. In particular, if |yp| approaches zero the above reduce to 





aia ies Ids 
(50) E, = wae cos 9, 
" : Ids ; 
(51) ho= A 3 sIn @¢, 
(52) n= tds sin ¢. 

47p 


These elementary forms are consistent with potential theory. 
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Another interesting check is to note that if, simultaneously, |yop| < 1 and 
|yp| >> 1, equations (47), (48), and (49) agree with (23), (24), and (26) for 
z=h=0. 

Apparently the field components H, and Hy, for z = h = 0 cannot be 
written in a concise form. The higher derivatives of J>K 9 can be reduced to 
a linear combination of the functions [9Ko, JoKi, 11Ko, and J;K,; where J, 
and K, are modified Bessel functions of order unity. These forms will not be 
written out. However, it is of interest to note that if |yp| >> 1 each of these 
products can be replaced by 1/(yp) where the error is of order 1/(yp)*. To 
this approximation 

Idssin¢@ . 


(53) FH, >- Typ = inca Es, 


(54) Hye Sate Pe, 
2rvp 1p 

The previous discussion of the quasi-static case was restricted to the case 

z = h = 0. The manner in which the fields vary with 2 or A in this range 

(i.e., |yop| K 1) is quite complicated. For example, the fields do not generally 

vary in a simple exponential manner with depth into the conductor. Never- 

theless, for the case z and h < 0, it is possible to express the fields in terms of 


modified Bessel functions. Setting yo = 0 it is found that 


- —yRo 2 3 
(55) i Ids fs ps2 (2 P_ an, an) 


~ 4e(otiew)L Ro y\ds az az 


and 


3 2 
(56) Ids (#8 =) 


= 24 (o+7€w) axdz’ dxdz 


where 


Poo guleth) 
(57) N J S Jo(Ap)dr, 
0 


u 


cog (eth) et™ 

r a sips: acca 
(58) I J > Jo(Ap)A dd R° 
The form of equations (55) and (56) may be readily verified by carrying out 
the indicating differentiations and comparing the results with equations (34) 
and (35). 

The integral N is known. From Foster (1931) or Magnus and Oberhettinger 
(1954) 


(59) N = Iol(y/2)[Rit(2+h)]]Kol(y/2)[Ri— (2+h)]]. 


This result can be checked by noting that equation (59) reduces to the integral 
in equation (38) when z+/ = 0, and by seeing that MN satisfies the wave 
equation (V?—~?)N = 0 when zg and hk < 0. 
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The electric field in the conducting half-space may now be found from 
E = —7’Il+grad div Il. 


Expressing these in Cartesian co-ordinates it is found, after some simplification, 
that 


Ids 4 oT ye ) {2p i 
0) Se atal- (245 Ro ‘7-4 ayant oz , 
ie aa ee >) 9 oN i. 
(61) i se Ee Ro rye Oxdydz 


; Ids ? ie p) |. 
2 4 =- secnanbedtis | antennas 
(62) £, 4 (o+iew) ‘e 2 f Ro or 


As a check on these results it may be noted that as |yRo| and |yR:| — 0, 
the electric field components can be expressed as gradients of a scalar potential. 
Thus 


ee SARS 
(63) E = 4a(ot+tew) grad] 2 erehy P 


This result is consistent with potential theory. 

Another interesting limiting case is when |yp| >> 1 and z+ < p. Under these 
conditions the Bessel functions J) and Ko may be replaced by the first terms in 
their respective asymptotic expansions. This leads to the asymptotic 
approximations 


(64) N = (1/yp) exp (y(z+h)) 
and ON/dz & (1/p) exp (y(z+h)). 











The field components are then given by 


” bas i eee @), 
(65) E,= 2a(o+iew) , ay p 
Ids Y fh 

.f at aa ae aa AACR) S) ; 
a) Ey = or (o-+iew) , dxdy \p 


Expressing these in cylindrical co-ordinates 


= > ~ _Lds cos @ yet) 
(67) iss = Pe é : 





a ee Ids sin e evtetn) 
) y= t(o+iew)p* 






These are in agreement with equations (23) and (24) under the conditions that 
lyop| K 1, |yp| > 1, andzandhA SO. 

The manner in which the exponential factor exp (y(z+h)) occurs is rather 
interesting. It is only in the integral N that this factor emerges. The integral 
P, within the approximation that |yp| > 1, gives a negligible contribution. 
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This exponential factor is usually obtained by somewhat heuristic arguments. 
In the quasi-static treatment it is established on a sound basis and the limits 
of its validity are given. 


ACKNOWLEDGMENTS 
I would like to thank Mrs. Eileen Brackett for her assistance in preparing 
the manuscript and Mr. Kenneth Spies for his careful reading of the paper. 
I also wish to thank F. J. Zucker for suggesting the problem. The research 
was sponsored by the U.S. Air Force Cambridge Research Laboratories under 
contract A58-40. 


REFERENCES 
Banos, A. and WEsLEy, J. P. 1953. Report from Marine Physical Laboratory, University 
of California. 
Fock, V. A. 1933. Ann. phys. 17, 401. 
Foster, R. M. 1931. Bell System Tech. J. 10, 408. 
Lien, R.H. 1953. J. Appl. Phys. 24, 1. 
Maanus, W. and OBERHETTINGER, F. 1954. Functions of mathematical physics (Chelsea 
Publ. Co. New York), p. 133. 
Moore, R. K. 1951. Ph.D. Thesis, Cornell University. 
Norton, K. A. 1937. Proc. I.R.E. 25, 1203. 
SOMMERFELD, A. N. 1926. Ann. Physik, 81, 1135. 
Walt, J R. 1951. Can. J. Phys. 29, 577. 
1952. J. Appl. Phys. 23, 497. 
1953a. J. Appl. Phys. 24, 958 (in oqrasen _ V should be replaced by U). 
1953b. Trans. Am.‘Geophys. Union, 34, 1 
1954. Can. J. Technol. 32, 1. 
1956. Can. J. Phys. 34, 586. 
1957. J. Research Natl. Bur. Standards, 59, 365. 
1959. Appl. Sci. Research, B8, 213. 
Ww AIT, r, J. R. and CAMPBELL, L. L. 1953. J. Geophys. Research, 58, 167. 
Watson, G. N. 1944. Bessel functions, 2nd ed. (Cambridge University Press). 


APPENDIX 
It is of interest to observe that the vertical component of the magnetic field 


may be expressed in closed form. The first step is to note that IT, evaluated at 
h = O and z = 0+ is given by 


Me _ 5 (Joe) , 
Co 0 rt 


Multiplying numerator and denominator by u1—1o it follows that 


(69) 


(70) TOF sisi | J da Jo(Ap) dd — ~{" puado(ne)ar |. 
Co i> 0 


From Watson (1944) 


(71) £ ui rJo(Ap) dd = pte” 
0 


and using 
(72) (82 2433 ‘) F09) = 


p op” dp 
it follows that 
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3) ee 2- ( 21a - shee _ io 
Co ¥vi-Yo - p op” ap p , p op” dp p 


ae 
= Gioyiypt [tree (L+-yop)e 





Then 
a)” Pie Beets 
Op 


= tas sin 9 (3+ 3yep+ rie" e™ — (8-437 0+7i6)e 
2m (vi—Yo)p 
This is an exact result valid when z = 4 = 0. In the limiting case for 
|y1p| > 1, it agrees with equation (26), and if |yop| « 1, it agrees with equation 
(49). It should be noted that strictly the integrals in equations (70) and (71) 
are divergent. The operations performed, using these relations, can be made 


rigorous by simply replacing 


yr .iibyte Oe 


a50 0 


ON MEASUREMENT OF THERMAL CONDUCTIVITY 
AT HIGH TEMPERATURES! 


M. J. Lausitz 


ABSTRACT 


A method is given for exact mathematical analysis of linear heat flow systems 
used in measuring thermal conductivity at high temperatures. It is shown that 
a popular version of such a system is very sensitive to the alignment of its 
components, which seriously limits the temperature range of its satisfactory 
use. 


INTRODUCTION 


For the past several years great interest has been shown in the thermal 
conductivity of solids at high temperature. Diverse methods for the deter- 
mination of thermal conductivity and the results obtained for various solids 
have been reported in the literature. However, the agreement among the 
experimental results, particularly for ceramics, has not been spectacular. The 
variation in values obtained for the thermal conductivity of ceramics has been 
generally ascribed to the nature of the specimen and its microstructure and 
impurities. The variation in the structure and composition of the specimen 
is undoubtedly responsible for a large part of the disagreement of the reported 
experimental results, but another factor, the straightforward uncertainty of 
the experimental determination, also plays an important role, usually grossly 
underestimated. The exaggerated accuracy sometimes reported for values of 
thermal conductivity at high temperatures is the result of an excessive con- 
fidence in the accuracy of thermocouples, and ignorance of the heat flow and 
temperature patterns developed in the experimental arrangements. Most 
experiments, particularly the steady-state ones, aim at the establishment of 
a simple heat flow pattern. The actual heat flow patterns do, of course, deviate 
from those assumed, but such deviations are usually considered negligible. 
This article is written to show the importance of those neglected deviations 
of the heat flow patterns, and to stress the necessity of precise mathematical 
analyses of the experimental arrangements used for measuring thermal 
conductivity. 

The steady-state methods of measuring thermal conductivity employ either 
of two heat flow patterns: radial or linear (the ellipsoidal method of Adams 
(1954) is essentially a variation of the radial heat flow method). Of the two, 
the linear method has lately been the more popular and will be discussed 
mathematically in the following pages, where first there will be given a 
general outline of the mathematical approach, and then a detailed analysis 
of one particular experimental arrangement. 

‘Manuscript received March 3, 1961. 
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MATHEMATICAL ANALYSIS OF THE LINEAR HEAT FLOW METHOD 
OF MEASURING THERMAL CONDUCTIVITY 
A general arrangement of a linear heat flow apparatus is shown in Fig. 1. 
(A) is the specimen of radius a and conductivity k,, (D) is the guard, of 
thickness b—d and conductivity ky. The space between the guard and speci- 
men (B) is filled with an insulating substance of conductivity k2, and outside 
the guard (C) there is insulation of conductivity k;. The outside diameter 
of the apparatus is 2c, and its length L. There are several heaters wound on 





> 2A (a) 





+2 rm(b)—~ 





2h, (b)— 





r=c r=b r=d r=a or 


Fic. 1. A general arrangement of a linear heat flow apparatus. 


the guard: the mth heater, of length 2),,(b) and positioned at z,,(b), dissipating 
power Q,,(b) per unit length. There are also heaters on the specimen. The /th 
of these, of length 2\,(a), is positioned at z,(a) and dissipates power Q,(a) 
per unit length. There are 7 thermocouples located at both the guard and 
the sample, at positions 2;. 

A few conditions must be assumed for a complete mathematical analysis 
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of the system. Primarily, all the conductivities are assumed independent of 
temperature. This is not a severe limitation, as indicated by the very good 
agreement between calculated and observed temperature distributions shown 
in Fig. 2. Then, radial symmetry is assumed, and to make calculation easier, 
the temperatures along z = 0 and z = L, and along r = ¢ are all assumed 
zero. These again are not very critical assumptions, as the surface tempera- 
tures have little effect on temperature and heat flow distributions in the 
inside of the furnace (see, for instance, Laubitz 1959). 

Under these conditions, the temperature in region X (X stands for A, B, 
C, or D) are given by, 


T(r,2) = QO [Fa(X)Io(anr) +G,(X)Ko(anr)] sin ag2 


where J) and Ko are the associated Bessel functions, an = nw/L, and F,(X) 
and G,(X) are constants for each region X. G,(A) is, of course, zero. The 
constants F,(X) and G,(X) are so chosen as to satisfy the continuity con- 
ditions at the interfaces of the regions A, B, C, and D. The temperatures of 
interest are those along r = 6 and r = a, and it can readily be shown that 


(1) Te,s)- > Qua) | (x) sin a, Xala) even a) dines 


+ LD See gale) . sim aghn(B) sin oy () sin ans 
m - 3 


n 


where x is a or b, and f,(x) and g,(x) are functions of the ratios of conductivities 
and the parameters of the apparatus. These functions are too complicated to 
be given here, but their explicit formulae are given in Laubitz (1961). The 
above equation shows that the temperatures at any point along a or b are 
given by equations linear in the power inputs. Once the coefficients of these 
equations are calculated, it is then a simple matter to study the effects of 
lateral temperature match or mismatch, of relative displacements of specimen 
and guard, and other possible conditions affecting the behavior of the apparatus. 

To check the stringency of the assumptions made in the derivation of the 
formulae, calculated and observed temperature distributions were compared 
for a furnace available in this laboratory. This furnace was 32 cm long and 
had an outside diameter of 25cm. The guard, of dense alumina, had seven 
heaters and six thermocouples. The specimen, also of alumina, had two 
heaters and also six thermocouples. Regions B and C were filled with powdered 
alumina. Two comparisons were made: one with the guard maintained at a 
uniform temperature 7) = 1000°C, and no heat input into the sample 
heaters, and the other with only the sample heaters energized. The calculated 
and observed relative temperature distributions, in excellent agreement with 
each other, are shown in Fig. 2. In the calculations, the summation of the 
infinite series was carried up to 12 terms. 
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Fic. 2. Comparison of observed and calculated relative temperature distributions: A, 
uniform guard temperature of 7», no power into sample heaters; B, no power into guard 
heaters. Lines are calculated temperatures, and dots observed temperatures. 


DETAILED DESCRIPTION OF A PARTICULAR FURNACE 

The usefulness of the mathematical analysis in the measurement of thermal 
conductivity using linear flow lies in the fact that accurate formulae can be 
obtained for two of the major sources of experimental error, lateral tempera- 
ture mismatch and relative specimen-—guard displacement. Unfortunately, the 
general formulae for these effects are too unwieldy to show the relationship 
between the pertinent factors; each apparatus has to be analyzed separately. 
As an example, one analysis is given here in some detail. It is pertinent to a 
design very similar to one originally described by Francl and Kingery (1954), 
with five heaters on the guard, two on the specimen, and four thermocouples 
located on each of the guard and specimen. In the calculation, the parameters 
assumed, given as ratios to the length L, were as follows: 


a = 0.058 2, = 0:326 

b = 0.174 Z2. = 0.442 

c = 0.465 23 = 0.558 

d=0101 24 = 0.674 
2,(b) = 2,(a) = 0.256 di(b) = A5(b) = 0.070 
zo(b) = 0.384 ho(b) = A3(b) = Aa(b) = 0.047 
23;(b) = 0.500 Ai(a) = dA2(a) = 0.064 


z4(b) = 0.616 
2,(b) = 2.(a) = 0.744 


The conductivities of guard and specimen were assumed equal (k;), and 
so were those of regions B and C (k.). 

Although the factors f,(x*) and g,(x) of equation (1) decrease as n~*, the 
convergence of the series is spoiled by the sines. To find the minimum number 
of terms required to give good results, the temperature distribution along 
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the specimen was calculated as a function of the number of terms summed 
for ki/ke = © and Q,(a) = Q2(a). The results are shown in Fig. 3, as devia- 
tions from the linear interpolation of temperatures between points 2, and zy. 
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Fic. 3. Calculated deviations from linearity in the temperature distribution along the 
specimen for ki/ke = ©, as a function of the number of terms in the summation: —— 6 
terms; — — — 12 terms; --- 16 terms; —+—-~- 24 terms. 


Figure 3 is typical of an approach made by a Fourier series to the function 
it represents: the series quickly reaches a level within about one per cent of 
the function, and then slowly settles down to its final shape. In this parti- 
cular case, summation over 12 terms gives temperatures well within one per 
cent of the ultimate ones, which is sufficient for the study of errors. It is 
important to note, however, that the ultimate temperature distribution 

.deviates from linearity, best shown by the effect of summation on the tem- 
perature of the specimen mid-point. 

As all calculations of thermal conductivity in such a system are based on 
assumption of linearity of temperature, any deviation therefrom is a source 
of error. There are two aspects to the deviation: its magnitude, which is con- 
veniently given by the magnitude of the deviation at the mid-point of the 
specimen relative to the specimen temperature, 7; and its “symmetry”, 

which can be described by the deviation of the ratio R from unity, where 


2(7 (a, 2.)— T (a, Z2)] 


R==— — — — ae 
? 7 (a, 2:)— 7 (a, 2) +7 (a, 3) —T (a, 24) 
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If the apparatus is used for comparative determination of thermal conductivity, 
R is the ratio of the conductivity of the unknown to that of the standard. For 
a system symmetrical in the z-direction, R is only a function of the ratio 
of the conductivities, k:/k2. The relative temperature deviation from linearity 
of the specimen mid-point depends not only on k;/k2, but also on the over-all 
specimen temperature gradient, AT = T(a, 2:)—T (a, 2s). Figure 4 shows both 
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Fic. 4. Relative deviation from linearity of the specimen mid-point temperature (curve 
A) and deviation from unity of R (curve B) as a function of k/ke. 


these deviations as a function of ki/k2, for AT/T>) = 0.05. It can readily be 
appreciated that the deviation from linearity of the mid-point temperature 
is an error in absolute determination of thermal conductivity, while the 
deviation of R from unity is an error in comparative measurements. 

The errors mentioned above occur in a perfectly matched system, that is 
one where 7'(b, z,) = T(a, z;). Lateral specimen-guarded temperature mismatch 
is another source of error. Its treatment within this analysis is elementary. If 
we define e; = 7(a, 2;)—T (0, z;), the algebraic mismatch of the ith set of 
thermocouples, then temperature deviation along the sample caused by these 
mismatches will be sums of the e; multiplied by coefficients which, for any 
particular apparatus, are only functions of k;/k, As an example, the error 
in R due to mismatch, E, for the particular system here considered is 

KE = — 0:98 (eres) + 1.94 (€2—es) pit 
a AT 1 2 ’ 

0.17 (€1— e4) + 0.57 (¢2—es) 
7 AT 
__ 0.03 (e1—e4) +0.09 (e2—e2) 


= AT es ki /k2 = 100. 




















ki /k» = 10, 
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The errors e; can arise from random mismatch, inaccurate thermocouple 
location or reading, or from relative specimen-guard displacements. 

An interesting feature of the apparatus here considered, and one that does 
not seem to have been appreciated heretofore, is that it contains more arbi- 
trary variables than restraining conditions. The variables are the seven 
independently variable heaters; the restraining conditions are the specimen 
temperature 7, the gradient (AT) along it, and the four lateral temperature 
matches, only six in all. There is a further seventh condition, but it is weak: 
power inputs cannot be less than zero. Within that last condition, there is a 
wide range of heater settings that will satisfy the previous six conditions. This 
in itself has certain advantages, as it facilitates the matching of the guard 
and specimen temperature. Unfortunately, it also results in an arbitrariness 
in the temperature distribution of the guard and specimen, which greatly 
increases the uncertainty of the experimental results. As an example, Fig. 5 
shows the relative deviation of temperature from linearity for the guard and 
specimen for various permissible heater settings. The limits on the settings 
are given by Q3(b) = 0, and Q4(b) = 0, and from Fig. 5 it is seen that, although 
the specimen temperatures are not seriously affected, the guard temperatures 
vary extensively. The variation in guard temperature is of no consequence 
if the guard and specimen are perfectly aligned. If, however, there is a relative 
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Fic. 5. Relative deviation from linearity of the temperature of the guard (curve A) and 
of the sample (curve B) for various allowable heater settings: 
—— Q3(b) = 0; +++ Q4(b) = 0; —--+ Qs(b) = Q4(d). 
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displacement between specimen and guard in the z-direction, temperature 
mismatches will occur the magnitude of which will depend on the heater 
settings. These will in turn result in errors in the experimental determination 
of conductivity. In the apparatus here considered, maximum error for any 
given specimen-—guard displacement will occur when either Q3(b) or Q4(b) equals 
zero. It turns out that to a first approximation, the error in R due to a relative 
displacement of guard and specimen of 6z, E(éz), is the same in magnitude for 
Q;(b) = 0 as for Q4(b) = 0, but of different sign. Its magnitude is given by, 
To 


|E (62) max| a Eo(ki/k2) ° AT Ag 


where Az = 2,—2n-1, is the thermocouple separation. 
E,(1) = 0.65; £,(10) = 0.138; and £,(100) = 0.01. 


For any particular heater setting, the error in R will be somewhere between 
+E(62)max and —E(6z)max, depending on the ratio of Q3(b) to Q4(b). Figure 6 
shows the possible error in R, |E(6z)max|, as a function of temperature, for 
dense alumina in insulating brick, as used by Francl and Kingery (1954), 
and dense alumina in powdered alumina insulation. Two curves are given 
for each case, one for a constant AT = 100°C, and one for a constant 
AT/T>) = 0.1. The relative displacement assumed, 6z/Az, was 1 mm in 25 mm, 
which is certainly very reasonable for a system constructed out of ceramics, 
unless very special precautions are taken during the lining up of specimen 
and guard. The errors shown in Fig. 6 are certainly not negligible at higher 
temperatures. They are, of course, the maximum possible due to the mis- 
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Fic. 6. Error in the comparative method of measuring thermal conductivity due to a 
misalignment of specimen and guard of 1 mm in 25 mm, for dense alumina in firebrick (curve 
A) and dense alumina in powdered alumina (curve B). AT = 100°C; AT/T,=0.1. 
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alignment, and it is unlikely that in the actual experiment either Q3(b) or 
Q,(b) would be very close to zero. On the other hand, it should be remembered 
that dense alumina has a high thermal conductivity for a refractory substance, 
and for most other ceramics the maximum error would be much larger. 

A very important feature of the error arising from the relative specimen— 
guard displacement, and one which is not obvious, is that it affects the 
specimen temperature intervals 7 (a, z:)—T7 (a, 22) and T(a,23)—7 (a, 24) in 
precisely the same way. This means that in a comparative measurement, 
where the end sections of the specimen are the standards and the middle 
the unknown, serious errors can arise in R, the ratio of the temperature 
gradient on the unknown to that on the standard, without any difference 
between the temperature gradients of the top and bottom standards. Thus, 
the agreement between the temperature gradients of the standards is no 
indication whatsoever of the errors in the determination of the unknown, 
contrary to popular opinion. 

The error due to misalignment can be minimized by choice of heater 
settings. For each particular temperature, there is an optimum ratio of 
Q3(b) to Q4(b), for which the misalignment error is very small. This ratio can 
be calculated, and the experimental settings brought close to it. The reduction 
in uncertainty achieved this way can outweigh the reduction achieved by 
precision of alignment, which is difficult at high temperatures. 

It may be mentioned here that radial displacements of sample relative to 
the guard are almost of no consequence. For the system here considered, the 
maximum change in temperature due to a relative radial displacement of 
0.5 mm is 0.2% for k,/ke = 1; 0.04% for ki/ke = 10; and 0.003% for k,/k: 
= 100. Furthermore, this maximum temperature change occurs at the mid- 
point between thermocouples, and is much smaller at the thermocouple sites. 

A few remarks may be made here about the effects of the features in which 
the real systems differ from the mathematical one. By comparison of calcula- 
tion with observation, the effect of the temperature variation of the thermal 
conductivity has been shown to have negligible effect on the temperature 
distribution. The outside boundary conditions are also not very important; 
they may change the magnitude of the temperature or its linear gradients, 
but will not seriously affect the temperature deviations from linearity. Some 
systems, like that of Francl and Kingery (1954), differ in two further aspects: 
their specimen is square in cross section, and its central regions have a smaller 
transverse dimension than its ends on which the heaters are wound. Placing 
a square specimen in a round guard is very much like the proverbial square 
peg in a round hole; it defies analysis, if not reason, and only vague statements 
can be made about its effects on the temperature distribution. There will, of 
course, be a component added to the sample temperature with a fourfold 


angular symmetry of the form 


dS Anlalanret sin aps. 
n 


It is likely that this component will have little effect on the temperature 
distribution for perfect alignment. For vertical misalignment, the effect on 
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the temperature distribution will be larger than for a cylindrical system with 
specimen diameter equal to the side of the square, but probably smaller than 
for a cylindrical system with specimen diameter equal to the diagonal of the 
square. Radial misalignment has very complex results, which may possibly 
still be negligible for perfect vertical alignment. The extension of the transverse 
dimensions at the specimen heaters is also very difficult to analyze accurately. 
However, relaxation calculations have been made for various ratios of heater 
diameter to specimen diameter, 7y/a. These showed that increasing the ratio 
from unity to two has very small effects on the specimen temperature dis- 
tribution. For instance, for k;/k, = ©, it reduces the deviation from linearity 
of the temperature of the specimen mid-point from —0.016 to —0.011. Its 
effect on the guard temperature in the central region is negligible. 

As a final example, a sample estimation of the uncertainty in the deter- 
mination of conductivity in a Francl and Kingery type system is given. The 
temperature, 7°, is assumed to be 600° C, and the over-all temperature gradient, 
AT, 50° C. The thermocouple separation is 25 mm. The specimen and guard 
are made of alumina, and the insulation is brick, giving a k,/ke of about 
20. The uncertainty in the thermocouple separation is of the order of the 
thermocouple diameter, about 1%, and the uncertainty in the temperature 
measurement is assumed to be about 0.5%. The possible misalignment of 
guard and specimen is of the order of 1mm, and the guard and specimen 
temperatures are matched to 1° C. The uncertainty in R, the ratio of the 
conductivity of the unknown to that of the standard, is then as follows: 

uncertainty due to thermocouple separation 
uncertainty due to temperature measurement 
error due to deviation of temperature from linearity 

for perfect matching 1.5%, 
uncertainty due to matching error 1.5%, 
uncertainty due to possible displacement of guard 

to specimen 5%. 

As these errors and uncertainties are not random, except for the matching 
error, the total uncertainty is of the order of the sum of the individual ones, 
about 11%. To that has to be added the uncertainty in the conductivity of 
the standard. 

One can also estimate the expected difference between the temperature 
gradients of the top and bottom standard. These are: 

uncertainty due to thermocouple separation .5%, 


MOT 


uncertainty due to temperature measurement To. 
uncertainty due to matching error 0%, 
error due to deviation of temperature from linearity 
for perfect matching 8.0%. 
The last error always increases the gradient of the hotter standard with 
respect to the colder. One would then expect the hotter standard to have 
a gradient 5 to 10% larger than the colder one. This is roughly in agreement 
with the sample data published by Francl and Kingery (1954, p. 81). There, 
after due correction has been made for the temperature variation of con- 
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ductivity of the standard, the gradient of the hotter standard exceeds that 
of the colder one by 5% for the run with the Pb sample, and by 23% for the 
run with the MgO sample (if one assumes that the separation of the thermo- 
couples on the standards was identical). The test temperature was about 
200° C. 


CONCLUSION 

The disagreement among the results of thermal conductivity measurements 
reported in literature is large, particularly for ceramics. Part of the disagree- 
ment is real, due to variation of the specimen tested, and as such yields useful 
information about the substances studied. A significant part of the disagree- 
ment, however, is due to experimental error, consistently underestimated. A 
clear insight into the source and nature of the errors can be gained by exact 
mathematical analysis of the apparatus employed in the experiments. This 
is particularly important for errors which will affect the accuracy of the 
measurements without affecting their apparent precision. 

The mathematical analysis of realistic experimental systems is usually 
within the limits of a desk calculator. In this paper, the method of analysis 
for systems using linear heat flow is indicated, and one particular system, 
that of Francl and Kingery (1954), is discussed in detail. 

The analysis of the apparatus for comparative measurements of Francl and 
Kingery brings to light a number of interesting features. It shows that the 
specimen temperature deviates from linearity even for perfect matching of 
the guard and specimen temperature. It further shows that the system is 
very sensitive to the vertical alignment of guard and specimen, and that the 
agreement between the temperature gradients of the top and bottom standard 
gives no indication of the accuracy of the determination of conductivity of 
the unknown. 

Evaluation of the uncertainty in a system like that of Francl and Kingery 
shows that its temperature range is very limited if the accuracy is to be kept 
within 10%. For alumina in insulating brick the upper temperature at which 
the uncertainty does not exceed 10% is about 600° C, and this upper limiting 
temperature rapidly decreases for the specimen whose conductivity is less 
than that of alumina. 
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THE CRYSTAL STRUCTURE OF Bi;,Te.Se'! 


J. A. BLAND? AND S. J. BASINSKI 


ABSTRACT 


The crystal structure of BizTesSe has been determined by X-ray single 
crystal methods. The unit cell is trigonal, space group R3m with hexagonal cell 
edges a = 4.28+0.02 A, c = 29.86+0.10 A, and there are three formula units 


Ms Sy 
per unit cell. The structure is of the C33 type and a comparison is made with 
other C33 structures, Bi2Se3, BiesTe;, and BisTe2S. An important common 
feature of these compounds is that the atom on the center of symmetry has a 
relatively large contact distance with each of its six Bi neighbors. 


There is considerable interest in the relationship of the fundamental elec- 
tronic properties, activation energies, and charge carrier mobilities of semi- 
conducting compounds to the short-range order existing in those compounds 
(i.e., co-ordination numbers, bond lengths, and bond angles). BizTe2Se is a 
semiconductor belonging to a series of substances with the C33 type of 
structure, namely, BioSe;, BizTe;, BisTesS, and it therefore appeared useful 
to determine its structure accurately so that a detailed comparison of inter- 
atomic distances in the binary and ternary compounds could be made. 

Crystals with the chemical composition Bi,Te.Se were obtained from a 
homogeneous ingot prepared by melting the constituents together in a sealed 
tube. They were in the form of thin plates with a metallic luster and they 
could be easily deformed. A single crystal was obtained and accurately set 
on a Buerger camera with the aid of an optical goniometer. The usual X-ray 
method of setting the crystal proved to be difficult on account of the high 
density of reciprocal lattice nets parallel to (0001), the basal plane of the 
crystal, which resulted in higher layers being present on zero level photo- 
graphs. These pictures showed, however, that [0001] is a trigonal axis. Reflec- 
tions in the zero, second, and third reciprocal layers perpendicular to the [0110] 
axis of the hexagonal cell were recorded using MoKa radiation and the system- 
atic absences were found to be consistent with the trigonal space group R3m; 
this choice was supported by the observation that the distributions of the 
intensities of the reflections in the three layers were identical. Indices of the 
reflections were assigned according to a hexagonal cell with dimensions a = 4.28 
+0.02 A and c = 29.86+0.10 A. The density calculated on the basis of three 
formula units per hexagonal cell is 7.93 g cm~* and the density obtained by 
the displacement method on a small quantity of crystals was shown to be 
greater than 7.4 gcm~. 

The structure was determined using the 000/ reflections recorded on the 
zero level Buerger picture perpendicular to [0110]. The intensities were 
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estimated by visual comparison with a standard scale and corrected for 
absorption of X rays in the specimen and for the Lorentz-polarization factor. 
The structure of tetradymite, BizTe.S (Harker 1934), was used as a starting 
point in the analysis, and refinement was carried out using Fy and Fo—F, 
one-dimensional syntheses along the c-axis of the hexagonal cell. A typical 
Fy synthesis is given in Fig. 1. The positions of the 15 atoms in the hexagonal 
cell are 


(000, 332, 338)+. 
Bi+00v, 
Te+00u, 
Se 000, 


where v = 0.396,+0.015 and uw = 0.211;+0.045. The limits of error were 
estimated ' from the gradients at the atomic positions on the final difference 


Bi 


c 
ed 


0 v-I/3 \/3-u /6c 


Fic. 1. Bi,Te2Se. One-dimensional Fy synthesis using 000/ reflections. 


synthesis. In the calculation of structure factors, a temperature factor of 
0.7 A~® was assigned to each atom; values of Fy scaled for comparison with 
F., are given in Table I. 

Bi:Te2Se is isomorphous with tetradymite, BizTeS, and may be considered 
as a layer structure in which each atom is co-ordinated by six neighbors; Bi 
is linked to three Se atoms and to three Te atoms, and Se is linked to six 
Bi atoms. In Table II(a), interatomic distances are compared with corre- 
sponding distances in other structures of the series, BigSe; and Bi,T3; (Semi- 
letov 1956) and BisTesS (Harker 1934). Although the accuracy of these 
structures is not known and detailed comparison is thereby made difficult, a 
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TABLE I 


Comparison of Fo and F, values of 
000/ reflections for Biz:Te2Se 
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TABLE IlI(@) 
Interatomic distances (A) in C33-type compounds 


Normal 
we Se —_—_———— valence 
I II Ill IV radius 
Bi.Se; BioTe; BioTeS BisTe2Se sum 


Structure No. 





.55 
.68 
.88 
34 
74 


Bi-Seo 3.05 Bi-Tey 3.24 Bi-S 3.06 Bi-Se 3.10+0.03 
Bi-Se 2.99 Bi-Te, 3.04 Bi-Te, 3.13 Bi-Te 3.03+0.05 
Ser-Se, 3.31 Te:-Te: 3.73 Te:-Te: 3.70 Te-Te 3.65+0.06 


lrpent 


feature which emerges from the study is that in each compound, the atom 
at the origin has a relatively large contact distance with each of its six Bi 
neighbors. This is apparent in Table II(b) where the Bi-M distances, with 


TABLE I1(d) 


Comparison of interatomic distances (A), Bi-M, in C33-type 
compounds tabulated as Bi-Mo and Bi-M;, 


M at origin, M in general position, 
Mo M; 

3.06111 

Se 3.05! 3.10!V 2.99! 

Te 3.2411 3.0411 3.131 3.031V 


The superscripts refer to the four structures in Table II (a). 


M representing S, Se, and Te, are set out in two columns, the first showing 
distances involving M at the origin, Mo, and the second showing distances 
with M in the general position M,. The effect is most marked in BizTe; where 
Bi-Teo is 3.24 A, which is considerably greater than the largest Bi-Te, dis- 
tance of 3.13A in BixTe:S. Similarly the Bi-Sep distances of 3.05 A and 
3.10 A in Bi.Se; and BisTe2Se are probably significantly larger than the Bi-Se, 
distance of 2.99 A in Bi.Se;. It is also probable that the Bi-So distance of 
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3.06 A in BizTe2S is relatively large because it has the same value as the 
mean Bi-Se distance in the other structures; the covalent radii of S and Se 
are 1.04 A and 1.17 A (Wells 1950) so that normally the Bi-S distance should 
be shorter than the Bi-Se distance, Table II(a). The relatively large size of 
the atom at the origin is also clear from a scrutiny of BizTe2Se. In this structure 
the Bi-Seo distance is just significantly longer than the Bi-Te, distance; the 
inequality is in the opposite direction to that expected from the values of the 
covalent radius sums given in the final column of Table II (a). 


We wish to thank Dr. E. Mooser and Dr. W. B. Pearson of the National 
Research Council of Canada, Ottawa, for suggesting the problem and pro- 
viding the crystals. We also thank Mr. R. C. McAdam of the Department 
of Mines and Technical Surveys, Ottawa, for the chemical analyses. 
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WAVES IN A RAREFIED IONIZED GAS PROPAGATED 
TRANSVERSE TO AN EXTERNAL MAGNETIC FIELD! 







TomiyA WATANABE? 






ABSTRACT 


Waves being propagated in a rarefied and fully ionized gas and transverse to 
an external magnetic field have been studied, particularly hydromagnetic 
waves. Three modes of waves, in which the perturbed magnetic field is parallel 
to the external magnetic field, are found to be propagated. In a high-frequency 
limit, they tend to electromagnetic waves, electron sound waves, and ion sound 
waves. In the condition that the Alfvén velocity is greater than the ion sound 
velocity but smaller than the light velocity, the last mode tends to a hydro- 
magnetic wave in the low-frequency limit. The other two modes of waves can 
be propagated only at frequencies higher than the critical frequencies, both of 
which almost equal the electron plasma frequency. The condition that hydro- 
magnetic waves should be attenuated severely due to collisions between electrons 
and ions has been derived. 















1, INTRODUCTION 


Waves propagated in a uniform conducting fluid under an external uniform 
magnetic field have been studied by many authors. For instance, we can 
mention van de Hulst’s work (1951), which employed the Navier-Stokes 
equation but with an additional term representing the ampére force and the 
ohmic law of electric conduction in the case of moving media: 









(1.1) j/o = E+ =uAH, 






where the meanings of notations are as follows: E, electric field; H, magnetic 
field; j, density of electric currents; u, velocity of the fluid; and a, electric 
conductivity. He reaffirmed the occurrence of hydromagnetic waves discovered 
by Alfvén (1950), which in an extremely conducting (¢ — ©) and non-viscous 
fluid can be propagated along the external magnetic field with the so-called 
Alfvén velocity defined by 










Vv, - = aie 
A WV (4np)’ 
where JJ» is the intensity of the external magnetic field and p is the density of 
the fluid. One of van de Hulst’s important and new contributions is that 
waves can be propagated not only along the external magnetic field but also 
transverse to it, if the fluid is compressible. In that case, the propagation 
velocity V is given by 


(1.3) V = V(VaAtQ) 


where C, is the sound velocity of the fluid. 
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WATANABE: WAVES IN RAREFIED IONIZED GAS 


In many interesting cases, the medium in which hydromagnetic waves 
should appear is a gas, especially a rarefied gas, which is ionized fully or 
partially. Ionized gases are good conducting media, and thus, hydromagnetic 
waves may be expected to occur in them. In that case, however, it is not 
clear a priori that the law of electric conduction (1.1) would still hold. This is 
an experimental law concerning metallic substances. The equation of motion 
is also a problem. Generally, an ionized gas is in a so-called plasma state. 
Namely, it consists of electrons and ions, number densities of which are equal 
to each other, in order to avoid the dispersion due to the electrostatic repul- 
sion. Since the ion mass is much greater than the electron mass, the ion gas 
does not necessarily move with the electron gas. We cannot describe such 
a situation with only one equation of motion, that is, the Navier-Stokes 
equation with an additional term due to the ampére force, but we need two 
equations of motion, one for the ion gas and the other for the electron gas. 
These problems must be reviewed from a microscopic point of view. Astrém 
(1951), Schliiter (1952), and Ferraro (1955) treated the problem of waves 
propagated along the external magnetic field in this manner. Astrém also 
studied waves propagated oblique to the external magnetic field as well as 
normal to it. But, he did not take into account the gas pressure, which may 
be of importance in oblique and, especially, normal propagation. Piddington’s 
work (1955) took into consideration the pressure of the electron gas, but 
he omitted both motion and pressure of the ion gas. In this paper, we shall 
try to study the problem, taking both motion and pressure of the ion gas 
into consideration, but we will restrict the problem to the case of normal 


propagation. 


2. FUNDAMENTAL EQUATION 

As mentioned above, the first problem is to formulate the equations of 
motion and of electric conduction rigorously. In order to do so, it may be 
best to start from the Boltzmann equation. Let us consider waves propagating 
in a fully ionized gas, which consists of electrons and ions. The mass of an ion, 
m,, is assumed much heavier than that of an electron, m,. In an undisturbed 
state, that is, in a state where no waves appear, the number density of elec- 
trons, V., is equal to that of ions, Vy, and it is expressed by V,. The electric 
charge of an ion is represented by e (>0). Accordingly, the electronic charge 
is —e. We shall employ the Gaussian system of units of electromagnetic 
quantities throughout the paper. 

Under these conditions, one can derive the equation of the electron gas 
from the Boltzmann equation as follows (Spitzer 1956): 


(2.1) oo( S84 (u,: v)u,) = —N, (p+! ue A H) — Vet —peVO+Pai. 


Pe = Nm, is the mass density of the electron gas. The permeability is assumed 
equal to 1. @ is the potential of the external and non-electromagnetic force 
per unit mass such as the gravity force, which we shall omit in this paper. 
u, is the mean velocity of electrons in a volume element, AV, and is given by 
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a 
Ue = Fra a Ve 


where V,’s are velocities of individual electrons. The quantity t, is the stress 
tensor, the Ym component is given by 


(2.3) Vim = Ay Le Ve We) 1(Vo— We) m 


the summation extending over the volume element. 

In the simplest case that the distribution of the random velocities v.—U, 
is isotropic, Ym vanishes unless / equals m, and the three diagonal components 
Wor, Vyy, Vez are all equal to each other and to the scalar pressure p (Spitzer 
1956). In this situation, we have 


(2.4) Vi — Vp. 


We assume this simplification throughout the paper. According to Spitzer 
(1956), this replacement of the stress tensor by a scalar pressure means we 
ignore the damping in wave propagations discovered by Landau (1946), and 
thus, the following analysis should be taken only as a model case neglecting 
this kind of damping. 

The quantity P,, is the total momentum transferred to the electrons per 
unit volume per unit time by collisions with ions. Schliiter (1950, 1951) has 
assumed the following form for this quantity, 


(2.5) Fes Se N.Niae(Ue—U)) 
where 


> Vv mem, 
(2.6) Oy = 
N,im.+m, 
and wu, is the mean velocity of ions. v_ is the collision frequency of an electron 
with ions. 
Under the two assumptions concerning the stress tensor f, and the quantity 


P.,, the equation of motion of the electron gas becomes 


(2.7) lex (Ue* v)u, ) + NNieu(te—u)) =— vp.—ev (E+! ue. A H) : 


be may be considered as the pressure of the electron gas. It is to be noted 
that this equation is the same as the macroscopic one except for the addition 
of the term due to the Lorentz force and that due to collisions. The equation 
of the ion gas is similarly given by 


(2.8) o( 22:4 (u,-v)u,) +0 Ne. —t,) = —vertew (E+ bu, A H) 


where 


(2.9) 
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Formulation of equations of motion of the electron gas and the ion gas into 
eqs. (2.7) and (2.8) tells us that the behavior of a fully ionized gas can be 
pictured as that of two kinds of fluids, electronic and ionic, with a mutual 
interaction. This ‘two fluid’”’ model is the standpoint on which Schliiter (1952) 
developed his theory. Irrespective of its macroscopic appearance, it has micro- 
scopic bases except for the assumptions about the stress tensor and the force 
due to collisions. 

In order to consider waves of small amplitude we linearize both the above 
equations, and then we have 









‘ 0 e T- 1 
(2.10) Peo at N5ees (u.—U;) = —vp-eN,( E+! ue. A Hi) : 





(2.11) Pio oy Nj1e(U— Uu,) = verteN,( E+! UuU;A Hh) . 






peo and pio are the mean densities of the electron gas and the ion gas respec- 
tively. 

Adding the above two equations, we have the equation of motion of the 
whole gas: 









ou 
(2.12) po _ = —Vppt IAM, 






The meanings of the 









where it may be noticed that vg = vie and ay = a 


notations are as follows: 
(2.13) 


(2.14) Uy = (pecUe+pi0U;)/po, the mean velocity of the mass center of a 
volume element, 


Po = Peo tpi & po, the total density of the ionized gas, 
















(2.15) Pp = Pet+pi, the total pressure of the ionized gas, 





(2.16) j = e.,(u,—u,), the density of electric currents, 






Ho, the intensity of external magnetic field. 








As a result of the above definitions, we have 






1 Pi0 
217 u, = uo- 
( ) e _ eN, po 












9 = Pees 
- ™ wi & pe 


Multiplying eq. (2.10) by —pi) and eq. (2.11) by peo and then adding the 
two equations, we have the following equation after some manipulation: 





4nr dj, j 1 pio—peo ; 
( —s = le Sate 
(2.19) Sats = E++ u, A Ho aa j A Ho 














tage (Pi VPe —PpeoVp ) 
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where w, is the electron plasma frequency; 
(2.20) 
and 


5 Vei- 


a 


Equation (2.19) is the Pa ohmic law to be employed instead of (1.1). 
It is written as follows to a good approximation: 


0 1 
ante a ) = sia uy A Ho— al A Hoty vg (PV Pe peoVP1). 


Wp i 
This equation, lacking the last term, is used by Schliiter ee in his study 
of waves propagating along an external magnetic field. 
In the steady state, eq. (2.22) reduces to 


j = cEj+o,E) toad A E:) ; 


where E’ is the electric field observed from the system moving with the 
ionized gas: 


(2.24) E’ = E++ — 4 


and oy are Pedersen and Hall conductivities defined by 
25) C4. 

2.26) 

We can easily prove that 

(2.27) a 


where a, given by eq. (2.21) may be called the total conductivity. Equation 
(1.1) is the first approximation. Equation (2.23) is the next approximation, 
extending the ohmic law (1.1) to an anisotropic case, and it has often been 
used as one of the fundamental equations. However, it is inadequate for the 
treatment of high-frequency waves since it imposes some drastic restrictions 
(as seen from the above derivations), which we shall avoid in the following 
treatment. 

The other fundamental equations are the equations of continuity, which 
can be derived rigorously from the Boltzmann equation, 


ye+-div (pee) = 0, 
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(2.29) $614 div mad at 


In linear approximations they can be written as follows: 


(2.30) 20+ pug dite 


(2.31) 7+ py div u, = 0. 


3. DERIVATION OF DISPERSION RELATION 
We shall consider plane waves, in which wave quantities of the first degree 
may be assumed to depend on time and space through only one factor, 
(3.1) ei@ttker), 
In that case, differential operators for wave quantities of the first degree can 
be replaced by algebraic operators as follows: 


0/dt > iw 

Ve — ike 
9 ie 
(3.2) div A — i(k +A) 
rot A — i(kA A) 


The state of the electron gas is assumed to change adiabatically, and thus, 


| 


we have 
(3.3) Vie = CeVpe. 


C, is the sound velocity of the electron gas and is given by 


* peo a 


‘ +2 0 kT 
(3.4) ee 
where ¥-, x, and 7, are the ratio of specific heat capacities, the Boltzmann 
constant, and the temperature of the electron gas. The similar assumption is 
made for the ion gas, and thus we have 


(3.5) Vp1 = CiVpr. 
C, is the sound velocity of the ion gas; 


(3.6) Cy 
my, 


We assume that the temperatures of both gases are equal to each other and 
that y. = y;. Then, we have 


Ce ™, _ pio 
C 


m e Ped 


(3.7) 


Making use of the relations (3.3), (3.5), (2.30), (2.31), (2.17), (2.18), (3.7), 


and (3.2), we have 
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(3.8) Vp = VPet+ Vb 
* poC i J, r ’ a “s \ 
=e 2k saa (k j)ky . 


Similarly we have 


1 1 < 7: pol pol : 

( Seine File neg CeO Ailes 
(3.9) —. (pioVPe—peoVPs) = ar is (k- u,)k+é 7 a. : (k+j)k. 
Therefore, the fundamental equations (2.12) and (2.22) become 
(3.10) iwpou, = *j A Hot+i Cige (k-u,)k— Cine a (k -j)k, 

Np 
: ‘ 1 Ho 
bis) iia Hor Q 
Wp c Hy 
ae es ees 4 ss es pol 8 
a (k u)k++7 No) (k-j)k. 


In the following, we shall consider plane waves propagated in the x direction. 
The gas is assumed so rarefied that collisions are negligible. The external 
magnetic field is taken along the z-axis: 


(3.12) H,: H, = 0, H, = 0, H, = Hy = const. 







Equations (3.10) and (3.11) tell us that waves of the mode in which the 
electric field and the currents are parallel to the external magnetic field 
(z-axis) can be found without any motion of the ionized gas (uw), = 0). Waves 
in this case are electromagnetic waves, for which the dispersion relation is 








given by 







(3.13) 






where V, the phase velocity, is, of course, given by 







(3.14) V = fk. 
The above dispersion relation is familiar in propagation of electromagnetic 
waves in the ionosphere. 

Taking the x and y components of eq. (3.10), we have 


as S odlllig) ie = an 3,) 
(3.15) Ups = eN, 1—2(C3/ V’) @ it+i= ’ 


« . . ] Wy. 
8) mm” “oN, oo” 








where w, is the ion gyro-frequency defined by 







eH 


mic” 





(3.17) a, = 
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The electron gyro-frequency is given by, of course, 


eH 0 


(3.18) “= = 


and thus, we have 


(3.19) @e  ™ 


®, Me 


On the other hand, eq. (3.11) gives us 


(3.20) it “i . +e two Wejy—i oe poli ,2 


Upe tt 7 


eN, (eN : a 


(3.21) it ok os oe Sd 


D 
From Maxwell's fundamental equations, we can seis the electric field E, and 
the perturbed magnetic field, h, in this mode as follows: 


E. oo 
E, = 

Ei, = 

h, 


9 
(3.23) 4 


Substituting eqs. (3.23), (3.24), (3.16), and (3.17) into eqs. (3.21) and 
(3.22), and then making some rearrangements, we have 


nig W jWe w "Ce 1 1 =(Ci/1 V 2-9) - Weed i (Ci/ ra zz 
(3.24) (14! a tat Py TAC a Jz +1 w 1-2(C, 


3.95) wee 1 (Cy iin (%- teense ONIN eae ), bt 
C59) “a 1-CUV)* “Na 1-7) Ia 
D 


The determinant formed from the coefficients of the two equations must be 
equal to zero, in order that j, and j, may not be zero identically: 


eae w \"( a, \’ (C;/V’) ) (3 wu, wo Ce 1-—(Ci =) 
oa 3) @ (14,64 V") we ' w ww V1—-2(C,/V>) 
(*,- ee " __ WWe 1 ) at 
w, 1—(c°/V") w 1- 2(Ci/V*) é 


which is the dispersion relation. 


vx 


Put 
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wm a(S), nS nF rod. 
Then we have 
(3.29) Se oe == 
and 
(3.30) Y = ph. 
Using these abbreviations, the dispersion relation is rewritten as follows: 
(3.31) F(x) = Ao(a)x*—A,(c)x?+A2(a)x—AA3(c) = 0, 
where 
Ag(c) = uo’, | 


Ax(c) = o[f (u+2)+Apn}o— {(2+7+2p) +ru}], 

(3.32)  A2(o) (LEA (u+2)}o%=((1+y+29)+Aut4+y-+2—)}0 | 
+(1+p)(o+2a), | 

A3(¢) = o?—(2+7+2p)o+ (1+ )?. ) 


4. WAVES IN THE ABSENCE OF EXTERNAL MAGNETIC FIELD 
In this case, the dispersion relation (3.31) is reduced to 
(4.1) [ox —A(o—1)][uox?— { (ut+2)o—2}x+(c—1)] = 0. 


The first mode of waves, for which the dispersion relation is given by 


(4.2) gee a(1 4 
o 


are electromagnetic (EM) waves. One can easily find that eq. (4.2) is identical 
with eq. (3.9). Two other modes of waves are given by 


/ 9 
ph ae oe oe tf. €V 8 
(4.3) x = 14H, Ne Vb(1 , +7. 


The dispersion relations are illustrated in Fig. 1. 

Waves of (4.3) with + sign have the following properties: that x — ~1, 
as¢—© (w >w,), that is, waves tend to sound waves in the ion gas in a 
high-frequency limit, and that x > 4, as 70 (w Kw»), that is, waves in 
this case are sound waves in the whole ionized gas including electrons and 


ions. Waves, where x = 3}, or 


(4.4) V = /2C; 

may be considered as sound waves in the whole ionized gas. In a low-frequency 
limit, the differential velocity between the ions and the electrons may be small, 
and thus, the pressures of both gases, which are equal to each other under 
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Fic. 1. The dispersion curves of waves propagated in a fully ionized hydrogen gas. Ordi- 
nates: propagation velocity in cm/sec. Abscissae: wave frequencies in units of the plasma 
frequency. No external magnetic field applied. Cj and C, are the propagation velocities of the 
sound waves in the ion and electron gases. C; = 3.105 cm/sec is taken as an example. The 
electron temperature is assumed to be equal to the ion temperature. c is the light velocity, 
equal to 3.10" cm/sec. uw is the ratio of the ion mass to the electron mass. It may be noticed 
that the propagation velocity of the IS wave changes abruptly at a frequency w = wp/Vu. 


our assumption, would contribute co-operatively to wave propagation. The 
_mass density of the whole gas is almost equal to that of the ion gas, and thus, 
the square of the sound velocity of the whole gas becomes twice that of the 
ion gas. It may be noticed that the propagation velocity falls off rapidly from 
V2C, to C, around ¢ ~ 1/p, that is, w ~ w/Vu. At high frequencies, the 
wave velocity tends to that of the ion sound velocity, as stated already. This 
may be due to the fact that the ions easily move independently of the electrons 
at high frequencies. One can name this mode of waves ion sound (IS) waves. 
Waves of eq. (4.3) with — sign have following properties: x — 1/u, as 
o¢—o, that is, V—> ~V/uC; = C,. Waves in this case may be called electron 
sound (ES) waves. On the other hand, x becomes negative in the region 
a <1 (w < w,). Viewed from the definition of x, there is no wave propagation 
in this frequency range. In the frequency range o > 1, the dispersion relation 
for ES waves is approximated well by 


(4.5) w? = wrtCek? 
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which is a well-known form. If the gas temperature is zero (C, = 0), we have 
non-propagating oscillations, with frequency equal to the electron plasma 
frequency. Such oscillations are called plasma oscillations. Even though the 
temperature is not zero, a similar situation appears when the wavelength is 
so long that the frequency becomes independent of the wavelength. 


5. WAVES IN THE PRESENCE.OF AN EXTERNAL MAGNETIC FIELD 
In a high-frequency limit, the dispersion relation (3.32) gives us 





(5.1) x=, 
z (ut2)t Yes |i, 
9 = ov 
(5.2) x om = 11 
a 


The first corresponds to electromagnetic waves propagated with the light 
velocity. The second means sound waves in the ion gas and the third also 
sound waves but in the electron gas. The situation is very similar to that in 
the former section. 

In a low-frequency limit, the dispersion relation (3.32) gives us 


as p+2r 

5.3 oe ae ee 

(5.3) f ibs 

If p <1, that is, Vi «c?, the above solution approximately equals 
(5.4) h/x & p+2h, 


which is equivalent to (1.3). Note that the appearance of such hydromagnetic 
waves has been proved with the use of the condition: Vi « c®. If p> 1, that 
is, V; > c?, we have x & 1; waves are propagated with the light velocity. If 
the ion sound velocity is much greater than the Alfvén velocity, the propa- 
gation velocity is equal not to the Alfvén velocity but to +/2 times the ion 
sound velocity. 

In order to illustrate the dispersion relation over the whole frequency range, 
we take an example as follows; 


V, = 3X10’ cm/sec, C, = 3X10 cm/sec. 


Accordingly, we have p = (V,/c)? = 10-*, y = pu = 1.836X10-', and 
A = (C,/c)? = 10-'°. The results are shown in Fig. 2. As shown by the dis- 
persion relation, the three modes of waves do not appear independently but 
are coupled with each other through the external magnetic field. However, 
EM and ES waves are scarcely affected by the magnetic field. As in the case 
of the absence of the external magnetic field, they cannot be propagated in 
a frequency range lower than the plasma frequency, or at least, cannot be 
propagated without attenuation. The phase velocities of EM and ES waves 
tend to infinity around go ~ 1+~7/y and go ~ 1—vy. IS waves, on the other 
hand, are greatly modulated in a low-frequency range, where o < p, i.e., 
w? K ww; the phase velocity tends to the Alfvén velocity there. 
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Fic. 2. The dispersion curves of the waves propagated in a fully ionized hydrogen gas 
and normal to an external magnetic field are shown. Ordinates: propagation velocity in cm/sec. 
Abscissae: wave frequencies in units of the plasma frequency. C; and C, are the propagation 
velocities of the sound waves in the ion and electron gases in the absence of the external 
magnetic field. Cj = 3.105 cm/sec is taken as an example. The electron temperature is assumed 
to be equal to the ion temperature. Va, the propagation velocity of the Alfvén waves, is 
assumed to be 3.107 cm/sec. c is the light velocity, equal to 3.10" cm/sec. wi and we are the 
ion and electron gyro-frequencies. u is the ratio of the ion mass to the electron mass. It may 
be noticed that the propagation velocity of the IS waves changes abruptly at a frequency 
w = Vwiwe and at another frequency w = we. 


6. ATTENUATION OF HYDROMAGNETIC WAVES DUE TO COLLISIONS 

Let us consider the effect of collisions between electrons and ions. Multi- 
plying eqs. (2.10) and (2.11) by u, and wu, respectively, and then adding 
them, we have 


2) : (Apeote+ bps) + (Uy: Vie) + (ur: Vp1) si Npter(U1—Ue)”— (j-E). 


The first term of the right-hand side of eq. (6.1) is written as j?/o with the 
definitions of the density of electric currents (eq. (2.16)) and of the electric 
conductivity o (eq. (2.21)). The second term of the right-hand side (j +E) is 
rewritten as follows, making use of the Maxwell's fundamental equations; 


a «deal 
(6.2) (jeE) = div S+5, (e2) ‘ 
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where S, the Poynting’s vector, is given by 
. C 
(6.3) S =—E Ah. 
4n 


One can write 


2 
(6.4) (ue: Vee) = 5 2 C es) 


Peo 
making use of eqs. (3.3) and (2.31). Finally, we have 


> 0 1 2 1 2 + 2 Pe 2 pL E’ h* ° 7 
(6.5) ~- (Apeottet+ ps0) + aCe ts Cie + ae = divS+-. 


This is an extended Poynting’s theorem. The three terms in the curly bracket 
represent the kinetic energy of the gas, the potential energy of the gas, and 
the energy of the electromagnetic field. The first term of the right-hand side 
represents the source of the electromagnetic energy which is available for 
propagation and is given by 


(6.6) divS = i5-k-E Ah) 


_ Jae Vv’ lc? .2 
=s G-erre”’ 
This relation is derived from Maxwell’s fundamental equations. The second 
term, the Joule’s heat loss, represents the dissipation of energy due to collisions 
between electrons and ions. It may be noticed that 1/¢ is proportional to the 
collision frequency, vq (eq. (2.21)). 
In hydromagnetic waves, C; K V = Va Kc, we have 


2 
A? 
2 


Ju 


(6.7) ldiv S| = ae 
@ 


and this quantity should be greater than the Joule’s heat loss in order that 
hydromagnetic waves may be propagated without severe damping. Equation 
(3.26) tells us that 


. <a Da 
(6.8) lj2| = - liy| 
and thus, |j,| is much smaller than |j,|, since hydromagnetic waves can exist 
only in the frequency range, w K ~V/wywe < we. Therefore, j?/c is almost equal 
to j,?/o, and the condition of slight attenuation is given by 


(6.9) a Va »+ = 


® c 


Wey K Fw We. 
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It may be noticed that the condition is not governed by Pedersen or Hall 
conductivities, but by the total conductivity. In the inverse case of 


We >> 4W We 
hydromagnetic waves are heavily attenuated. 


7. CONCLUSIONS 


Three modes of plane waves, in which the perturbed magnetic field is parallel 
to the external magnetic field, can be propagated in a collision-free and fully 
ionized gas and transverse to the external magnetic field. They are coupled 
to each other through the magnetic field. In a high-frequency limit, the 
couplings are loosened, and the three modes of waves behave as electro- 
magnetic (EM) waves, electron sound (ES) waves, and ion sound (IS) waves. 
In the condition that the Alfvén velocity is sufficiently smaller than the light 
velocity and sufficiently greater than the ion sound velocity, the first two 
modes of waves are scarcely influenced by the magnetic field. Their phase 
velocities tend to infinity around the plasma frequency, and they cannot be 
propagated at all, or, at best, they propagate with attenuation in frequencies 
lower than the plasma frequency. The last mode of waves, however, can be 
propagated without attenuation over the whole frequency range, and are 
greatly modulated by the magnetic field at frequencies sufficiently lower than 
a critical frequency, which is given by the geometrical mean of the ion and 
electron gyro-frequencies. In that frequency range, they become hydromagnetic 
waves; they are propagated not with the ion sound velocity but with the 
Alfvén velocity. The hydromagnetic waves are heavily attenuated if the colli- 
sion frequency between ions and electrons is so large that the product of the 
collision frequency and the wave frequency is much greater than the product 
of the ion and electron gyro-frequencies. 

Waves in which the electric field and the currents are parallel to the external 
magnetic field are neither influenced by the magnetic field nor coupled with 
any other wave of the above modes. They can be propagated as electromagnetic 
waves at frequencies greater than the plasma frequency. 
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THE K-LL AUGER SPECTRUM OF Dy"! 
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ABSTRACT 


The K-LL Auger spectrum of ¢sDy'* has been studied using the Chalk River 
iron-free B-ray spectrometer at an instrumental resolution of ~0.07%. After cor- 
recting for the estimated intensity of the L; 43.81-kev conversion line (96% 
M1+4% E2), which coincides with the second (K-L,L2) Auger line, the following 
energies and intensities are obtained for the seven Auger lines: 

Energy 35.492 35.963 36.369 36.732 36.806 37.195 38.002 

(kev) +0.006 +0.007 +0.007 +0.007 +0.010 +0.006 +0.009 

Relative 1.00 1.31 0.22 0.83 0.28 2.39 1.07 

intensity +0.03 +0.05 +0.05 +0.05 +0.03 +0.05 =+0.08 
The energies are from 20 to 50 ev higher than the non-relativistic predictions. 
The relative intensities are consistent with experimental results at other atomic 
numbers. 













INTRODUCTION 


In previous papers (Ewan et al. 1959, 1960; Ewan and Merritt 1960) measure- 
ments of the K-LL Auger spectra in this laboratory at Z = 62, 64, and 94 
have been reported. The results confirm, qualitatively, the theoretical pre- 
dictions of Asaad and Burhop (1958), but are not in quantitative agreement 
either as to energies or relative intensities. However, the experiments have 
confirmed one important feature predicted by their theory, namely, the 
existence of the 1,L;(°P;) satellite line. This supports their assumption that 
the coupling between electron-shell vacancies is neither pure j—j nor pure 
L-—S since more than six lines are observed. On the assumption that the coupling 
is intermediate between L—S and j-j, they predict that the K—LL Auger 
spectrum should have nine lines, i.e., six main lines plus three satellite lines 
of low intensity. The predicted intensities and energies of the latter are such 
that only the L,L;(*P;) satellite should be easily detected experimentally. 
The other two are predicted to be less intense and closer to their partners, 
hence much more difficult to observe. 

In this paper we report measurements of the K-LL Auger lines in e¢gDy™. 
These results were obtained during a detailed study of the internal conversion 
spectrum following ¢;Tb'*' beta decay (Graham et al. 1961), which will be 
described fully elsewhere. The results again confirm the existence of a L;,L3(*P1) 
satellite line observed in earlier experiments. The relative intensities at 
Z = 66 are consistent with those at neighboring values of Z. The systematic 
variation with Z, indicated by other measurements, is in reasonable agree- 
ment with the theoretical relativistic intensity ratios at Z = 80 (Asaad 195)). 




































EXPERIMENTAL PROCEDURE AND RESULTS 
Tb’ was produced according to the sequence Gd!(n,y)Gd!"8- — Tb!*, 
Enriched Gd! was irradiated in the NRX reactor at a flux of 6 X10" n/sqcm/ 


1Manuscript received April 13, 1961. 

Contribution from General Physics Branch, Atomic Energy of Canada Limited, Chalk 
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sec for 30 days. Following irradiation the two elements were separated by 
elution from an ion exchange column with ammonium a-hydroxy-isobutyrate 
(Smith and Hoffman 1956). To purify the product from the eluting agent 
the Tb'® fraction was acidified and absorbed on a very small Dowex-50 
column (0.03 cm?X5 cm); the bed was washed with 30 column volumes of 
0.2 N HCl, and the Tb'* was eluted with 6 N HCl. Polythene ware was 
used throughout the procedure to minimize contamination by material from 
container walls. The final product was sublimed from a boat-shaped tantalum 
filament through a slit 0.030 in. wide onto 800 ug/cm? thick aluminum foil. 
Several different sources were prepared in this way; the results reported here 
were obtained from one having an initial strength of ~ 100 uc and an estimated 
surface density of <50 ug/cm*. 

The Auger spectrum was studied in the Chalk River double-focusing 8-ray 
spectrometer (Graham et al. 1960) at a resolution setting of ~ 0.07%. The 
focused electrons were detected in a side-window continuous-flow proportional 
counter having a thin (100 wg/cm*?) VYNS plastic film window. The 
window was supported by a fine mesh metallic grid (~90% open area) and 
the methane filling was maintained at a pressure of ~8 cm Hg absolute. The 
detector efficiency under these conditions is essentially constant, and equal to 
the open area of the window grid, for the electron energy region of interest 
here. 

The K-LL Auger spectrum is shown in Fig. 1, plotted on a linear scale. The 
data were taken automatically and recorded both by an electric typewriter 
and on punched tape. Six successive sweeps were made of this region of the 
Tb!" spectrum using 40-second counting intervals. The six runs were added 


\ 


es ae egy ae 

'e 'p, 's °P, *p 'D ) *p) 

L 4581 —~ cP (So) (Pull Pe) ( 2 ( 2 
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ELECTRON MOMENTUM IN GAUSS CM 


Fic. 1. The K-LL Auger spectrum in gsDy" following 8-decay from ¢;Tb'. The shape of 
the background from the 8-continuum and tails of higher-energy intense conversion lines is 
shown by the thin solid line. The components of the complex lines are shown by broken lines 
soe Seca in the text. The L,L;(*P2) satellite line is predicted theoretically by Asaad and 

urhop. 
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and corrected for variations in counter background rate (air-borne A“! activity 
discussed by Graham ef al. 1960) directly from the punched tape record 
using IBM accounting equipment. The Auger peaks are all wider than the 
instrumental peak shape for two reasons, i.e., natural width and source thick- 
ness. The former effect accounts for the difference in width between the L 
conversion lines of y 43.81 and the Auger lines; the latter is indicated by the 
pronounced asymmetric low-energy tail on all lines. 

The measured energies and relative intensities of the Auger lines are sum- 
marized in Table I. The energy uncertainties listed are intended to represent 


TABLE I 


Comparison of experimental and theoretical energies and relative intensities of the K-LL 
Auger lines in g¢Dy!™ 


Relative intensities 
Energy in kev 











— Non-rel. Experimental 
Non-rel. theor.6 -—— — 

Transition Exptl. theor.? Z = 66 sey! 625m82° 
L,L\('So) 35.492+0.006 35.455 1.00 1.00+0.03 1.00+0.05 
L,L('P)) \ 35.9630 .007 { 35.938 1.27 1.310.054 1.39+0.06 
L,L(3Po) 35.972 0.047 

L2L2(4So) 36.369+0.007 36.350 0.17 0.22+0.05 0.26+0.06 
L,L3(8P}) 36.732+0.007 36.698 1.90 0.83+0.05¢ 1.18+0.06 
Lj L3(3P2) 36. 806+0.010 36.763 0.24 0.28+0.03¢ 0.29+0.09 
LoL3('De) 37.195+0.006 37.148 4.43 2.39+0.05 3.14+0.13 
L3L3(*Po) \ { 37.904 0.25 | 

L;L,(3P2) { 38.002+0.009 \ 37.950 2.03 f 1.07+0.08 1.39+0.06 





“Using the binding energies of Bergvall and Hagstrom (1960). 
’Asaad and Burhop (1958). 

“Ewan et al. (1960). 

4After subtraction of the L3 conversion line of the 43.81-kev y-transition as shown in Fig. 1. 
*The total intensity of this partially resolved doublet is 1.11 +0.04. 


standard deviations. They are a combination in quadrature of the possible 
errors in absolute calibration and those due to statistical fluctuations. The 
spectrometer was calibrated using sources of Bi? (ThB), Au’, and Cs'*, 
Corrections were made for a small asymmetry in source position relative to 
the dowel mounting pins. This was determined by rotating the source mount 
ring through 180° and remeasuring the positions of several intense conversion 
lines with the source in the inverted position. Throughout these runs the 
degaussing system currents were adjusted automatically to compensate for 
fluctuations in the earth’s magnetic field, using the electrical signals from an 
improved version of the Serson three-component magnetometer (Serson 1957). 
This control system now maintains the mean value of the vertical field compo- 
nent at the spectrometer to <+1X10~ gauss even during severe magnetic 
storms. Errors in line positions due to uncertainty in degaussing are now 
thought to be negligible, i.e, <+0.01 gauss cm or ~+1 ev for the Auger 
lines in question. Errors due to lack of reproducibility in source positioning 
are also negligible since the source rings are machined accurately to fit snugly 
on a pair of locating dowels; the position is reproducible within ~ 0.003 cm 
(Graham et al, 1960), which corresponds to an error in Auger-line energy of 
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<0.6 ev. Using the K line of the 411.770-kev transition in Hg! as our basic 
standard (momentum = 2222.39+0.14 gauss cm) the more intense conversion 
lines in Dy’ show an internal consistency of <2 ev if we employ the binding 
energies of Bergvall and Hagstrom (1960) including the work function, zp. 
The consistency expected for the Auger lines is somewhat poorer for two 
reasons: (a) the statistical accuracy is poorer, and (b) there were small drifts 
in current settings during the long automatic run. The current settings were 
measured before and after the six sweeps and showed a drift of ~2 parts in 
10°. The positions of the Z; and LZ» 43.81-kev lines were measured during 
subsequent calibration runs and hence provide a secondary calibration for 
the K-LL Auger lines. Their energies agree with the later measurements to 
within +4 ev. The relative energies are more accurate than the errors indicated 
in Table I by perhaps a factor of 2. 

The intensities of the various peaks have been deduced from the line areas. 
In order to calculate areas it is necessary to know the shape of the spectrum 
under each line including the magnitude and shape of the tails of the lines at 
higher energy. For purposes of analysis it was assumed that all Auger lines 
had a shape the same as that of the prominent K—Z2Z; line. The height and 
slope of the background line shown in Fig. 1 is that found empirically which 
satisfies this criterion. The increasing slope at high momentum is due to the 
tail of the very intense L, 48.91-kev conversion line at ~686 gauss cm. The 
matching of line shapes was carried out graphically using semilogarithmic 
plots as described in earlier papers (Ewan et al. 1959; Graham et al. 1960). 
The intensity errors listed in Table I are intended to represent standard 
deviations although the method of analysis is somewhat subjective, i.e., the 
results depend upon the shape of the tails and the choice of the background 
shape which is subtracted. The errors quoted are those from counting statistics 
combined in quadrature with an arbitrary 10% error assigned to all estimated 
tail areas and an error ranging from +30 counts to +100 counts per 240 
seconds for the background line. 

The energy and intensity one deduces for the K-L,L, Auger line depend 
upon the contribution assumed for the unresolved L; 43.81-kev component. 
The 43.81-kev transition has been assigned as 96% M1+4% £2 and is 
identified as a transition de-exciting the first rotational state based on the 
5/2+ ground state of Dy'™ (Graham et al. 1961). On this basis the intensity 
of the L; line should be equal to that of the LZ, 43.81-kev line as shown by 
the broken line in Fig. 1. The K—Z,Z, line intensity error quoted in Table I 
includes allowance for a +20% uncertainty in the intensity of the ZL; 43.81-kev 
line; the error quoted for its energy includes an allowance for uncertainties 
in both the position and the intensity of the L; 43.81-kev line. 


DISCUSSION OF RESULTS 
(a) Line Energies 
In column 3 of Table I are listed the energies predicted theoretically for 
Z = 66 by the formulae of Asaad and Burhop after insertion of the binding 
energies of Bergvall and Hagstrom (1960). The measured energies are from 
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19 to 49 ev higher, i.e. by amounts too large to be ascribed to experimental 
errors. The discrepancies are comparable to those found previously for Z = 62, 
64, and 94 (Ewan et al. 1959, 1960; Ewan and Merritt 1960) and are here 
mainly attributable to the inadequacy of the non-relativistic energy formulae 
of Asaad and Burhop (1958), which includes semiempirical correction factors 
based on the measured K—LL energies at Z = 83 (Mladjenovic and Slatis 
1954). Relativistic energy predictions are not available for Z = 66. In his 
recent paper, Asaad (1959) does, however, make relativistic energy pre- 
dictions at Z = 80. The discrepancies between these and the measured 
energies of Bergstrom and Hill (1954) at Z = 80 are as large as 200 ev, which 
suggests that the present theoretical techniques are not capable of predicting 
Auger-line energies with the precision now obtainable experimentally. 


(b) Relative Intensities 

As was discussed in an earlier paper (Ewan et al. 1960) there is always the 
possibility that the apparent Auger-line intensities may be altered by low- 
intensity conversion lines close to or coincident with the lines in question. 
The L,L2('P:)+JL; 43.81 unresolved doublet is a good example. The presence 
of the 43.81-kev transition in Tb'™ decay became certain only after a detailed 
study of the entire conversion electron spectrum had been made and the 
Dy'®* level scheme was firmly established. Even after substracting the pre- 
dicted L; 43.81-kev contribution, the remaining intensity assigned to the 
K-L,Ll, Auger line is significantly higher than for Z = 62, as is also the 
K-L,L, line. The other Auger-line intensities are, however, in acceptable 
agreement with the Z = 62 values. 

In Fig. 2 we compare the relative intensities reported here with the theo- 
retical predictions and experimental values at other values of Z. The intense 
K-L.L;('D>») line has been selected as the reference rather than the less intense 
K-L,L,(4So) line. The solid lines are the non-relativistic predictions of Asaad 
and Burhop. The relativistic predictions at Z = 80 are indicated by the 
arrows. The broken lines show the trends of the experimental intensity ratios 
with Z and agree well with the relativistic predictions at Z = 80 except for 
the weak K—L,L.2(!S)) line where the prediction is too low. The new results 
are consistent with the trend indicated by the other experimental results 
within errors. This comparison again points up the inadequacy of the non- 
relativistic theory. The latter does, however, predict with reasonable accuracy 
intensity ratios of the three lines involving the Z; shell and also those for the 
three lines not involving the LZ; shell but, at high Z, is in error by a factor of 
3 to 4 for the ratio between these two groups. 

The present results at Z = 66 once again provide convincing evidence for 
the K-L,L;(*P2) satellite line and confirm the previously reported evidence 
for it at Z = 62, 64, and 94. The statistical accuracy of the present results 
permits a more accurate estimate of the intensity of this line than was possible 
in the earlier work at this laboratory although it is not fully resolved. Its 
intensity is consistent with the earlier results at Z = 62 and 64 but the experi- 
mental ratio for *P2 to *P; is about a factor of 2 larger than predicted in the 
non-relativistic theory. The corresponding ratio is not explicitly available 
from Asaad’s relativistic calculations at Z = 80. 
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Fic. 2. Comparison of K-LL Auger-line intensities relative to the intense K—L2L;('D2) 
line. The non-relativistic predictions of Asaad and Burhop (1958) are shown by the solid 
lines. The broken lines are guides to the eye through the experimental points. The relativistic 
predictions of Asaad (1959) at Z = 80 are shown by the arrows. The experimental points 
are from the following sources: Z = 47, Johnson and Foster (1953); Z = 55, Graham et al. 
(1961b); Z = 62, Ewan, Graham, and Grodzins (1960); Z = 64, Ewan and Merritt (1960); 
Z = 66, present work; Z = 78, Ewan (1957); Z = 79, Hill (1953); Z = 80, Bergstrom and 
Hill (1954); Z = 83, Slatis (1955); Z = 84, Gray (1956); Z = 94, Ewan ef al. (1959). The 
error bars for the Z = 55, 62, 64, 66, and 94 results indicate standard deviations; the results 
at other Z values are from photographic intensities and have been assigned errors similar to 
those estimated by Ewan (1957) for his photographic results. A value of (K-L2L2)/(K-LiLs) 
= 0.08 has been assumed for the purposes of plotting the Z = 55 results. 


(c) Absolute Intensities and Fluorescent Yield 

The sum of the intensities of the K—LL Auger lines relative to that for K 
conversion lines in Dy' is 0.0375+0.004. An accurate measurement of the 
intensity of the K-LM, K-MM, etc. lines has not been attempted in this 
work since these spectra are much more complex and difficult to interpret. 
Moreover, they are partially masked by conversion lines. However, the 
intensity ratio [((K-L Y)+(K-—XY)]/|K-LL] is a slowly varying function of 
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Z (see, for example, Gray 1956 and Wapstra et al. 1959). Taking 0.6+0.1 as 
the value for this ratio at Z = 66 one deduces that the fluorescent yield at 
Z = 66 is w, = 0.943+0.007. This is to be compared with the value of w, = 
0.926 for Z = 66 computed from the semiempirical formula quoted by Wapstra 
et al. (1959). 





CONCLUSIONS 

The present experimental results again provide convincing evidence for 
the most prominent satellite K—L,L;(?P:2) line predicted by the non-relativistic 
theory of Asaad and Burhop and confirm earlier evidence from this laboratory 
of its occurrence at Z = 62, 64, and 94. It is therefore clear that the final 
electronic configuration of the atom must be described in terms of a coupling 
between the two L-shell vacancies which is intermediate between j—j and L-—S 
even at the highest atomic numbers. The six-line spectrum predicted for 
extreme j—j (or L—S) coupling is not consistent with our experiments. 

The relative intensities of the K—LL Auger lines are shown to be consistent 
with five of the six intensities predicted by Asaad (1959) from his detailed 
relativistic calculations at Z = 80. The K—L:L2('So) lines observed experi- 
mentally are, however, appreciably stronger than this theoretical prediction. 
A fully relativistic calculation at lower Z values, taking intermediate coupling 
into account, would provide a more rigid test of the validity of Asaad’s assump- 
tions than is presently possible. 
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A NOTE ON THE DIFFRACTION OF A SCALAR WAVE 
BY A SMALL CIRCULAR APERTURE! 







R. A. Hurp 







ABSTRACT 


The power series solution for the diffraction of a scalar plane wave incident 
normally on a small circular aperture in a hard screen is extended to include 
terms of order a! in the aperture field, and to terms of order a” in the trans- 
mission coefficient (a = 2 X radius of aperture/wavelength). To do this, a set 
of recurrence relations is developed, and a relatively simple method of obtaining 
the transmission coefficient is devised. 









INTRODUCTION 






In the course of some recent work on diffraction theory, a power series 
solution was obtained for the diffraction problem of a plane wave normally 
incident upon a small circular aperture in a hard screen. It is believed that 
the power series (in a = 2ra/\, where a is the radius of the aperture) has 
been carried somewhat further than formerly. Thus the series for the aperture 
distribution function includes terms of O(a?) while the transmission co- 
efficient 7 contains terms of O(a!°). In addition, recurrence formulae have 
been developed which allow (in principle) the calculation of the aperture 
field to any order in a. Also, a relatively simple expression for the transmission 
coefficient of the aperture has been found. 












FORMULAE AND RESULTS 


An integral equation derived by Bazer and Brown (1959) forms the starting 
point of the analysis: 






h(y) dy, 


(1) h(x) = cosh ote f — y) 





where 8 = —i/x. The (suppressed) time dependence is e~'. The function 
h(x) in equation (1) is essentially an aperture distribution function and is 
related to Bazer and Brown’s function f2(x) by A(x) = (r/a)fo(x). 

A solution to equation (1) was found by assuming an expansion of h(x) of 
the form 


(2) h(x) = : ah, (x), 


n= 
















{n/2] 


(3) h(x) = D0 S30. 


m= 






The notation [pf] in equation (3) means the greatest integer <p. 
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TABLE I 


The coefficients ©, 
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: 8 8 ;, 64°55 7 7? 22 88 .;, 16 
fo = 3905 8 +9 g+5 B+1288, $2 = Tes B+ 735 B+ 3 8 


7 1 1 os 7 
4 = 519 8 +75 8 t= 






1 
3520 6 
112, 1856 a 746 


1 8 hios 8 
=~ >> B+ B+ f 6 £8 as 48 432 
33,075 3 675 ° 4 P+2568", 3 14,175 8 6° 3 ea 





















5 2 2 4 pe Se 2 —_ 1 
= g7gF t758) s8 = T5698» Ss = 79,390" 


4 ar’. dale: Gekko a 512 .,. 64 > 
g4 ee gts g°-4——— B7-+5128", a z 7 8’, 
515 °* 35,721 ° * 2905 8 tg 8 +9126 = 79,815 a+; 567 ae +735 8 +3 8 


“~ 


9 
[4 = 


= 23 RD ike 4 i ae se 1 ; 1 
[3 —— nae & +: - ’ f6=s- ron C3 = SO. fan 
aOR Sap) B+o55 8 +758 6 8505 b+a35 Bg, 3 181,440 8 


10 2852 2, 78,224 14,272 oy 
[se roe" a ( 02 
sg 826,875 , +997,575 B+ 675 3°+2568"+ 10243 


10 3844 9, 22448 4 1184 6 198 5 OT 
& = 396195 © +g0'505 8 +5g5 +3 8 | = Tah +g 66 ti5 


10 os Be! 10 ca ae 2 _t8 tee 1 
" sam "tai a, te 90,720°* = 34598 800° 


8 509,368 73,856 .5 , 36,992 ; ( 
rc ae sae 34 09,36, ; 73,856 -% 36 2 i+ _ a 4 20488", 
5 


5,175" | 22,325,625" “35,721” " 675 


1 74 152,224 .. 20,864 5, 896 7, 256.5 
fo = 8 B+ — 3°4——— | - p° 
$2 1,408,325 .% 165, 125 8 14,175 8 145 Og! 


ets 74 223 | , 184 16 | 37 
{= Too1a75 Ot aT +ou5 O +45 















os dBi iat ely. 70s alee inet ei 
sigs 623,700 b+75 215 8 +35 6° $8 598,752 a0 8°, $19 = 19,958,400 * 




















446,024 623,024 1,910,528 55,808 2 264 10 
: Bon gt tw t0028 96 11,264 5104 i963" 
* 7,620,840,375 ° 5,740,875 297,675 ® + 405 8 ss + 


12 
0 


“yy 




















6016 , 
135 | 


512 10 

a+ 7a" 8", 

2 6029 oo. eke gt+—t 
me 199, 500 8,505 315 


12 _ _ 214,682 238,328 24, 18,464 
a ed 581,875 9 +1,013,025 a+ 1,725 B+ 























21,433 848 124 


f= 79,116,375 8 +30, 375 9 +55 e+e 6" ? 


om 








A at gt a ae. Bite aera 
8,981,280" ' 22,680" ! $10" 9,979,200" * $12" 479,001,600 ° 
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TABLE II 


The coefficients p; 


aE 
Pn n° 
1 _ 28 
Px = =. 
2_ 1 2 y a 4 2 . 2 5s ee Se 
aS gt 48 ’ p= pts B’, Ps “ ore a. 7 = 18° 7 B - py = To a 

5 ‘ = 2B oS + « Se Si 
pi = 5 8+88", Pi = igs B+ 8’, i 

2 19 8,3 oo ee 

pt = 739 B+z 8, ps = 301 a+s 8°. 
ie ‘ne vies 26 16 5 4 
a = ote t 5?4168", p= rea tag’ +5 

i Re 16 44 ae See. apy oe 

_ neti? +5 B, = aE; t 199% +7 8 
5 _ 233 > 6 _ 239 , 196 33 = a 
Pt = 3799 BAB +328", wi = Fag BG g5 Oty 8 

s 2383 am s 11,857 3° 

Pe = 94500 °* 315 38 +7 5 8 P? = 593 700 8 a 8 +4 

61 6 ’ 2 72 

of = 4 SO +5 5 t+ 648", ph me hg TOE py STE 0 5 SO ot 


‘eas 1350 6480 | 85,050 135 3 


6 1 12,737 2 , 136 44, 64 


ot = 7930+ 141,750 8 + G3 +5 O 
vi = spr gon tapas tg A +1286", ok = ato ng tag Ot GS ets 
ph = Sao tag Ot aes a+ Ge at , r 6°. 
A= sap set ae 8 +223! 3 Bi4 ~ 6°+2568", 
pa = 4 eis | Saas B+ re B+ ae +5 +s 
A= seerenmatirenase ot ag Ot gO tSI20" 
* 760,913 100,453. 2 , 28,996 5, 5696 27, 512.» 
Ps * 1,886,068,800 ° "1,071,630" "8,505 135 ° 3 8: 
pL = saunas ties! +pegeag #4 + aos ot +2 ot +1024". 
1 _ 8,084,789, 1,107,979 4 , 464,026 5 , 4624 oe 6°+-20388" 


Pi 207 467,568,000 °° 30,618,000” * 178,605 75 


When equations (2) and (3) were substituted in the integral equation it 
was found that the £2, satisfied the following relation: 


((n—1)/ bs ) 1 is 
mes Si ——— aon ee 2 ! 
(4) F2m 2B au 2m (2r+1)!”? 2r— Sw titée ak) /n ’ 


in which 
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TABLE III 


The coefficients T2,_ 


6.2 
=e a * 2°+328"'. 


1136 , 3872 256 
= 33,0757: 2025 s +> B*+1288°. 


19,046 _ _, 36,336 3? 608 a 


a —— D / 
= 3,465,125 178,603 45 P tg O° t5128". 


TS 


658,136 1,280,408 161,408 56,576 6 aus 


fies 0 W,J9tU 5 10 
7 = 17§90,840,375t 13,395,375 ° + 35,721 8 + 675 +g A+ 20888". 


2) 
m=0 t+2m S2ms 

and where 62,,, = 1 if mn = 2m and be, = 0 otherwise. A derivation of these 
equations is given in Appendix A. 

Equations (4) and (5) constitute a set of recurrence relations from which 
the solutions to the integral equation (1) may be built up. Table I shows 
values of the coefficients (2, up to and including nm = 12. It will be observed 
that the results agree with the values given by Bazer and Brown (1959), who 
give the expansion of f(x) up to and including a®. In Table II the coefficients 
p’ are given. As is shown in Appendix B, the p{ are useful in the computation 
of the transmission coefficient 7, the relations being: 


(6) t= Dd arm, 


n=) 


(5) pi = 


(7) rm = —46 > - Son Seman 


m=( (2m 


The coefficients 72, were calculated up to and including 2” = 10, and the 
values are shown in Table III. These results agree with Bazer and Brown, 
who give terms including 7.6. 


TABLE IV 


Comparison of the exact and power series values of 7 


Exact Power series Exact Power series 


2 


883 . 88480 
. 900 . 90122 
917 .91777 
932 93413 
945 . 95066 
958 . 96895 

. 968 . 99297 
976 1.03061 


8106 . 81057 
.8118 .81184 
8156 . 81563 
$219 .8219] 
8306 $3058 
8415 $4150 
854 $5442 
. 868 86905 


>~wmOOWS 
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= 
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A comparison of the power series with the exact results is given in Table IV. 
The exact values were taken from Bouwkamp (1954, p. 84). It will be observed 
that 1% accuracy is attained up to a ~ 2.5. 
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APPENDIX A 
DERIVATION OF THE RECURRENCE RELATIONS 
If the power series representations of cosh ax, sinh a(x—y), and h(x) are 
employed in equation (1), one obtairis 


! - & fox)" fas Se atia-y)™ $ ’ 
(A.1) Y, a"ha() > (nyt? od 2 Gm-d)! 2, 2"htp(y). 


After rearrangement of the power series in (A.1), the powers of a may be 
equated throughout to give 


n vl [(m—1) /2] ) 2m 
5 £ (x—y)  ta—2m—1(Y) 
A9 Viaar aoe , 
(A.2) Iin(x) = Yn +8) 4 x (2m+1 yr ; 
where y, = 1 if 2 is even, y, = 0 if m is odd. The binomial expansion 
Pp 
(A.3) (x—y)? = (—)? DY (—x)? (?) ' 
q=0 ¢ 


along with equation (3), is introduced into equation (A.2): 
[n/2] _ Yk ((m—1)/2) 1 [(n—2m—1)/2) ; el 
7 2m ne s%— 2m— 
(Ad) SS ee = Me Do . = he 
=i 


m= m=( (2m + l )! p=0 
2m 2m ‘ 
2b Ege ey 
q 


ql 


The integration with respect to y is now carried out. Note that odd powers 
of y vanish upon iit The result is 


[n/2] , Vax [(m—1)/2] 1 {(n—2m—1) /2] = ‘ie 
(A5) >> tine = 2-428 Y= ze. te ae 


a mao (2m-+1)! $=6 eno Nog J 
on 2¢ 
2p+2q+1° 
The abbreviation 
t2l 


2n 


(A.6) 3 = 


sno (+2 ° 
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is put into equation (A.5), with the result 


sO 2 x" (ne? > 2m 2m—1._2m—2 
A.7 oo = 419 n—2m—1,.2m—2¢ 
~) 2 S tm ia 2, ari q20 \2g rants 
The substitution » = m—q reduces the above equation to 

{n/2] [(n-1)/2] - om : 
A.8 ‘ en = +2 qnitiainieasincainie ere x”? 
oe 2, sa ae 2, ae ! 2, nro 
If like powers of x are sean one finds 

((n—1)/21 /¢ 
x iis 2m 1 Sn 

A9 4, = Sense ( —s soma 
‘ini - 428 20 55) Gm-tiyi me 


This is the desired formula. 
APPENDIX B 
THE FORMULA FOR To, 


The transmission coefficient 7 is defined by 
(B.1) 7 = —(2/aa)Im[A (0)], 
where A(0Q) is the amplitude of the far field, 


1 
(B.2) A(0) = af cosh ax h(x) dx. 
a 
Thus, from equation (6) one has 
co ‘ 1 
(B.3) > ar, = —— cosh ax Im[h(x)] dx. 
n=0 TAS} 


Since fo(x) = 1 it is easy to see from equation (A.2) that h,(x) is imaginary 
only when 1 is odd. This being the case, one has 


oo ie 9 al 0 x 2m oo S 
(BA) Lama 2 ax d Sy Ea imfhayss(e)} 
n=0 —_ m=( 2m p=0 


When the sums in (B.4) are rearranged, and like powers of @ are equated, 
it is found that 


9 el n Pei 
(B.5) fe = 2 { x . Qn — 9p)! ; Im[heop+1(x)]. 
J _ p=0 (on 
But 
(B.6) Im[hopy1(x)] = —7 > petty “ae 
q=0 


Putting (B.6) into equation (B.5), rearranging the series, and integrating 
gives the result 


a di : oo 2p+1 
(B.7) Ta, = e e 3 (2n— ae 2p+1- 


With the substitution m = n—p, one now obtains the desired formula equa- 
tion (7). 
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ESTIMATION OF THE EFFECT OF SURFACE ON THE DIPOLAR BROADENING 
OF MAGNETIC RESONANCE LINES IN CRYSTALS* 


G. C. BENSON AND P. BALKt 


The broadening of the nuclear magnetic absorption lines of a crystalline 
powder by dipolar interaction should depend on the state of subdivision of the 
material since the interactions of nuclei in surface regions differ from those of 
nuclei far removed from the surface. The present note provides a simple 
method of estimating the order of magnitude of this surface effect. 

In the following consideration, it is assumed that the crystalline powder is 
an ionic solid having a NaCl structure in which the distance between nearest 
neighbors is equal to a. The sample contains N positive ions and the same 
number of negative ions; ” pairs of these ions lie in the crystal faces which are 
assumed to be of the {100} type. It is further assumed that the lattice structure 
is rigid and perfect throughout the individual particles of the powder and that 
although (n/N), the fraction of the ions located in surfaces, may be significant, 
the crystals are not small enough for other structural features (edges etc.) 
to be important. If the powder has a specific surface, 4(m? g~!), it can be 
shown that 


(1) n= 10'ap AN 


where p (g cm~*) is the density of the crystalline material. 

The second moment, S2, of the line associated with transitions between 
states of the positive ion with different nuclear spin is given by van Vleck’s 
formula (van Vleck 1948), which in the present case can be written 


(2) S, = N™ w3| § 1. Ie+Vet D rd+75 L(1_+1)g¢2 Eri]. 


In this equation yo is the value of the nuclear magneton; J and g denote, 
respectively, the nuclear spin number and gyromagnetic ratio of the ionic 
species indicated by the subscript + or —. The distance between two nuclei 
i and j is r,;. The first sum in equation (2) is taken over all possible pairs of 
positive ions; the second sum is over all pairs of unlike ions. These sums are 
closely related to other sums occurring in calculations of the cohesive energy 
and the surface energy of the crystal, thus 


(3) dry = a °3(NA¢—nBe] 


*Issued as N.R.C. No. 6329. 
t+National Research Council of Canada Postdoctorate Fellow 1956-58. Present address: 
International Business Machines Corporation, Research Center, Yorktown, New York. 
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and 





(4) DD ry = a “(NA—nBe] 





where Ag and Aj, have been evaluated by Jones and Ingham (1925) and 
Bg and Bs by van Zeggeren and Benson (1957). 

Combining equations (1) to (4), the second moment can be written in the 
form 


(5) Sz. = St+AS2 








where 

(6) S: = 438 1.(I4 +1) gf: Ae+e 4 tnetai| 

and 

(7) AS: = —10'a- ‘pla iE LLANE BY +s r+ieBi |. 





The negative sign in equation (7) indicates that the line is narrowed by 
increasing the specific surface of the sample. 

In the case of NaCl itself @ = 2.820 10-* cm and taking the other necessary 
data from a review article by Pake (1956) it follows that 


(8) S2(NaCl) = 0.542 — 1.17X10~ A gauss?. 









The value of S} in this equation is very close to the experimental results 
obtained by Otsuka, Oshio, Kobayashi, and Kawamura (1959) with a single 
crystal of NaCl oriented in two different ways in the external magnetic field. 
For a powder with a specific surface of 100 m? g—, the present theory predicts 
a decrease in the second moment of about 2%. This change is of the same 
order as the experimental error in current line width determinations. Hence 
it appears that such a surface effect would not be detectable until the specific 
surface of the sample exceeded 100 m? g~'; this means that individual particles 
would contain less than 10° molecules. 

Recent theoretical studies of the configuration in the surface region of an 
ionic crystal have shown that the equilibrium positions of the ions in this 
region differ slightly from the regular lattice sites. Thus a refinement of the 
present treatment is to correct AS». for this lattice distortion. Adopting the 
one-layer model examined by Benson, Balk, and White (1959) and using the 
same notation as these authors, the effect of surface distortion is to add a 
correction 
















(9) AS: = 10'0*sipa ° et DS! Gro 





+5 LAE AS! Ul U,[6] (<4 £-+22)+ 2? ual6le—24)"* | 
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to AS». After substituting numerical values for NaCl, taken from the preceding 
reference, equation (9) gives 


(10) AS3 = 0.44710 A gauss’. 


This correction tends to decrease the surface effect estimated previously. It 
may be concluded that if any surface effect is observed in cubic crystals it 
must be due to interactions of other than a dipolar nature. 

All the calculations in this note have been described with specific reference 
to the NaCl structure but other crystal types could be treated in a similar 
way. There is, however, no reason to expect that the dipolar surface effect 
would be appreciably different in magnitude for other structures. 
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ANGULAR VARIATION OF AXIAL SPIN-HAMILTONIAN EIGENVALUES IN 
STRONG MAGNETIC FIELD 


H. A. BUCKMASTER 


Tables of eigenvalues and matrix elements of transition probabilities are 
required in the design of those solid-state masers which are based on the 
three-energy-level scheme proposed by Bloembergen (1956). Chromium- 
doped ruby has been widely exploited as the active material in masers because 
of the simplicity of its energy level scheme, its chemical stability, its high 
thermal conductivity, and the ease with which it may be machined. It is not 
surprising that the above-mentioned tables have been calculated for chromium 
in ruby, i.e., for S = 3/2 in an axial crystalline electric field by Weber (1959), 
Schulz-Du Bois (1959), and Stahl-Brada and Low (1960). More extensive 
tables will be required when other satisfactory active materials for masers are 
found. Such tables are also useful whenever it is necessary to locate magnetic- 
resonance fine-structure transitions corresponding to AM = +1, +2, +3, etc. 
where M is the strong magnetic field quantum number. This is particularly 
true when the applied external magnetic field B is at an arbitrary orientation 
to the symmetry axis of the crystalline electric field. It is convenient to describe 
the orientation of B with respect to the symmetry axis (z-axis) by the spherical 
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co-ordinates (6, @). The eigenvalues of an axial spin-Hamiltonian for arbitrary 
spin with the external magnetic field at an arbitrary orientation with respect 
to the symmetry axis can be expressed in a closed form when the external 
magnetic field is large compared to the splittings due to the crystalline electric 
field. The terms in these expressions, involving the spin S and the components 
of S,, can be tabulated for the commonest values of S facilitating the calculation 
of the eigenvalues. 
The axial spin-Hamiltonian has the form 


(1) M= 4 BB,S,+¢,8(B,S:+B,S,)+D[S:—45(S+1)] 


where the symbols are in accordance with convention. It is convenient, in 
strong magnetic fields, to choose a new set of axes such that the Zeeman 
splitting term is diagonalized. This is equivalent to defining a new term in g 
oriented at an angle y to the z-axis and defined by 


. cosy = (g/g) cos @ sin Y = (g,/g) sin @ 


when the external magnetic field B is oriented in the xz plane (@ = 0°) at an 
angle @ to the z-axis, and by 


cos Y = (g/g) cos 6 sin y = —(g,/g) sin 6 
when B is oriented in the yz plane (¢ = 7/2) at an angle 6 to the z-axis where 
g’ = g, cosO+g;, sin’. 
Right-handed rotation transforms (1) into 


2 2 
(2) w= £9BS,+D| St (#) cos 0+ Si («:) sin’@ 


— (214) : sin 6 cos 0(.S,S,+ SS.) — 4 S(S+1 | 


in the xz plane where (S,, S,) — (S,, S,) and 
© 2 9 © 2 9 
(3) WH wn ebs.+| S:(t!) cos'o+ si() sin’6 


a (2 i)’ sin 6 cos 0(S)Su+ SSe)—}% S(S+1 | 


in the yz plane where (S,,.S,) — (S,, S,). It is then found by second-order 
perturbation calculations that both spin-Hamiltonians yield 


t 


g 
4 


2 
p*. (#82) sin’a cos'o( 4). (45(54+1) 84" 1 


‘ =) 
at) otal 4) pec. o1s?__11% 
:) sin‘6| [7 )-[25(S+1)-2M*—1] 


) “\" . : 
(4) Ey = gaBM +5 E@) cos 6 | . (M*—4S(S+1)] 
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if G>D, where S, = M, M—1,..., —M+1, —M and G = gBBy = hy 
is the energy of the Zeeman splitting of the ground state in the absence of a 
crystalline electric field. 

The general expression for the energies of transitions between levels with 
strong field quantum number M and M—m, i.e., AM = +m is 


(5) Ey—Ey_-n = gsbm+P| a(t) cos'6— 1 lomc2ar—m 
2 
—D’. (2%) -sin’@ cos’0- (".) -[4S(S+1) -—24M(M—m) —8m’?—-1] 


f 
+p*-(t) sin'e-(%,) [25(S-+1) 6M (M—m) —2m*—1] 


which becomes for 
(a) Mo M-—1,i.e. AM = +1 


(6) Exy—Ey—1 = g8B+D(M—-1/2 )-[3(E8)" coset | 


an (4s) 
. 


F 4 
+p*-(t) -sin’ 6: 


(b) Mes M—2, i.e. AM = +2 


-sin’6 cos’ 6: (+) -[4S(S+1)—24M(M—1)-—9] 


| 
3h were a me 
8G [2S(S+1)—-6M(M—1)-—3]. 


on. 


(7) Ey—Exy ,; = 249-+20(a—1)-[ a(t) cos’ 6— | 
— p?.{ S82 ae, eee, : YRS ne =O 
D Fi ‘sin 6 cos 6: IG -[4S(S+ 1) —24M(M— 2) —33] 


° - 4 ! 2 
+p%-(&) sin'e-( 2.) .2s(s41) ~6M(M—2)—9].* 


(c) M+ M-3, i.e. AM = 43 
(8) Ey—Eu—s = 3g8B+3D(M-3 »-[ a(#) cosé— | 
-p*. (taf) -sin’@ cos'e-( 2.)  tus(s-+1) —2uar(ar—3) —73 


A ‘ 
+p (&) sin's-(3.).p2s(s+1) -~6M(M-—3)—19]. 


_*These expressions are the corrected forms of equation (9.7) of Low (1960) and equation (5) 
of Bleaney (1951), 
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If the spin-Hamiltonian contains a term corresponding to the angular 
momentum operator O§ = (S$. ote a, which has the same transformation 
properties as the spherical harmonic Yeo, then, if S = 3, transitions corre- 
sponding to AM = +6 are allowed. 


mh | 


SS SS SS et 


|] NOS SS 


+8, —7 
+7, —6 
+6, —95d 
+5, —4 
+4, ae 
+3, —2 
+2, —1 
+1, 0 
M 
+15/2, —13/2 
+13/2, —11/2 
+11/2, —9/2 
+9/2, 7/2 
+7/2, —5/2 
+5/2, —3/2 
+3/2, —1/2 
+1/2, +1/2 
+8, —7 
% —6 
+6, at 
+5, —4 
+4, —3 
+3, —2 
+2, —1 
+1, 0 
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TABLE la 
[4S(S-+1)—24M(M— 


1)—9] 





s 
1/2 3/2 5/2 7/2 «9/2 11/2 13/2.—«15/2 
924 
672 —612 
~460 —408  —348 
288 —244 —192 —132 
156 —-120 —76 24 +436 
hk i 0 +44 496 +156 
-12 48 43 #472 +16 +168 4228 
0 +12 +32 +460 +96 +140 +192 +252 
1 2 3 4 5 6 7 8 
~ 1065 
~793  —729 
Jae | le tea 
~369 —321 —265 —201 
~217 -177 -129 —-73 —9 
i a =a -2ak eee. eee 
-33 —-9 +423 +463 +111 +167 +4231 
-1 415 +39 +71 +111 +159 +4215 +279 
TABLE Ib 
[25(S+1) -6M(M—1)—3] 
Ss : 
1/2 3/2 5/2 72 9/2 11/2 13/2.—«15/2 
~168 
— 120 —90 
—80 sm Aan 
—48 —26 0 +36 
a +16 +42 +72 
J: ae aa 446 +72 +4102 
6 6 660 la ea 
0 46 +16 +390 44 ##+70 496 « +128 
2 3 4 5 6 7 8 
~195 
~143 111 
—99 
~63 —39 
2 oe ee 
ae ~ 421 445 
3 a aa ee ee 
+t ° 4 [ap See ae ae 








| 





an 
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(9) Exy—Ey—s = 6g8B+6D(M— 3)| 3 (£1) cos's— 1] 


—D*° (ti s) sin’@ cos'd( ©) uss 1) —24M(M—6) —289] 


8G 


+D° (é :) sin ( &)iescs-1 )—-6M(M-—6) —73]. 


TABLE Ila 
Ree 24 nm M—2)—33] 
































S 
M 1/2 3/2 5/2 7/2 9/2 11/2 13/2 15/2 
+15/2, —11/2 —768 
+13/2, —9/2 540 —480 
+11/2, —7/2 —352 —300 —240 
+9/2, —5/2 — 204 — 160 —108  —48 
+7/2, —3/2 —96 —60 —16 +36 +96 
+5/2, —1/2 —28 0 +36 +80 +132 +192 
+3/2, +1/2 0 +20 +48 +84 +128 +180 +240 
1 2 3 4 5 6 7 8 
+8, —6 —897 
+i, —5 —649 —585 
+6, —4 —441 —385 —321 
+5, —3 — 273 — 225 —169 —105 
+4, —2 —-145 —105 —57 =m +63 
+3, —] —57 —25 +15 +63 +119 +183 
2, 0 —9 +15 +47 +87 +135 +191 +255 
+1, +1 —1 +15 +39 +71 +111 +159 215 +279 
TABLE IIb 
eerers _— M —2)-—9] 
s 
M 1/2 3/2 5/2 7/2 9/2 11/2 13/2 15/2 
+15/2, —11/2 —129 
+13/2, —9/2 —87 —57 
+11/2, ~¢/2 —53 —27 +3 
+9/2, —5/2 —27 —5 +21 +51 
+7/2, —3/2 —=¥ +9 +31 +57 +87 
+5/2, —1/2 +1 +15 +33 +55 +81 +111 
+3/2, +1/2 +3 +13 +27 +45 +67 +93 +123 
aa 1 2 ' 3 ae ; 5 6 “ha a ae 
+8, —6 — 153 
“ies —5 —107 —7% 
+6, —4 —69 —41 —9 
+5, —3 —39 —15 +13 +45 
+4, —2 —17 +3 +27 +55 +87 
+3, —1 —3 +13 +33 +57 +85 +117 
+2, 0 +3 +15 +31 +51 +75 +103 +135 
+1, +1 +1 +9 +21 +37 +57 +81 +109 +141 
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Equations (7), (8), and (9) are useful expressions in the approximate 
analysis of the data for the angular spectrum of the AM = +2, +3, and +6 
transitions observed in dilute gadolinium ethyl sulphate by Buckmaster 
(1956a, 6). No attempt has been made to find expressions in closed form for 
the transition probabilities when the magnetic field is arbitrarily oriented 
with respect to the crystal axis. In general, the intensity of the AM = +2 
and +8 transitions are weaker than that of the AM = +1 transitions by a 
factor of (D/B)? and (D/B)‘ respectively and vanish only if the external 
magnetic field is parallel to the crystal axis. 

Tables I to IV give the values of the terms involving only S and M that 
appear in equations (6), (7), (8), and (9) for the commonest values of S. 


TABLE IIa 
[4.5(S+1)—24 











M 1/2 3/2 5/2 7/2 9/2 11/2 13/2 15/2 
+15/2, —9/2 —628 
+13/2, —7/2 —424 —364 
+11/2, —5/2 —260 -208 —148 
+09/2, —3/2 -196 —92 —40 +20 
47/2, 1/2 —52 —16 +28 +80 +140 
+5/2, +1/2 —8 +20 +56 +100 +152 +212 
43/2, +3/2 —4 +16 +44 +480 +124 +176 +236 

1 2 3 4 5 6 7 8 
+8, —5 —745 
+7, —4 —521 —457 
+6, =f —337 -21 —2i7 
+5, —32 -19 -145 -89 —25 
+4, | -89 —49 af +55 +119 
+3, 0 —28 +7 +47 +95 +151 +215 
> +1 -1 +2 +55 +95 +143 +199 +263 
TABLE IIIb 
[2S(S+1)—-6M(M—3)—19] 

M 1/2 3/2 5/2 7/2 9/2 11/2 13/2 15/2 
+15/2, —9/2 —94 
4138/2, ~—7/2 —§8 —2 
+11/2, —5/2 —30 4 +26 
40/2, -3/2 ~10 +12 +38 +68 
47/2, —1/2 +2 +20 +42 +68 +98 
+5/2, +1/2 +6 +20 +38 +60 +86 +116 
+3/2, +3/2 +2 +12 426 +44 +66 +92 +4122 

1 2 3 4 5 6 7 8 
+8, —5 —115 
+7, —4 -7%5 —43 
+6, -3 —43 —15 +17 
+5, =% —19 +5 +33 +65 
+4, = <3 +17 +41 +69 +4101 


3, 0 +5 © +21 +41 +65 +93 +125 
+2, = ai! +5 +17 +33 +53 +77 +105 +137 





; 
; 
; 


| 
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TABLE IVa 
[4.5(S+1) —24.M(M—6) —289] 
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S 
M is «(3A 5/2 7/2 9/2 1/2 13/2 15/2 
+15/2, —3/2 —304 
+13/2, —1/2 -172 —112 
+11/2, +1/2 —80 —28 +32 
+9/2, +3/2 —28 +16 +68 +128 
+7/2, +5/2 —16 20 +64 +116 +176 
2 3 4 5 6 7 8 
+8, -2 —385 
+7, <1 —233 —169 
+6, 0 -121 -6 -1l 
+5, +1 —49 -1 +55 +119 
+4, +42 -17 +23 +71 +127 +191 
+3, +3 —2% +7 +47 +95 +151 +4215 
TABLE IVb 
[2S(S+1)-6M(M—6)—73] 
S 

M 1/2 3/2 5/2 7/2 9/2 11/2 13/2 15/2 
+15/2, —3/2 —13 
+13/2, —1/2 +5 +35 
+11/2, +1/2 +15 +41 +71 
+9/2, +3/2 +17 +39 +65 +95 
+7/2, +5/2 +11 +29 +51 +77 +107 
1 2 3 4 5 6 7 8 

+8, 2 —25 
+7, -l -3 +29 
+6, 0 +11 +39 «+71 
+5, +1 +17 +41 +69 +101 
+4, +2 +15 +35 +59 +87 +119 
+3, +3 +5 +21 +441 +65 +93 +125 
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THE 1.2-MEV BETA DECAY OF [8 
K. G. McNEILL AND J. D. PRENTICE 


The decay scheme of I|'8 has been extensively studied by both beta- and 
gamma-ray spectroscopy (for example, Benczer et al. 1956; Gupta and Jha 
1956; Wapstra et al. 1953). The Kurie plot of the beta spectrum obtained by 
Benczer et al. (1956) showed the presence of three components. By subtraction, 
the end-point energies and relative intensities of the two highest-energy groups 
were established and the presence of a weak lower-energy group (83) was 
indicated. Owing to the inherent inaccuracy of the successive subtractions 
neither the end point nor the-shape of the Kurie plot for 83; could be obtained. 
By subtracting from the end-point energy of the ground-state beta spectrum 
(Es, = 2.12 Mev) the measured energy (0.98 Mev) of the gamma ray from 
the second excited state of Xe’, they derived a value of 1.14 Mev for the 
maximum energy of 83. 

By observing the beta spectrum in coincidence with the 0.98-Mev gamma 
ray the present authors have confirmed the proposed decay scheme and have 
measured the end-point energy and spectrum shape of 8; directly. 

A Nal (TI) crystal, irradiated with neutrons, was used as both source and 
detector for the beta particles thus providing a zero thickness source and 
eliminating self-absorption, counter window, and scattering corrections. A 
1 in. diameter by 1 in. long Nal crystal was irradiated for 40 minutes in a 
water tank surrounding a Sb—Be neutron source. It was partially shielded 
from the gamma rays from the source by a 1-in. thickness of lead which 
served to reduce sufficiently the long term phosphorescence which otherwise 
appeared. After irradiation the crystal was mounted on a DuMont 6292 
photomultiplier connected by conventional electronics to a 100-channel pulse- 
height analyzer (100 PHA). A 12 in. diameter by 2 in. long Nal (TI) crystal 
was placed near the irradiated one to detect the emergent gamma rays from 
the latter. The output from the second detector was analyzed by a single- 
channel analyzer set to accept pulses from the photopeak of the 0.98-Mev 
gamma ray emitted by the second excited state of the daughter nucleus, 
Xe!**, The output from this analyzer gated the 100 PHA so that the latter 
only accepted pulses, in the irradiated crystal, in coincidence with 0.98-Mev 
gamma rays. As separate gamma-gamma coincidence experiments with an 
[8 source showed that no other gamma rays were in coincidence with those 
at 0.98 Mev, the spectrum obtained from the 100 PHA, after suitable cor- 
rection, was the spectrum of the beta particles emitted in the transition 
from the ground state of I'8 to the second excited state of Xe!*8. Corrections 
were applied for random coincidences between the 0.98-Mev gamma rays and 
beta particles from the higher-energy branches of the beta decay and also 
for true coincidences due to the decay of Na*‘ produced by the irradiation. 
Pulses from the Compton distribution of the higher-energy gamma rays 
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emitted by Mg*™ following the beta decay of Na*™ could appear in the 0.98- 
Mev gate and cause spurious real coincidences. The Na** has, however, a 
half life of 15 hours and for this reason and because the capture cross section 
of Na* for neutrons is small, the Na* activity was weak. The long half life 
also permitted an experimental correction to be made for the Na* contribution 
by observing its spectrum after the I) had decayed away. 

The Kurie plot of the spectrum finally obtained is shown in Fig. 1. It is 
seen to have a straight-line shape and to intercept the energy axis at 1.20 
Mev. The estimated error in this result is +0.05 Mev and is largely caused 
by statistical errors in the counts obtained. 





® ENERGY — Mev 





® ENERGY— Mev 


Fic. 1. The beta spectrum of I in coincidence with 0.98-Mev gamma rays. N = number 
of beta particles per channel. G = modified Fermi function = (p/W)F. W = total energy 


of beta particles in units of m c?. 
Fic. 2. The beta spectrum in coincidence with 0.54-Mev gamma rays. 


A similar experiment with the gamma-ray gate set on the 0.54-Mev cascade 
transition from the second to the first excited states of Xe!’ yielded a straight- 
line shape and a value of 1.13+0.04 Mev (Fig. 2). We, therefore, conclude 
that the best value for the end-point energy of 8; is 1.16+0.05 Mev and 
that the spectrum shape is consistent with its being an allowed transition. 
This value agrees with the value of 1.14+0.02 Mev deduced by Benczer et 
al. (1956) and together with their result that 2% of all transitions proceed 
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by this branch leads to a log ft value of 6.6. The shape and log ft value suggest 
an allowed transition and tend to confirm the 1+ and 2+ assignments for the 
ground state of I’ and second excited state of Xe!” respectively. 

In connection with the usefulness of the experimental technique of making 
the source in the crystal it is perhaps of interest to note that it was also 
possible to confirm that a 0.75-Mev y-ray observed in the gamma-ray spectrum 
belonged to a competing electron capture decay branch. In this case the Te 
K X rays following electron capture, absorbed in the radioactive crystal, were 
used to gate the spectrum of gamma rays detected in an external crystal. A 
gate set to accept the K X rays inevitably accepted many pulses from low- 
energy §-particles in coincidence with all the other gamma rays as well 
but when the gate was moved to either side of the X-ray peak the 0.75-Mev 
gamma ray disappeared, the remainder of the coincidence spectrum remaining 
unchanged. Thus the 0.75-Mev y-ray was clearly shown to be in coincidence 
with Te K X rays but not with @-particles. 


3ENCZER, N., FAKRELLY, B., Koerts, L., and Wu, C. S. 1956. Phys. Rev. 101, 1027. 
Gupta, R. K. and Jua, S. 1956. Nuclear Phys. 1, 2. 
Wapstra, A. H., VERSTER, N. F., and BOELHOUWER, M. 1953. Physica, 19, 138. 
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INELASTIC SCATTERING OF NEUTRONS FROM 
ROTATING LIQUID HELIUM II* 


A. D. B. Woops 


This note describes an attempt to see if macroscopic rotation of liquid helium 
I] changes the character of the energy/momentum dispersion curve for the 
Landau excitations; no effect was observed. The excitations have been ex- 
tensively studied in recent years using inelastic neutron scattering (Yarnell, 
Arnold, Bendt, and Kerr 1959; Palevsky, Otnes, and Larsson 1958; Henshaw 
1958; and Henshaw and Woods 1960, 1961). The interesting properties of 
liquid helium II! moving at high velocities (see, for example, Walmsley and 
Lane 1958 and Dowley and Hallett 1960) suggested that it might be desirable 
to determine whether the excitations, particularly those in the region of the 
roton minimum, are affected by rotation of the helium. In this experiment 
changes in either the energy of the scattered neutron group corresponding to 
an excitation near the roton minimum or in its width were looked for, An 
increase in width of the neutron group on rotation would imply enhanced 
roton/roton interaction and, hence, shorter roton lifetimes when the helium 
was rotating. 

*Issued as A.E.C.L. No. 1262 
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Accordingly the following experiment was performed. Neutrons with an 
initial wavelength of 4.06A were scattered from liquid helium through an 
angle of 814° and the scattered neutrons were energy-analyzed using a time- 
of-flight rotating crystal spectrometer (Brockhouse 1960) at the Chalk River 
NRU reactor. The helium was in a chamber with an inside diameter of 2.6 in. 
which contained a bucket 2} in. in diameter. This bucket was rotated at a 
speed of 60 r.p.m. with the helium at a temperature of 1.54° K. At this tem- 
perature the roton line width is known to be perceptibly broadened (Larsson 
and Otnes 1959; Henshaw and Woods 1960), so that further increase on 
rotation should have been detectable. The runs lasted for almost 1 hour and 
rotation was begun about 15 minutes before counting. 

The results obtained are shown in Fig. 1 where the counting rate per time 
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Fic. 1. Neutron groups observed for stationary and rotating conditions of inner bucket. 
The broken vertical line represents the position of zero energy transfer. The position of the 
peak corresponds to an excitation in the region of the roton minimum. 


channel (with background scattering subtracted) is plotted for both the 
rotating and stationary conditions of the inner bucket. The scales have been 
adjusted so that the peak heights are equal. No significant difference between 
the patterns is observable. It is estimated that any shift in peak position is 
less than 0.3° K and any increase in width less than 0.5° K. 

The author would like to thank Dr. D. G. Henshaw tor the loan of the 
cryostat. Mr. E. A. Glaser and Mr. J. R. Freeborn gave valuable technical 


assistance. 
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THE EFFECT OF DISSOLVED IRON ON THE THERMOELECTRICITY 
OF SILVER AT VERY LOW TEMPERATURES* 


W. B. PEARSON AND I. M. TEMPLETON 


Thermoelectric measurements made in this laboratory, as well as earlier 
measurements of Borelius et al. (1932) and Keesom and Matthijs (1935), 
show that iron causes a very large characteristic thermoelectric power at low 
temperatures (say <10° K) when dissolved in very dilute concentration in 
Cu or Au. These ‘giant’? thermoelectric powers may be as much as 1000 
times larger than values predicted for diffusion thermoelectricity by the 
original Wilson theory (1936). No comparable measurements appear to have 
been made for Ag—Fe alloys, and indeed Hansen (1958) indicates that the 
solid solubility of Fe in Ag is very small. We have therefore prepared an 
Ag-Fe alloy wire which we quenched rapidly from 940° C by ejecting it for- 
cibly from the annealing furnace into a bath of iced water. The wire had a 
residual resistance ratio of 75.810-* indicating that some iron was indeed 
retained in solid solution by the quench. The absolute thermoelectric power 
of this alloy (Fig. 1) shows that iron also gives rise to a very large absolute 
thermoelectric power at low temperatures when dissolved in silver; further- 
more, the alloy had a minimum in its electrical resistivity. 

The low solid solubility of Fe in Ag, compared to that of Cu or Au, perhaps 
explains why it is relatively easy to obtain pure silver with a positive thermo- 
electric power at low temperatures. Both Pullan (1953) and Pearson (1960) 
report measurements on “‘pure’’ Ag whose positive absolute thermoelectric 


*Issued as N.R.C. No. 6386. 
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TEMP. °K 


Ag 42x10 


Ag-Fe 
75.8 x 107% 





Fic. 1. Absolute thermoelectric power of Ag-Fe, Cu-Fe, and Au-Fe alloys and “‘pure’’ 
Tadanac Ag as a function of temperature. The figures give the residual resistance ratios, 
Ry 2° k/( Ree? k —Ri.2° k), for the alloys. 


power increases with slightly less than the square of the absolute temperature 
in the temperature range from 2° to about 18° K. The sample of silver of 
Johnson, Matthey and Company reported to have, at very low temperatures, 
a small negative almost constant S (MacDonald et al. 1958) probably con- 
tained a trace of dissolved Fe, since the same silver also showed a minimum 


in its electrical resistivity (Pearson 1959), 


We thank Mr. J. W. Fisher for preparing and heat-treating the Ag—Fe wire. 
Boretius, G., KEEsom, W. H., Jonansson, C. H., and Linpg, J.O. 1932. Proc. Acad. Sci. 


Amsterdam, 35 (1), 15, 25. 
HANSEN, M. 1958. Constitution of binary alloys (McGraw-Hill Book Co, Inc.), p. 20. 
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THE K-LL AUGER SPECTRUM OF ;;Cs'%3 * 
R. L. GRAHAM, F. Brown, G. T. EWAN, AND J. UHLERT 


The measurements of the K—LL Auger lines reported here were made 
primarily to test the non-relativistic theory of Asaad and Burhop (1958) at 
the lowest value of the atomic number, Z, that has so far been studied with 
the Chalk River iron-free 8-ray spectrometer (Graham et al. 1960). Discre- 
pancies exist between the measured intensity ratios of the K—LL lines and 
the non-relativistic predictions but become smaller as Z decreases, as might 
be expected (Graham and Merritt 1961). Also the theoretical treatment 
suggests that there should be three satellite lines in addition to the six main 
lines of the K—LL Auger spectrum as a consequence of the nature of the 
coupling between the two final electronic Z state vacancies. This coupling is 
expected to be close to a pure j—7 coupling at high Z (six-line spectrum) and 
become intermediate between j—j7 and L-—S coupling at lower Z (nine-line 
spectrum). The predictions suggest that the satellites should become more 
prominent and better separated from their partners at lower Z and hence 
more easily detected experimentally. In our experiments to date, at Z values 
>62, it has been possible to detect in each case only one of the three predicted 
satellite lines. In this note we report measurements of the energies and 
intensities of the K—LL Auger lines at Z = 55. We again observe the K-L,L; 
(®P.2) satellite line and report preliminary experimental evidence for the other 
two satellite lines. 

The source material, Xe™’, was prepared by neutron irradiation of natural 
Xe for 6 days in a neutron flux of 2.5610" cm-? sec—! in the NRU reactor. 
We are indebted to Dr. R. O. Kelly for supplying an iron capsule and for 
arranging the irradiation. An isotopically pure source of Xe* was prepared 
using the newly installed mass separator (Brown et al. 1961). The Xe! 
activity was confined to an area 20 mmX1 mm by means of a defining slit 
through which passed about half the ions of mass 133. The ions were driven 
_about 8 yg/cm? into an 800 ug/cm? aluminum-backing foil by virtue of their 
30-kev kinetic energy (Davies 1961). 


*Issued as A.E.C.L. No. 1264. 
{Visitor from the Nobel Institute of Physics, Stockholm, Sweden. 
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For the measurements reported here, the spectrometer was adjusted to have 
an instrumental resolution of ~0.07% in momentum. Four sweeps of the 
spectrum were made using 100-second counting intervals. The net counting 
rate after subtraction of the y background rate measured by a separate counter 
is shown in Fig. 1. The error bars indicate standard deviations of the gross 
counts at each current setting. The experimental line width of ~0.15% is 
mainly attributable to the straggling of electrons as they emerge through the 
aluminum foil from the penetration depth of the Xe atoms. 

The energies and relative intensities are summarized in Table I. The relative 


TABLE | 


Energies and relative intensities of K-LL Auger lines in 5;Cs' 


Line energy (kev) Relative intensity 





Transition Theory? Experiment Theory? Experiment 
L,L,('So) 24.34 24.39+0.02 1.00 1.00+0.06 
L,L('P) 24.70 24.75+0.02 1.44 | a = 
LiL? Po) 24.73 0.05 Lae 
LoL 2('S¢) 25.01 Or < 0.20 5 \ 
LLi(2P;) aoa 25 .08+0.02 20 | 1.58+0.19 | 1 9949.07 
L,L3(3P2) 25.08 25.14+0.03 0.24 0.42+0.19 } 
ene 25.38 25 .43+0.02 4.56 4.23+0.07 
3L3(3 Po) 25.67 0.22 . 
L3L3(*P») 25.72 25.79-+0.02 1.91 } 1.91 £0.06 


/ 


@Asaad and Burhop (1958). 


intensities were obtained from the line areas. For purposes of matching line 
shapes the net counting rates, after subtraction of the 6-continuum rates, 
were plotted on a semilog scale in the manner discussed in earlier papers (see, 
for example, Graham et al. 1960). The K-ZL2L;(!D2) line was used as the 
“‘standard”’ shape since it is most intense and is free from satellites. The shape 
of the intense K 80.99 internal conversion line shown at the right in Fig. 1 
was also used as a guide in estimating the shape of the low-energy tails for 
each of the Auger lines which is not well determined statistically. The arrows 
in Fig. 1 show the positions of the six Auger lines and their three satellites as 
predicted by the formulae of Asaad and Burhop shifted ~50 ev to match 
the K-L,L; and K-—L;L; peaks. The discrepancies between theoretical and 
experimental line energies listed in Table I are similar to those found at other 
Z values, see, for example, Graham and Merritt (1961). The central Auger 
line has three components, the K—LeL.('S9) and K-L,L,(°P1) being indis- 
tinguishable because of the close energy spacing. The high-energy shoulder, 
however, is caused by the K—L,L;(*P») satellite line which has been seen at 
other Z values. The reader is referred to the paper by Graham and Merritt 
(1961) in this issue for a discussion of systematic trends in relative intensities 
with Z. 

Note that the K—L,L2 and K-L;L; lines are not consistent with the standard 
line shape and have excess counts at the locations predicted for their respective 
satellites. The K—L,L;, line, on the other hand, is consistent within errors. We 
believe the Lil: and L3L; discrepancies are too large to be attributable to 
experimental error and that they are real. The Z = 54 Auger lines due to the 
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NOTES 1089 


K 233.2 line in the Xe 2.3-day isomer were [1% of the Z = 55 Auger lines 
at the time of our experiments and hence were too weak to account for these 
excess counts. Thus there is qualitative evidence for the additional two satellite 
lines expected theoretically. The use of a stronger source, prepared with a 
de-accelerating system to slow down the Xe ions and reduce their penetration, 
should make it possible to investigate these effects quantitatively. 

The absolute coefficient for K-—LL Auger yield can be deduced from the 
data of Fig. 1 since the K-shell vacancies arise almost entirely (>99%) from 
the K 80.99 conversion line shown at the right. This gives the ratio 
2(K-LL)/K 80.99 = (6.58+0.10)%. From an additional measurement of the 
K-LM and K-LN Auger lines we deduce that 2(K-LY)/2(K-LL) = 0.48 
+0.03. Using a value of 2(K—X Y)/Z(K-LL) ~ 0.07 (Wapstra et al. 1959), 
one obtains a value for the total Auger yield of (10.2+0.2)%. The resulting 
fluorescent yield of wc = 0.898 is to be compared with the value of wx = 
0.877 predicted by the semiempirical formula of Hagedoorn and Wapstra 
(Wapstra et al. 1959) for Z = 55. 


AsaaD, W. N. and Buruop, E. H. S. 1958.. Proc. Phys. Soc. 71, 360. 
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Cone. R. L., Ewan, G. T., and GeiGcer, J. S. 1960. Nuclear Instr. and Methods, 9, 
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LETTERS TO THE EDITOR 


Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


The Lifetime of the 6.25-kev Level in Ta‘*!* 


Recently Muir and Boehm (1961) have shown the existence of a 6.25-kev level in Ta!*!, 
and have determined from the M-subshell conversion ratios that it decays by El gamma- 
quantum emission about 2% of the time. 

The Weisskopf estimate for the mean life of a 6.25-kev electric dipole transition in Ta is 
10~* sec. Making a correction for the known conversion coefficient (44+7) and applying an 
inhibition factor of 10° which is typical for low-energy El transitions in heavy elements 
(Wilkinson 1960), we obtain 2.5 usec as an estimate for the mean life of the level. 

We have observed the transition and have measured the level lifetime. 

A source of W'*! (140 days) was prepared by employing 10-Mev protons from the Chalk 
River Tandem accelerator to initiate the (p,m) reaction in a 0.0005-in. tantalum foil. The 
tungsten was then separated chemically! and deposited on a thin teflon backing. 

The low-energy gamma-ray spectrum was observed in a krypton- and CO>-filled proportional 
counter having a thin Be window. This spectrum (Fig. 1A) shows the LZ X rays of tantalum 
and the 6.25-kev gamma ray. Figure 1A also shows a spectrum taken in delayed coincidence 
with the Ta K X rays (57.8 kev), the latter being detected in a Nal crystal, 1.5 in. diameter, 
2 mm thick. The gate on the proportional-counter spectrum was opened 3.5 usec after the 
detection of a K X-ray pulse and held open for 10 usec. This coincidence spectrum shows clearly 
the 6.25-kev transition, as well as random coincidences with the L X rays. The area under the 
6.25-kev coincidence peak has been measured for gate openings delayed by 3.5, 10, 15, and 
140 usec, the latter serving to provide a measurement of the random spectrum. The net counting 
rate in the remaining three delayed-coincidence spectra shows an exponential decrease with 
time corresponding to a half-life of 6.10.7 X10~* sec. This result is shown in Fig. 1B. 

Another measurement of the half-life has been made employing a time sorter originally set 
up to measure the lifetime of the 14.4-kev level in Fe” (Clark 1961). The same detectors were 
employed, start pulses for the timing being derived from K X-ray pulses from the 2-mm 
Nal crystal, and stop pulses from the 6.25-kev peak in the proportional-counter spectrum. The 
result is shown in Fig. 1C, and gives a value 74/2 = 7.0+0.4X10~ sec. 

Our average value is thus 71/2 = 6.8+0.4X10~ sec, which implies a mean life 9.80.5 X 10~* 
sec, in reasonable agreement with the estimate made for an E1 transition. 

The value obtained for the lifetime implies an experimental width of 1.6X10~ cm/sec in 
a Méssbauer-effect experiment. The gamma-ray energy, combined with the Debye temperature 
of tantalum (~230° K) implies that the fraction of recoilless transitions at room temperature 
will be as high as 90%. Two of the main problems in the investigation of the recoilless radiation 
from the Ta!*! 6.25-kev level will be those of extraneous vibration and temperature variations. 


The assistance of Mr. G. C. Dixon is gratefully acknowledged. 


CrarK, M. A. 1961. To be published. 
Mutr, A. H. and Boru, F. 1961. Phys. Rev. To be published. 
Wirkinson, D. H. 1960. Nuclear spectroscopy, B, 873, edited by Fay Ajzenberg-Selove (Academic Press). 
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Fic. 1. A. Direct spectrum of low-energy gamma rays in the decay of W!8!, The spectrum in delayed coincidence 
with Ta K X rays is also shown. 

B. Delayed-coincidence intensity of the 6.25-kev gamma ray in Ta!®! (10 ysec gate). 

C. Delayed-coincidence intensity of the 6.25-kev gamma ray using the time sorter (1 wsec per channel). 


ANNOUNCEMENTS 


Frequency Measurement of Standard Frequency Transmissions'’? 


Measurements are made at Ottawa, Canada, using N.R.C. caesium-beam frequency 
resonator as reference standard (with an assumed frequency of 9 192 631 770 c.p.s.). Frequency 
deviations from nominal are quoted in parts per 10. A negative sign indicates that the 
frequency is below nominal. 





GBR, 16 kc/s 


MSF, 63-hour 24-hour 
60 kc/s average* averaget 





—140 —150 — 146 
— 152 —152 — 148 
—157 —152 —148 
—156 —155 —149 
—159 —157 — 156 
— 152 — 153 —149 
N.M. —147 —149 
N.M. — 139 —141 
N.M. —147 —138 
—154 — 152 —143 
—161 — 146 —147 
— 154 N.M. —145 
— 160 — 149 —149 
—159 — 152 — 152 
—148 —144 — 147 
—152 — 143 — 142 
— 156 —147 —142 
—152 —146 — 143 
—150 —151 —142 
— 156 —149 —140 
N.M. — 156 —151 
— 156 —157 — 152 
—155 —154 — 156 
—153 —147 —149 
N.M. — 138 — 143 
— 153 — 150 — 148 
—153 —149 —145 
—157 —152 —147 
— 138 — 143 —144 
— 136 —141 —142 


Average — 153 —149 — 146 


Midmonthly 
mean of WWV —147 


Note: N.M. no measurement. 
*Time of observations: 11:30 to 18:00 hours U.T. 
tTime of observations: 15 to 15 hours U.T. 


RECEIVED May 12, 1961. 
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MANUSCRIPTS 


General.—Manuscripts, in English or French, should be typewritten, double spaced, on 
paper 83X11 in. The original and one copy are to be submitted. Tables and captions 
for the figures should be placed at the end of the manuscript. Every sheet of the manuscript 
should be numbered. Style, arrangement, spelling, and abbreviations should conform to the 
usage of recent numbers of this journal. Greek letters or unusual signs should be written plainly 
or explained by marginal notes. Characters to be set in boldface type should be indicated 
by a wavy line below each character. Superscripts and subscripts must be legible and carefully 
placed. Manuscripts and illustrations should be carefully checked before they are submitted. 
Authors will be charged for unnecessary deviations from the usual format and for changes 
made in the proof that are considered excessive or unnecessary. 

Abstract.—An abstract of not more than about 200 words, indicating the scope of the work 
and the principal findings, is required, except in Notes. 

References.—References should be listed alphabetically by authors’ names, un- 
numbered, and typed after the text. The form of the citations should be that used in current 
issues of this journal; in references to papers in periodicals, titles should not be given and 
only initial page numbers are required. The names of periodicals should be abbreviated in 
the form given in the most recent List of Periodicals Abstracted by Chemical Abstracts. All 
citations should be checked with the original articles and each one referred to in the text by 
the authors’ names and the year. 

Tables.—Tables should be numbered in roman numerals and each table referred to in the 
text. Titles should always be given but should be brief; column headings should be brief and 
descriptive matter in the tables confined to a minimum. Vertical rules should not be used. 
Numerous small tables should be avoided. 


ILLUSTRATIONS 


General.—All figures (including each figure of the plates) should be numbered con- 
secutively from 1 up, in arabic numerals, and each figure referred to in the text. The author’s 
name, title of the paper, and figure number should be written in the lower left corner of the 
sheets on which the illustrations appear. Captions should not be written on the illustrations. 

Line drawings.—Drawings should be carefully made with India ink on white drawing 
paper, blue tracing linen, or co-ordinate paper ruled in blue only; any co-ordinate lines that 
are to appear in the reproduction should be ruled in black ink. Paper ruled in green, yellow, 
or red should not be used. All lines must be of sufficient thickness to reproduce well. Decimal 
points, periods, and stippled dots must be solid black circles large enough to be reduced if 
necessary. Letters and numerals should be neatly made, preferably with a stencil (do NOT 
use typewriting) and be of such size that the smallest lettering will be not less than 1 mm 
high when the figure is reduced to a suitable size. Many drawings are made too large; originals 
should not be more than 2 or 3 times the size of the desired reproduction. Whenever possible 
two or more drawings should be grouped to reduce the number of cuts required. In such groups 
of drawings, or in large drawings, full use of the space available should be made; the ratio of 
height to width should conform to that of a journal page (4}X7}in.), but allowance must 
be made for the captions. The original drawings and one set of clear copies (e.g. small 
photographs) are to be submitted. 

Photographs.—Prints should be made on glossy paper, with strong contrasts. They 
should be trimmed so that essential features only are shown and mounted carefully, with 
rubber cement, on white cardboard, with no space between those arranged in groups. In 
mounting, full use of the space available should be made. Photographs are to be sub- 
mitted in duplicate; if they are to be reproduced in groups one set should be mounted, 
the duplicate set unmounted. 
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A total of 100 reprints of each paper, without covers, are supplied free. Additional reprints, 
with or without covers, may be purchased at the time of publication. 

Charges for reprints are based on the number of printed pages, which may be calculated 
approximately by multiplying by 0.6 the number of manuscript pages (double-spaced type- 
written sheets, 84 X1lin.) and including the space occupied by illustrations. Prices and 
instructions for ordering reprints are sent out with the galley proof. 
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